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Abstract 


A  number  of  definitions  of  zeros  of  linear  time -invariant  multi- 
variable  systems  have  appeared  recently.  This  work  surveys  selected  litera¬ 
ture  on  these  zeros.  Two  questions  are  addressed  here.  First,  how  are  zeros 
defined  and  how  are  these  definitions  interrelated.  Second,  how  can  they  be 
calculated . 

NThe  definitions  of  zeros  are  considered  for  three  system  representa¬ 
tions:  1)  the  transfer  function  matrix,  2)  the  state  space  representation 
in  the  frequency  domain,  and  3)  the  state  space  representation  in  the  time 
domain.  The  definitions  of  zeros  for  transfer  function  matrices  are  shown  to 
be  (mostly)  equivalent.  However,  several  different  sets  of  zeros  are  defined 
for  state  space  representations.  The  interrelationships  between  all  of  these 
definitions  is  discussed  in  detail. 

It  turns  out  that  the  calculation  of  zeros  directly  from  the  defini¬ 
tions  is  not  always  tractable.  The  properties  of  zeros,  however,  provide 
several  algorithms  for  calculating  zeros.  These  properties  are  reviewed  and 
a  brief  summary  of  the  algorithms  for  calculating  zeros  is  given. 

^Finally  ,  a  new  algorithm  for  the  calculation  of  invariant  zeros  is 
introduced.  It  is  based  on  the  geometrical  properties  of  linear  time  invariant 
systems.  This  algorithm  is  applicable  to  the  most  general  class  of  systems 
(A,B,C,D). 
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CHAPTER  1 
INTRODUCTION 

The  deftni.ci.on  of  a  zero  of  a  scalar  transfer  function  is  well 
known.  Indeed,  the  properties  of  zeros  are  very  important  for  describing 
the  open  and  closed  loop  behavior  of  dynamical  systems.  There  is  a  natural 
interest,  then,  in  extending  this  concept  to  linear  time- invariant  multi¬ 
variable  systems.  However,  it  is  not  clear  Just  how  this  should  be  done. 

The  approach  generally  taken  is  to  define  zeros  for  multivariable  systems 
so  chat  these  zeros  retain  some  property  of  zeros  of  a  scalar  transfer 
function.  As  it  turns  out,  the  zeros  so  defined  also  have  other  properties 
which  can  be  considered  generalizations  from  the  scalar  case. 

Zeros  defined  for  multivariable  systems  have  been  of  considerable 
interest  recently.  Macfarlane,  et  al,  have  developed  the  theory  extending 
the  classical  Nyquisc-Bode  and  root  locus  techniques  Co  linear  time -invariant 
multivariable  systems.  In  these  generalized  techniques,  Che  zeros  have  a 
role  analogous  to  their  role  in  the  classical  theory  for  single  input-single 
output  systems.  See  [1],  for  example.  Zeros  play  a  major  role  in  the 
construction  of  minimal  order  inverse  systems  [2],  Che  construction  of  re¬ 
duced  order  models  [3],  decoupling  theory  [4],  and  servomechanism  design 
[5],  [6].  More  recently,  they  have  been  used  in  relating  the  structure  and 
coefficients  of  the  quadratic  weighting  matrices  to  the  resulting  eigen- 
structure  of  the  optimal  state  regulator  [7]  and  in  the  stability  of  the 
optimal  state  regulator  using  high  gain  feedback  [8].  In  a  more  theoretical 
context,  zeros  have  proven  useful  in  describing  equivalence  classes  of 
linear  time- invariant  systems  under  the  action  of  the  group  of  state,  input, 
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and  output  space  transformations,  state  feedback,  and  output  injection  [9]. 
They  also  appear  in  such  diverse  applications  as  the  factorization  of  poly¬ 
nomial  matrices  [10],  and  many  other  areas. 

The  first  part  of  this  paper  is  devoted  to  a  survey  of  the  existing 
literature  on  zeros  for  Linear  time- invariant  multivariable  systems.  The 
representation  of  a  multivariable  system  has  three  forms:  the  transfer 
function  matrix,  state  space  system  in  the  time  domain,  and  state  space 
system  in  the  frequency  domain.  Correspondingly,  the  definitions  of  zeros 
has  been  extended  to  each  of  these  representations.  In  what  follows,  these 
definitions  are  given  and  their  interrelationship  is  explored.  This  in¬ 
cludes  Che  relationship  between  different  definitions  of  zeros  for  the 
same  system  representation  and  the  interrelationship  of  zeros  defined  for 
different  system  representations.  It  turns  out  chat  for  reachable  and 
observable  systems,  the  zeros  defined  from  each  representation  coincide.  It 
is  only  when  unreachable  or  unobservable  modes  occur  that  differences  in 
definitions  appear. 

Once  zeros  are  defined,  it  is  of  interest  to  develop  algorithms  for 
their  efficient  calculation.  As  is  turns  out,  neither  Che  definitions  nor 
their  elementary  properties  are  well  suited  to  either  hand  calculation  or 
calculation  on  a  digital  computer.  Therefore,  several  properties  of  zeros 
are  explored  which  lead  to  algorithms  for  computing  zeros.  Several  of  these 
properties  in  some  way  generalize  the  properties  of  zeros  of  scalar  transfer 
functions . 

Using  these  properties,  several  algorithms  for  calculating  zeros 


have  appeared  in  the  literature.  A  brief  review  of  the  algorithms  is  included 


quency  domain  properties  which  folLow  from  the  definitions  in  Chapcer  4. 

Also  introduced  are  geometrical  properties  of  zeros.  These  properties 
essentially  originate  in  the  time  domain  when  the  state  space  representation 
is  used.  Although  these  properties  are  also  a  consequence  of  the  definitions 
of  Chapter  4,  they  can  be  developed  independantly  by  a  geometrical  analysis 
of  system  properties. 

Chapcer  6  gives  a  brief  survey  of  the  algorithms  to  calculate  zeros 
which  have  appeared  in  Che  literature  Co  date.  A  brief  outline  of  the  theore¬ 
tical  basis  of  each  algorithm  is  given  along  with  a  discussion  of  the  type 
of  zeros  it  calculates  and  of  its  numerical  limitations. 

Chapter  7  presents  a  new  algorithm  for  calculating  zeros  based  on 
the  geometrical  properties  of  state  space  systems  given  in  Chapter  5.  The 
algorithm  actually  calculates  a  canonic  form  of  the  system  which  explicitly 
displays  fundamental  subspaces  closely  connected  to  zeros.  From  this  canonic 
form  the  zeros  are  easily  computed.  Several  examples  are  given  to  illustrate 


its  use. 
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CHAPTER  2 

MATHEMATICAL  PRELIMINARIES 


2.1.  Introduction 

Zeros  are  defined  for  both  continuous  and  discrete  systems.  In 
fact,  it  is  possible  to  discuss  zeros  for  both  types  of  systems  at  the 
same  time.  This  is  because  once  the  system  is  transformed  into  the  fre¬ 
quency  domain  (by  Laplace  or  z-transform)  it  assumes  an  algebraic  represen¬ 
tation.  In  this  form  the  analysis  techniques  are  the  same  for  both  systems. 
However,  it  is  pedagogically  convenient  to  exploit  this  connection  as  some 
properties  of  zeros  lend  themselves  nicely  to  continuous  time  interpretation 
while  others  arise  naturally  in  a  discrete  setting. 

In  this  chapter  terminology,  notation,  and  certain  preliminary  re¬ 
sults  are  introduced.  This  includes,  in  Section  2,  the  continuous  and  dis¬ 
crete  systems  to  be  discussed.  The  connection  between  these  system  repre¬ 
sentations  in  the  frequency  domain  is  discussed  and  certain  differences 
relevant  to  this  work  are  pointed  out.  Section  3  presents  a  brief  develop¬ 
ment  of  the  matrix  techniques  needed  to  analyze  these  systems,  namely  the 
Smith  form  and  the  Smith-McMillan  form.  These  forms  play  a  fundamental 
role  in  defining  zeros.  These  results  are  directly  applicable  to  some 
properties  of  zeros  and,  in  some  cases,  the  method  of  proof  will  be  used  in 
later  sections. 

The  notation  in  this  work  is  as  follows: 

The  field  of  real  numbers  will  be  denoted  R  and  the  field  of 
complex  numbers  by  C.  Some  of  the  results  below  can  be  generalized  to 
arbitrary  fields  of  real  and  complex  numbers ,  however,  this  work  will 
consider  only  the  fields  of  real  and  complex  numbers. 


b 


The  ring  of  polynomials  In  the  Indeterminant  a  is  the  set  of 
polynomials  in  s  with  finite  degree  and  coefficients  in  R.  Let  a(s)  be  a 
polynomial.  If  the  coefficient  of  the  term  with  largest  degree  in  a(s)  is  1, 
the  polynomial  is  said  to  be  monic .  Let  b(s)  be  another  polynomial.  The 
notation  "a(s)|b(s)"  means  "a(s)  divides  b(s)".  A  rational  function  is  the 
ratio  of  two  polynomials.  The  set  of  rational  functions  forms  a  field. 

Matrices  shall  be  denoted  by  capital  letters  and  their  elements  by 
small  letters.  Let  A  be  a  mxn  matrix  of  real  numbers.  The  notation 


diag[a1,a2> . . .  ,ap]  ,  p-mln(m.n) 
shall  mean  a  matrix  whose  elements  are  given  by 


/ 


i  "  J 


ij 


0  i  i  j 


i  ■  1,2 . m 

j  ■  1,2, .  . .  ,n 


1 

Note  that  it  is  not  assumed  that  m*n.  The  transpose  of  A  will  be  denoted 
T 

by  A  (not  A1).  The  rank  of  A  is  the  order  of  the  largest  minor  which  is 
not  identically  zero. 

A  matrix  whose  elements  are  polynomials  is  called  a  polynomial 
matrix .  Let  N(s)  be  a  polynomial  matrix.  The  local  rank  of  N(Q  for  any 
complex  number  z  is  the  rank  of  N(z)  (each  element  evaluated  at  s  ■  z)  and  is 
denoted  p [ N (z ) ] .  The  normal  rank  of  N(')  is  defined  as 


max  p  [  N  (z )  ]  ‘1  P  [  N ( • )  ] 


Alternatively,  the  normal  rank  of  N(s)  can  be  determined  by  considering  the 
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elements  .is  members  of  the  ring  of  polynomials  In  s  and  using  the  definition 
of  the  rank  of  a  matrix.  Suppose  that  N(s)  Is  square.  Then  N(s)  la  said 
to  be  unlmodu lar  if  Its  determinant  la  a  non-zero  constant.  This  Implies 
that  the  matrix  Is  invertible,  l.e.  n'^s)  is  again  polynomial  matrix. 

A  matrix  M  la  said  to  be  a  common  left  divisor  of  the  matrices  N 
and  P  if  there  exists  matrices  N  and  D  such  that 

N  -  MN,  D  -  MD 

The  matrices N and  D  are  said  to  be  right  multiplies  of  M.  A  matrix  L  Is  said 
to  be  a  greatest  common  left  divisor  of  N  and  D  if  It  la  a  common  left  divisor 
of  N  audD  and  L  Is  a  right  multiple  of  every  other  common  laft  diviaor  of  N 
and  D.  Suppose  that  these  matrices  are  polynomial  matrices.  If  the  greatest 
common  left  divisor  L  of  N  and  D  is  unlmodu lar,  then  N  and  D  are  said  to  be 
left  coprime.  Right  coprime  matrices  are  defined  similarly. 

For  vector  space  analysts,  Roman  letters  will  denote  matrices  or 
the  maps  which  they  represent.  Script  letters  denote  R-vector  spaces.  The 
range  space  of  a  map  B  is  denoted  ^(B)  «*J i;  the  null  space  is  denoted  a[B]  . 

The  dimension  of  a  vector  space  *  Is  written  as  d(£). 

Let  A:  and  take  a  subspace  V^X.  Define 

A  V  ^  Ly  I  y  •  Ax  and  x  € tri 

A  V  •  lx  €  £  I  A  x  €  V  J 

Tilt  aubspace  V  1*  said  to  be  an  Invariant  subspace  of  A  if  In  addi¬ 
tion,  let  B:  and  F:  X  “•  Now  V  Is  said  to  be  an  (A ,  B)  -  Invariant 

aubspace  If  there  exists  a  matrix  F  such  that  (A  +  BF)V<-V»  •'  Is  an 


Invariant  subspace  of  A  +  BF. 

Let  >1,  J  be  subspaces  of  Z.  Define 

ft+a'  ■  tx  I  r  +  s  •  x,  rt'i,  s 

If  »!  +  «/■£  and  if  i"V-0,  then  »i  and  ^  form  a  direct  sum  decompos ition  of  % 
and  it  is  denoted 

Let  A:  X  -Z  and  Let  y  be  a  subspace  of  X.  CaLl  the  vectors  x, 
y€X  equivalent  mod  V  if  x-y€y.  Define  the  factor  space  Zly  as  the  set 
of  all  equivalence  cLasses 

x  -  i.y  I  y  gX,  x  -  y  €VJ ,  x  £y 

To  turn  Z/y  into  a  linear  vector  space,  define 

x1  +  x2  -xL  +  x;,  C  xL 

for  x^,  € Z,  end  c  a  rea  1  scalar .  Note  that  Z/y  is  not  a  subspace  of  X, 

Now  suppose  y  is  an  invariant  subspace  of  A  such  that  d(^-)  -j  and 
a+p“d(%).  Then  chose  a  subspace  ft  of  Z  such  that  Chose  a 

bas is  t-Tj!  J  ■  l, . . .  ,pi  for  ft  and  a  bas is  1  ;  i  -  1 . jJ  for  V.  In  the 

basis  lr^,...,vcrJ  the  matrix  A  has  the  form 


The  factor  space  Z/y  is  isomorphic  to  >i  and  Apxp  is  the  map  induced  in  Z/y 
4 

by  A.  The  matrix  Aj^j  is  called  the  restriction  of  A  to  and  is  denoted 
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(A  I  V)  ■  • 

For  complete  detalLs  of  these  ideas,  see  [11].  For  a  short  re¬ 
view  of  the  geometric  concepts,  in  particular  the  (A, B>- invariant  subspace, 
see  [12] . 

2.2.  Representations  of  Linear  Systems 

Continuous  and  discrete  time  systems  are  introduced  here.  Their 
interconnection  is  discussed  to  motivate  the  geometric  presentation  below. 

It  will  be  assumed  that  the  continuous  time  systems  are  linear 
time -invariant  systems  of  the  form 


x(t)  =  Ax(t)  +  Bu (t)  (2.2,1a) 

y(t)  =  Cx(t)  +  Du(t)  (2.2,1b) 

where  x£R  ,  u  €  Rm,  y  €  Rr  and  (A,B,C,D)  are  real  constant  matrices  of  the 
appropriate  dimension.  The  frequency  response  characteristics  of  this 
system  can  be  studied  by  taking  the  Laplace  transform  of  (2.2,1).  Assuming 
x(0)s,0,  this  yields 

sx(s )'  »  Ax(s)  +  Bu (s)  (2.2,2a) 

y(s)  =  Cx(s)  +  Du(s)  (2.2,2b) 


It  is  convenient  to  represent  (2.2,2)  in  the  form 


si  -  A  -B 

x(s) 

x(s) 

0  """ 

-  P(s) 

3 

_C  D_ 

_u(sj_ 

u(s) 

y(s) 

(2.2,3) 
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The  matrix  P(s)  is  known  as  the  system  matrix  [13].  Under  the  assumption 
of  zero  initial  conditions,  it  contains  all  of  the  information  about  the 
internal  behavior  of  the  system  that  (2.2,2)  contains.  If  only  the  input- 
output  behavior  of  the  system  is  of  interest,  the  state  can  be  eliminated 
and  (2.2,2)  can  be  solved  for  y(s)  in  terms  of  u(s): 

y(s)  -  [C(sl-A)  ^B  +  D]u(s)  =G(s)u(s)  (2.2,4) 

The  matrix  G(s)  is  called  the  transfer  function  matrix.  The  transfer 
function  matrix,  the  system  matrix,  and  the  system  (2.2,1)  will  be  the 
starting  points  for  studying  zeros  of  continuous  time  systems. 

Consider  instead  of  (2.2,1),  the  discrete  system 

Xi+1  *  Axi  +  Bui  (2.2,5a) 

yL  -  cx.  +  Du.  (2.2,5b) 

where  x£Rn,  u  £  Rm,  y  g  Rr  and  (A,B,C,D)  are  real  constant  matrices  of  the 
appropriate  dimension.  By  using  z-transforms  in  an  analogous  way  to 
Laplace  transforms,  the  system  (2.2,5),  assuming  zero  initial  state,  can 
be  transformed  to 


zx(z)  =  Ax(z)  +  Bu(z) 
y(z)  *  Cx(z)  +  Du( z) 
As  above,  this  can  be  written  as 


zI-A  -B 

x(z) 

-  P(z) 

x(z) 

m 

0 

C  D 

u(z) 

u(z) 

y(z) 

—  _ 

—  — 

_  __ 

_ 

(2.2,6a) 

(2.2,6b) 


(2.2,7) 


In  a  similar  fashion,  the  transfer  function  matrix  of  (2.2,6)  is  calculated 
as 

y(z)  -  [C(zl-A)  ^8  +  D]u(z)  -  G(z)u(z)  (2.2,8) 

Comparing  (2.2,7)  and  (2.2,8)  to  (2.2,3)  and  (2.2,4),  respectively, 
it  is  seen  that  the  functional  form  is  exactly  the  same.  Hence,  any  de¬ 
finitions  of  zeros  based  on  P(s)  or  G(s)  will  apply  equally  well  to  their 
counterparts  in  discrete  systems.  Much  of  the  rest  of  the  presentation 
is  done  in  terms  of  a  transformed  variable,  usually  written  s.  However,  this 
variable  should  be  interpreted  as  either  the  Laplace  transform  variable,  or 
the  z-transform  variable,  which  ever  is  appropriate.  In  a  few  instances, 
the  topic  under  discussion  will  be  limited  to  either  a  continuous  or  discrete 
system.  This  will  be  pointed  out  in  the  text. 

It  would  be  desirable  to  extend  this  analogy  to  systems  represented 

in  the  forms  (2.2,1)  and  (2.2,5).  To  do  this  it  is  necessary,  however,  to 

explicitly  note  the  distinction  between  a  reachable  state  and  a  controllable 

state  of  a  system.  For  the  system  (2.2,5),  let  ^(i;iQ,x^)  be  the  transition 

function  from  the  stace  x  at  time  i  to  the  state  x->>(i;i  ,x  )  at  time  i. 

o  o  o  o 

A  state  x  is  controllable  if  there  exists  an  integer  i  such  that 

<?(i;0,x)  -  0  (2.2,9) 

a  state  x  is  reachable  if  there  exists  an  integer  i  such  that 

$(0;i,0)  =■  x  (2.2,10) 


Corresponding  definitions  can  be  given  for  continuous  time  systems.  The 


sec  of  controllable  (reachable)  states  form  a  subspace  in  state  space. 

For  continous  systems,  the  reachabLe  space  and  the  controllable  space  are 
the  same.  This  is  not  true  for  discrete  systems;  however,  it  is  true  that 
the  reachable  space  is  always  contained  in  the  controllable  space.  Similar 
remarks  hold  for  Che  observable  and  detectable  subspaces.  It  will  be  shown 
that  if  the  reachable  and  observable  subspaces  are  employed  in  characterizing 
zeros,  then  results  obtained  from  the  system  representations  (2.2,1)  and  (2.2,5) 
are  the  same.  Hence,  in  this  work  only  the  reachable  and  observable  sub¬ 
spaces  will  be  used. 

2.3.  Smith  Form  and  Smith-McMillan  Form 

Two  of  the  three  representations  from  which  zeros  will  be  de¬ 
fined  are  matrices  whose  elements  are  either  polynomials,  as  in  the  system 
matrix  P(s)  (2.2,3),  or  rational  functions,  as  in  the  transfer  function 
matrix  G(s)  (2.2,4).  The  analysis  of  the  properties  of  P(s)  or  G(s)  is  done 
here  via  two  canonic  matrix  forms  called  the  Smith  form  (for  polynomial 
matrices)  and  the  Smith-McMillan  form  (for  matrices  of  rational  functions). 
This  section  will  present  a  brief  development  of  these  two  forms.  The 
interested  reader  is  referred  to  [11]  and  [13]  for  more  details  and  other 
properties . 

Consider  a  nxm  polynomial  matrix  A(s)  of  rank  r.  To  obtain  the 
Smith  form  of  A(s) ,  define  the  following  row  operations  on  A(s): 

1)  Multiplication  of  any  row  by  a  non- zero  constant 

2)  Addition  to  any  row  of  any  other  row  multiplied  by  an 
arbitrary  polynomial. 

3)  Interchange  of  any  two  rows 
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In  a  slmiliar  fashion,  define  elementay  column  operations. 

Now  select  the  element  of  A(s)  which  has  least  degree.  Suppose 
that  it  is  aij(s).  Then  for  any  other  element  a^(s) 

ahk(s)  "  +  rhk(3)  (2.3,1) 


where  q(s)  is  a  polynomial  such  that  the  degree  of  r^k(s)  Is  zero  or  the 
degree  of  rhk(s)  is  less  than  the  degree  of  a^(s).  If  a^s)  doesn't 
divide  every  element  of  A(s),  then,  by  using  elementary  row  and  column 
operations,  A(s)  can  be  replaced  by  a  matrix  A^(s)  which  contains  an 
element  of  degree  lower  than  a^Cs),  say  r^^(s).  Now  continue  this  pro¬ 
cess  using  rhk(s)*  After  a  finite  number  of  steps  X,  the  matrix  A.(,s) 
will  contain  an  element  which  divides  every  other  element  of  A.(s).  Again 

M 

using  elementary  row  and  column  operations,  this  lowest  degree  element  can 
be  placed  in  the  upper  left  hand  corner  and  the  other  elements  in  the 
first  row  and  column  reduced  to  zero.  Further,  this  element  can  be  chosen 
to  be  monic.  Denote  this  matrix  by  B(s).  Then  B<,s)  has  this  form: 


B(s) 


bu(s) 


(2.3,2) 


where  b^(s)  is  a  monic  polynomial  which  divides  every  element  of  B^s). 

Now  repeat  the  process  with  B(s).  Eventually  this  procedure 
must  terminate  as  A(s)  has  finite  rank.  The  matrix  so  obtained,  call  it 
S(s),  is  said  to  be  the  Smith  form  of  A(s)  and  its  structure  is  given  by 


S(s)  -  diagf^Cs),  €2(*>  ,  •  -  -  ,€r(s>  ,  0 . 0) 


(2.3,3) 


where  by  construction 


€i-l(s)  ^  6i(s>  i *  2 , . . . ,r 


(2.3,4) 


The  polynomials  €^(s) ,  i *  1, . . . ,r  are  called  the  invariant  factors  or 
invariant  polynomials  of  A(s). 

Each  of  the  elementary  row  and  column  operations  given  above 
can  be  represented  by  pre-  and  post-  multiplication  of  A(s)  by  the 
appropriate  matrices,  respectively.  For  instance,  the  interchange  of  two 
rows,  the  ith  and  jth  say,  is  accomplished  by  pre-multiplication  of  A(s) 
by  a  permutation  matrix  obtained  from  an  identity  matrix  which  has  the 
ith  and  jth  rows  interchanged.  It  should  be  noted  that  each  matrix  re¬ 
presenting  an  elementary  operation  is  a  non-singular  unimodular  matrix. 

It  therefore  follows  that  the  matrices  A(s)  and  S(s)  are  related 


by 


S(s)  =  L(s)A(s)R(s) 


(2.3,5) 


The  nxn  matrix  L(s)  is  a  product  of  matrices  which  represent  elementary 
row  operations.  Since  each  matrix  is  unimodular,  L(s)  is  also  unimodular. 
Similar  remarks  apply  to  R(s).  Two  matrices,  such  as  A(s)  and  S(s),  re¬ 
lated  as  in  (2.3,5)  are  called  equivalent  matrices. 

Since  L(s)  is  a  unimodular  matrix,  its  local  rank  for  any  value 
of  s  is  always  equal  to  its  normal  rank.  The  same  is  true  of  R(s).  It 
follows  that 


P  C  S (s ) ]  =  P[A(s)] 


(2.3,6) 


for  all  s. 


The  special  form  of  S(s)  makes  this  a  very  useful  property. 


1  5s  +2 

_s  +  l  (s  -  2)^_ 


In  terms  of  equation  (2.3,5),  the  relationship  is  given  as: 


L6 


The  Smith  conn  has  many  important  and  interesting  properties. 
However,  only  the  relationship  between  the  invariant  polynomials  and  the 

’  *  *  ’  Sc. 

deteraencal  divisors  will  be  of  interest  here.  Lee  A  ' 


j  1»  *  *  *  *  *  ■} 


;s 

k  I 


denote  the  k-ch  order  minor  of  A(s)  wnich  is  formed  by  deleting  all  rows 
but  i^.i^,..-,!^  and  all  columns  but  j^, . . . ,  The  minors  of  S(s) 

are  related  to  the  minors  of  A(s)  by  the  Binet-Cauchy  formula  [11,  d.  9]  : 


^1 ‘ i2 ’ *  *  * * 

SI  1  *  ?:  s 


j  1>  j?  •  •  *  ‘ 


P  ' 


a  Sn 
P 

3  sm 
P 


Vi¬ 


la  at  ^ 

-.51  Aj  P;S 

®l>“  •  ,Qp 


31’ 


1’ 


•’pp 
•  >  3, 


;s  (2.3,7) 


?  »  1,2, . . . ,min(n,m) 

Here  the  sum  is  over  all  possible  permucations .  Let  D^(s),  j  »  1,2 . r 

denote  the  greatest  common  divisor  of  all  j-th  order  minors  of  A(s).  Define 

D0(s)  a  1.  For  Che  moment,  denote  the  greatest  common  divisors  of  the  jth 

* 

order  minors  of  S(s)  by  D^Cs).  It  follows  from  (2.3,7)  that 


D,(*)ID  (•)  j  *  1,2, . . . ,r 


(2 .3,8) 


(The  minors  of  orders  higher  than  r  are  rero  for  3(s).).  However,  L(s) 
and  R(s)  are  invertible.  Therefore,  A(s)  and  S(s)  can  change  places  in 
(2.3,7)  when  L(s)  and  R(s)  are  replaced  by  the  appropriate  matrices.  It 
now  follows  chac 


Dj(»)iDJ(s)  J  *  1,2 , . . . , r 


(2.3,9) 
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When  Che  decermencal  divisors  are  required  to  be  monic,  ic  is  seen  from 
(2.3,8)  and  (2.3,9)  that 

Dj(s)  -  D*(s)  j  ■  1,2 , r  (2.3,10) 

Ic  follows  from  Che  special  form  of  S(s)  chac 

v*>  "61(S)*  D2(3)  ■€1(s)€2(s),...,Dr(s)  -€1(s)...gr(s)(2.3,ll) 
Noce  chac  if  S(s)  is  square  and  has  full  rank  Chen 

Dr(s)  -  dec [S (s ) ]  (2.3,12) 

This  implies  chac  for  a  square  matrix  of  full  rank,  say  A(s),  chac 

dec[A(s)]  -  ®€l(s)€2C«)...  €r(«)  (2.3,13) 

where  ®  is  a  non-zero  conscanc. 

Example  2 

Lee  A (s )  be  defined  as  in  Example  1.  By  inspection,  the  greatest 
commom  divisor  of  all  single  elements  is  1.  Hence, 

<=L(s)  -  1 

The  determinant  of  A(s)  is  easily  calculated  as 

det[A(s)]  *  Ss“  +  22s  -4 
It's  monic  greatest  common  divisor  is 

s"  +  ll/4s  -  ~  (s) 

4»  4 


IS 


This  agrees  with  the  Smith  form  calculated  in  Example  1.  Q 

The  analogue  of  the  Smith  form  for  matrices  whose  elements  are 
rational  functions  is  called  the  Smith-McMillan  form  (or  McMillan  form). 

Consider  now  an  rxm  matrix  G(s)  whose  elements  are  rational 
functions.  Let  the  normal  rank  be  k.  The  matrix  G(s)  can  be  written  as 

G(S)  "  d(7T  N(S)  (2.3,14) 

where  d(s)  is  the  monic  least  common  demoninator  of  all  elements  of  G(s) 
and  N(s)  is  a  rxm  polynomial  matrix.  Now  compute  the  Smith  form  of  N(s), 
call  it  Ns(s).  Since  N(s)  and  Ns(s)  are  equivalent  matrices,  they  are 
related  as 


N(s)  -  L(s)N#(s)R(s)  (2.3,15) 

where  L(s)  and  R(s)  are  unimodular  matrices.  The  matrix  G(s)  can  now  be 
expressed  as 

G(S)  “d^T  L(s)Ns(s)R(s)  "  L(s)M(s)R(s)  (2.3,16) 

where 


M(s) 


_1 _ 

d(s) 


12.3,1 7) 


and  all  possible  cancellations  have  been  carried  out  in  M(s).  The  matrix 
M(s)  is  then  the  Smith-McMillan  form  of  G(s).  It  has  the  form. 


M(s) 


diag 


€x(s)  c2(s) 

’^(S)  ’  »2(s) 


(2.3,18) 
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By  construction,  S^(s)  and  1i(s),  are  coprime.  Recall  that 

€i_i<s)  l6i(s>,  1-2,. ...k.  It  follows  that  l(s)  I  '^^(s) ,  i»2,...,k. 

It  should  be  noted  here  that  although  €^(s)  and  1^(s)  are  coprime,  £^(s) 


and  t.(3)  may  not  be  coprime  for  i t j . 
Example  3 


Let  G(s)  be  given  by 


(s+1) (s+2) 
2 

s  +s-4 


(s+1)  (s+2) 


(s+1) (s+2) 
2 

2s  -s-8 
(s+1) (s+2) 

s  +  1 


(s+1) (s+2)  s2+s-4  2s2-s-8 


2  2 
s  -4  2s  -8 


The  Smith  form  of  N(s)  is  easily  calculated  as 


N  (s)  *  0  (s+2) (s-2) 

S 


Now  the  Smith -McMillan  form  of  G(s),  M(s),is  found  to  be 


M(s)  => 


(s+1)  (s+2) 


0 


0 


CHAPTER  3 


Zeros  Defined  from  the  Transfer  Function  Matrix 

3.1,  Introduction 

Several  authors  have  defined  zeros  from  the  transfer  function 
matrix  [ L3] [ 14] [ 15] .  The  definition  of  zeros  for  a  transfer  function 
matrix  is  a  generalization,  in  some  sense,  of  the  zeros  defined  for  a 
scalar  transfer  function.  Thus,  the  properties  of  zeros  defined  for  the 
transfer  function  matrix  will  be  analogous  to  the  properties  of  zeros  of 
a  scalar  transfer  function.  In  particular,  the  zeros  are  certain  fre¬ 
quencies  at  which  a  non-zero  input  produces  an  identically  zero  output, 
an  analogy  to  the  definition  of  zeros  of  a  scalar  transfer  function. 

In  general,  the  transfer  function  matrix  G(s)  can  be  factored 
into  polynomial  matrices  such  chat  G(s)»N(s)D  *(s).  Then  the  zeros  are 
shown  to  cause  a  loss  of  rank  in  the  "numerator"  of  G(s),  N(s)  [14].  If 
G(s)  is  used  to  define  an  "inverse  system"  the  zeros  of  G(s)  are  found  to 
be  the  poles  of  the  inverse  system  [14].  In  a  recent  development  [15], 
of  generalized  root-locus  and  Nyquist-Bode  techniques  for  multivariable 
systems,  the  zeros  of  G(s)  are  shown  to  play  an  analogous  role  to  their 
part  in  the  classical  theory.  Clearly,  zeros  for  multivariable  systems 
have  properties  which  are  analogous  to  properties  of  zeros  of  scalar 
transfer  functions. 

It  turns  out  that  none  of  their  properties  provide  a  very 
attractive  computation  method  for  obtaining  these  zeros.  However,  the 
Smith-McMillan  form,  which  provides  the  theoretical  basis  for  this  investi¬ 
gation,  does  provide  a  discription  of  zeros  which  is  somewhat  more  managable 
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computationally  than  the  other  properties. 

The  organization  of  this  chapter  is  as  follows;  Section  2 
provides  a  motivation  and  definition  of  zeros.  This  is  accomplished  by 
the  Smith -McMillan  form  of  G(s).  This  same  property  is  again  described 
in  Section  3  from  a  slightly  different  mathematical  approach.  It  is  then 
shown  that  the  same  set  of  zeros  is  being  described.  In  Section  4  the 
zeros  of  G(s)  are  described  as  the  poles  of  an  inverse  system.  Section  5 
provides  a  characterization  of  zeros  in  terms  of  the  minors  of  G(s). 
Sections  2-5  depend  fundamentally  on  the  Smith -McMillan  form  of  G(s). 
Section  6,  then,  briefly  describes  the  zeros  of  G(s)  in  terms  of  complex 
variable  theory.  (This  is  related  to  generalized  Nyquist-Bode  and  root- 
low  criteria  for  multivariable  systems).  Section  7  concludes  with  a 
summary  of  the  material  presented.  Examples  to  illustrate  the  theory  arc 
provided  throughout  the  chapter. 

3.2.  Definition  of  Transmission  Zeros 

Consider  a  r xm  matrix  of  rational  functions,  G(s).  Without 
loss  of  generality,  let  the  normal  rank  of  G(s)  be  min(m,r).  The  transfer 
function  matrix  G(s)  has  a  Smith -McMillan  form 

G(s)  -  L(s)M(s)R(s)  (3.2,1) 

where 

^i*'3  V 

M(s)-diag[  Vt(3)]  i  -  l, . . .  ,min(m,r)  (3.2,2) 

and  L(s)  and  R(s)  are  unimodular  matrices.  Since  L(,s)(Rfs))  has  a  non¬ 
zero  constant  determinant,  its  columns  (rows)  are  linearly  independent 
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for  all  s.  Let 


i)  ^(s)  i  p  be  the  columns  of  L(s) 


and 


11)  r^(s)  1  p  be  the  rows  of  R(s) 


where  p-min(m,r).  Using  this  notation,  G(s)  can  be  expressed  as 


P  <=,(3) 

G(s)  l.( s)  ^  rt(s) 


(3.2,3) 


Since  G(s)  relates  the  transformed  input  and  output  vectors  by 


y(s)  =  G(s)u(s) 


(3.2,4) 


it  follows  that 


€,(s) 


y(s)  -  Z^  it(s)  ^  rt(S>u(.) 


(3.2,5) 


The  relation  (3.2,5)  shows  explicitly  the  relationship  between  the  input 
and  output  for  any  specific  frequency.  Suppose  that  u(s)  is  chosen  such 
that 


Tj(s)u(s)  -  5ij^(s) 


(3.2,6) 


where  a(s)  is  a  non-zero  polynomial  and  6^  is  the  Kronecker  delta.  Then 
relation  (3.2,5)  shows  that 


i)  y(s)  is  zero  if  and  only  if  ^(s)  is  zero 


y(s)  is  infinite  if  and  only  if  ^(s)  is  zero 


ti) 


This  argument  qualitatively  describes  the  main  idea  that  provides  the 

★ 

physical  motivation  to  the  following  definitions  of  zeros  and  poles.  That 
is  to  say,  if  a  complex  frequency  is  a  root  of  <^(s)  for  some  i,  then 
transmission  of  that  frequency  through  the  system  is  blocked.  If  a  complex 
frequency  S2  is  a  root  of  t^s)  for  some  i,  then  the  output  becomes  in¬ 
finite  at  that  frequency.  These  properties  are  similar  to  properties  of 
zeros  and  poles  for  a  scalar  transfer  function. 

Therefore,  let 

P 

z(s)  -  tt  €  (s)  (3.2,7) 

i-l 


and 


P 

p(s)  -  rr  *  (s)  (3.2,8) 

i-l 


The  polynomial  z(s)  will  be  refered  to  as  the  zero  polynomial  and  p(s) 
will  be  refered  to  as  the  pole  polynomial  [17].  In  general. 


z(s)  -  (s-Zl)(s-z2)...(s-zk)  (3.2,9) 

where  among  the  complex  numbers  z^,  i»l,...,k  any  of  them  may  be  repeated. 
Definition  1  [ 13] 

The  transmission  zeros  of  G(s)  are  the  complex  roots  of 

z(s) . 

□ 


This  development  appears  in  [14], 
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That  is  to  say,  the  transmission  zeros  of  G(s)  are  the  complex  numbers 


z^,  i-l,...,k.  The  zero  polynomial  can  be  written  as 


2(3) 


t  k 

n  (s  -  zt)  ,  kt  >  0,  k 


i-1 


t 

E  k 

i-1 


(3.2,10) 


where  the  complex  numbers  z^,  i-l,...,t  are  required  to  be  distinct.  In 
this  case,  G(s)  is  said  to  have  a  transmission  zero  of  order  at  z^. 

In  the  sequel,  it  shall  be  useful  to  talk  about  the  poles  of 
G(s).  Although  the  poles  of  a  system  are  not  the  primary  topic  of  this 
paper,  a  definition  appears  here  for  future  reference.  The  comments  above 
about  the  multiplicities  of  zeros  apply  also  to  poles. 

Definition  2  [17] 

The  poles  of  G(s)  are  the  complex  roots  of  p(s).  a 

Example  4 

Consider  the  matrix  G(s)  given  in  Example  3.  It  has  the  Smith- 
McMillan  form 


u 


D 


M(s)  - 


1 


(s+1)  (s+2 ) 
0 

0 


s-2 

s+1 


Hence , 


z(s)  =  s-2 


P ( s )  -  (s+1)  (s+2) 


0 


This  system  has  a  transmission  zero  at  s  *  2  and  poles  at  s  ■  -1  and  s  ■  -2.  The 
pole  at  s  • -1  has  a  multiplicity  of  2.  Q 

3.3.  Another  Characterization  of  Transmission  Zeros 

In  the  heuristic  argument  in  the  last  section,  a  transmission 
zero  of  a  transfer  function  matrix  G(s)  was  shown  to  be  a  complex  frequency 
whose  transmission  through  the  system  was  blocked.  This  was  the  physical 
motivation  for  defining  the  transmission  zeros  of  G(s)  in  [14]  and  [15]  as 
well.  However,  the  mathematical  approach  was  slightly  different  in  chose 
references . 

Again  consider  a  rxm  matrix  of  rational  functions  G(s)  which  has 
full  rank.  It  is  well  known  [14]  that  G(s)  can  be  factored  (non-unique ly) 
as 

G(s>  -  D]1(s)Ni(s)  -  Nr(s)D‘L(s)  (3.3,1) 

where  N.(s)  and  N  (s)  are  rxm  polynomial  matrices,  D  (s)  is  a  mxm 

"  4  IT 

polynomial  matrix,  and  D^(s)  is  a  rxr  polynomial  matrix.  Furthermore, 

Na(s)  and  D^(s)  are  left  coprime  and  Nr(s)  and  Dr(s)  are  right  coprime. 

The  matrices  D_(s)  and  D  (s)  are  required  to  be  non-singular.  The  assump- 

L  At 

tion  that  G(s)  has  full  normal  rank  implies  that  both  N  (s)  and  N  (s)  have 

XT  JL 

full  normal  rank. 

Consider  the  case  where  r^m.  Suppose  there  exists  a  complex 


number  z  such  that 


and  that  z  is  not  a  pole  of  G(s).  Then  in  [14],  it  is  shown  that  there 
exists  a  nonzero  input  vector  g  and  an  appropriate  set  of  initial  conditions 
such  that  for  the  input 

u(t)  -  ge2tU(t)  (3.3,3)* 

the  response  of  the  system  is  identically  zero  for  all  t>  0.  A  similar 
result  holds  when  m  >  r  and  there  exists  a  complex  number  z  (not  a  pole) 
which  satisfies  (3.3,2). 

These  observations  suggest  that  the  complex  number  z  described 
above  is  related  to  the  transmission  zeros  of  G(s).  To  see  this  connection, 
consider  the  Smith -McMillan  form  of  G(s) 
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rx(s)  *  diagf  (s) , . . . Cs) , 1 , -  -  - , L]  (3.3,8) 

If  p<  m,  the  I's  are  placed  on  the  diagonal  to  make  T^(s)  invertible.  The 
mxm  polynomial  matrix  rr(s)  has  a  similar  form.  Using  this  factorization, 
G(s)  can  be  expressed  as 

G(s)  * L(s)r^1(s)E^(s)R(s)  *  L(s)Er(s)f"1(s)R(s)  (3.3,9) 

Comparision  of  (3.3,9)  to  (3.3,1)  shows  that  (3.3,9)  is  just  a  special  case 
of  (3.3,1)  where 

VS)  -  ^(s)L'L(s)  ;  N*(s)  A  Ei(s)R(s)  (3.3,10) 

D*(s)  A  R'1(s)rr(s)  ;  N*(s)  A  L(s)Er(s)  (3.3,11) 

Since  R(s)  and  L(s)  are  unimodular  matrices,  if  for  some  complex  number  z 

P[N*(z)]  <  Pn[N*(-)]  (3.3,12) 

this  implies 

P[E/z)]  <  PntE<fc(-)]  (3.3,13) 

Now,  from  the  special  structure  of  E^(s)  (see  equation  (3.3,7)),  it  is 
seen  that  (3.3,13)  is  true  only  if  z  is  a  root  of 

P 

z(s)  =■  ^  6i(s)  (3.3,14) 

The  same  argument  holds  for  N*(s).  Finally,  in  [15]  it  is  shown  that  if 
N^(s)  and  N^(s)  are  any  two  matrices  which  satisfy  the  conditions  for 
factorization  in  (3.3,1),  then  N(s)  and  N  (s)  are  equivalent.  In  particular, 
all  matrices  N^(s)  which  satisfy  the  factorization  conditions  in  (3.3,1)  are 
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equivalent  to  N^(s)  defined  in  (3.3,10).  By  an  argument  similar  to  the 
one  above,  it  is  seen  that 

P£NX(*)J  <  P[NX(-)]  (3.3,15) 

only  if  z  is  a  root  of  (3.3,14).  Similar  remarks  hold  for  matrices  Nr(s) 
which  satisfy  the  factorization  conditions  of  (3.3,1).  From  (3.3,10)  and 
(3.3,11)  it  is  clear  that  those  numbers  z  which  satisfy  (3.3,15)  will 
also  satisfy 


P[Nr(z)]  <  P[Nr(- ) ]  (3.3,16) 

and  conversely.  Thus,  the  following  theorem  has  been  proved; 

Theorem  1  [ 14] 

A  complex  number  z  is  a  transmission  zero  of  G(s)  if  and  only 
if 

P[Nx(z)]  <  Pn[NiC)], 

or,  equivalently, 

P[Nr(z)]  <  Pn[NrC)l, 

where  N  (s)  (N  (s))  is  obtained  from  any  factorization  of  G(s)  given  by 
(3.3,1).  3 

Two  remarks  are  in  order.  First,  note  that  Theorem  1  does  not 
provide  a  way  of  calculating  the  order  of  a  particular  zero;  i.e.  it  is 
not  possible  to  calculate  the  numbers  k^  in  equation  (3.2,10).  Hence,  it 
is  said  that  Theorem  1  does  not  allow  for  multiplicities. 


Secondly,  suppose  that  none  of  the  transmission  zeros  is  also  a 


pole  of  G(s).  Then  it  is  seen  from  equation  (3.3,9)  that 


P[G(z)]  <  Pn[G(-)] 


(3.3,17) 


when  z  is  a  transmission  zero.  However,  if  a  pole  and  a  transmission  zero 

ie 

coincide  ,  then  it  is  not  possible  to  assert  (3.3,17).  Indeed,  it  is  not 
even  clear  how  to  define  the  rank  of  G(s)  as  some  terms  of  G(s)  are  going 
to  infinity,  as  s  approaches  the  transmission  zero  [18]. 


Example  5 


Consider  the  matrix  G(s)  given  in  Example  4.  There  it  was  shown 
that  z  =  2  is  a  transmission  zero  of  G(s)  by  obtaining  the  zero  polynomial. 
This  may  be  checked  by  Theorem  1.  It  is  easily  seen  that 


and  clearly 

1»P[G(2)]  <  Pn[G(s)]  =  2 

Hence,  z=2  is  a  transmission  zero  of  G(s)  by  Theorem  1.  S 

★ 

That  this  is  possible  follows  from  the  remarks  on  the  Smith -McMillan 
form  in  Section  2.3  and  the  definition  of  the  pole  and  zero  polynomial. 
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3.4.  Transmission  Zeros  and  Inverse  Systems 

The  proof  of  Theorem  l  in  the  last  section  provides  another 
characterization  of  transmission  zeros.  Consider  again  (3 .3 , 1) , 

G(s)  -D'l(s)N  ,(s)  -H  (aJD"1^)  (3.4,1) 

x  x  r  r 

where  G(s)  is  a  transfer  function  matrix  and  the  matrices  N  (s)  and 

JL 

D  ,i>)  are  a  factorization  of  G(s)  described  in  Section  3.3.  In  the  same 
x 

way  it  is  snown  that  any  two  matrices  N  .(s)  and  N  (s)  are  equivalent,  it 

JL  a. 

can  be  shown  [15]  that  all  matrices  D^(s)  are  also  equivalent.  Recall 
the  special  factorization  of  G(s)  obtained  from  (3.3,9): 

G(s)  -  L(s)r"l(s)E,(s)R(s).  (3.4,2) 

X  * 

That  is. 


D*(s)  -  ' 1  ,  N*(s)  «E  (s)R(s), 


(3.4,3) 


Hence,  all  matrices  D  is)  are  equivalent  to  D.(.s),  i.e. 

A  A 


D  ( s )  -  D(s)D  (s)V(s) 

X 


(3.4,4) 


wnere  U(s)  and  V(s)  are  unimodular  matrices.  It  follows  that 

-1  p 

det[D,(s)]  -det[U(s)r(s)L  (s)V(s)]«a  n  *  (s)  (3.4,5) 

1  i-1 

wnere  3  iS  a  non-zero  constant.  This  shows  chat  the  roots  of  Che  deter¬ 
minant  of  D,(s)  are  the  poles  of  G(s)  (see  Definition  2).  This  leads  to 
another  parallel  between  transmission  zeros  defined  for  G(s)  and  zeros 
defined  for  scalar  transfer  functions;  the  zeros  of  a  system  are  the  poles 


of  Che  inverse  svscem 


To  make  this  precise  define  a  reflexive  generalized  inverse  [19] 


of  a  matrix  X  as  Che  matrix  X  wnich  satisfies 


The  reflexive  generalized  inverse  G(s)  of  G(s)  can  be  obcained  from  (3.4,2) 


by  inspection  as 


where  E,(s)  is  obtained  from  E  (s)  by  transposing  E.(s)  and  inverting 


Che  non-zero  elemencs.  Using  equations  (3.4,2)  and  (3 


This  gives  the  following  theorem 


Theorem  2  [ 14] 


The  transmission  zeros  of  G(s)  are  the  poles  of  G  (s)  where 


G+(s)  is  the  reflexive  generalized  inverse  of  G(s) 


Example  6 


Consider  a  rxm  transfer  function  matrix  of  full  rank  where 


rSm.  Then  there  exists  a  right  pseudo-inverse  of  G(s),  G  (s)  such  that 


Clearly,  G  (s)  satisfies  (3 


6).  Therefore,  the  zeros  of  G(s)  are  the 


A  similar  result  holds  for  m£r 


where  d(s)  is  the  monic  least  common  denominator  of  all  elements  of  G(s), 


and  N(s)  is  a  polynomial  matrix.  Let  the  Smith  form  of  N(s)  be  N  (s) 


Let  the  Smith -McMillan  form  of  G(s),  M(s),  be  given  by 


of  all  j-square  minors  of  N(s)  where  D  is  defined  to  be  1.  As  was  shown 


It  follows  that 


That  is  to  say,  d(s)  is  the  monic  least  common  denominator  of  all  1  by  1 
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D2(s)  sl(s)s2(s)  6l(s)62(s) 
d2(S)  d2(S)  ^(3)^(8) 


(3.5,7) 


In  the  right  hand  term,  cancellations,  if  they  occur,  are  between  '^(s) 
and  £>(s).  This  is  because  of  the  dividing  properties  of  the  invariant 
factors  and  the  fact  chat  £^(s)  and  ^(s)  are,  by  definition,  coprime. 

Let  ;t,(s)  be  the  polynomial  obtained  from  ’^(s)'^(s)  after  all  possible 
cancellations  have  been  made.  Since  D0  (s)  is  the  greatest  common  divisor 
of  all  2x2  minors  of  N(s),  it  follows  that  "^(s)  is  the  least  common 
denominator  of  all  2  x2  minors  of  G(s).  However,  all  of  the  factors  in 
(s)  appear  in  7(s).  It  follows  that  ^(s)  '^(s)  is  the  least  common  de¬ 
nominator  of  all  lx  land  2x2  minors. 

By  continuing  this  argument  in  the  obvious  way,  it  is  found  that 
p(s),  the  pole  polynomial,  is  the  least  common  denominator  of  all  minors  of 
all  orders  of  G(s).  It  follows  similarly  that  the  zero  polynomial  is  the 
greatest  common  divisor  of  the  numerators  of  all  minors  of  maximum  order  of 
G(s)  provided  these  minors  have  been  adjusted  to  have  p(s)  as  their  denomina¬ 
tor.  These  results  are  summarized  in  the  following  theorem: 

Theorem  3  [17] 

The  poles  of  G(s)  are  the  complex  roots  of  the  least  common 
denominator  of  all  non-zero  minors  of  all  orders  of  G(s). 

The  transmission  zeros  of  G(s)  are  the  roots  of  the  greatest 
common  divisor  of  the  numerators  of  all  non-zero  minors  of  maximum  order  of 
G(s)  provided  they  have  been  adjusted  to  have  p(s)  as  their  denominator.  D 

From  the  above  development  it  is  also  clear  that  the  order  of  a 
transmission  zero  can  be  identified  by  this  characterization. 


Example  8 


Let  G(s)  be  the  transfer  function  given  in  Example  7.  Using  the 
notation  for  minors  introduced  in  Section  2.3,  Che  transmission  zeros  of 
G(s)  can  be  calculated  by  Theorem  3  as  follows; 


G 


-  3j3:?) 

(s+1)2 (s+2) 


/2,3  \ 


„  liiiD 

(S+1)2 (3+2) 

3s (s-2 ) 
(s+1)2 (s+2) 


Now  it  is  easily  seen  that 

P (s )  -  (s+1)2 (s+2) 

and 


z(s)  -  (s-2) 


which  agrees  with  the  previous  results . 


a 


3.6.  Zeros  obtained  From  Complex  Variable  Theory 

Macfarlane  and  Postlethwalte  [16],  in  the  course  of  deriving  a 
generalized  Nyquist  stability  criterion  for  multivariable  systems,  develop 
another  characterization  of  zeros  of  a  transfer  function  G(s).  This 
characterization  is  valid  only  for  continuous  time  systems.  The  development 
will  only  be  outlined  here.  The  interested  reader  is  referred  to  [16], 
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Consider  a  mxm  matrix  of  rational  functions  G(s)  which  has 
full  rank.  For  a  specific  value  of  complex  frequency,  say  s,  the  matrix 
G(s)  has  m  eigenvalues,  q^Os)  i*l,...,m.  It  is  desired  to  express  these 
eigenvalues  as  functions  of  frequency.  Proceeding  in  an  obvious  way,  let 

det[q(s)Im-Q(s))  ^A(q,s)-0  (3.6,  i) 


Let 


d(q,s)  -qm(s)  +a1(s)qm-I(s)  +  ...+  am(s)  -0  (3.6,2) 

where  the  a^(s)  are  rational  functions.  Let  b^(s)  be  the  least  common 
denominator  of  a^(s),  i*l,...,m.  Then 

bo(s)^(q,s)  -bo(s)qm  +  bL(s)qm"l+...+  bmfs)  -0  (3.6,3) 

where  now  the  b^,  i*0,l,...,m  are  polynomials  in  s.  In  complex  analysis 
these  functions  are  called  algebraic  functions.  These  functions  have  many 
interesting  and  important  properties,  and  they  are  central  to  the  Nvquist 
criterion.  The  properties  relevant  to  this  discussion  shall  simply  be 
stated.  The  interested  reader  is  refered  to  [16]  and  [20]. 

It  can  be  shown  that  q(s)  *0  when  )  -0  and  q(s)  -*  when 
b^(s)-0.  This  suggests  these  definitions  of  poles  and  zeros  of  G(s): 
Definition  J  [  16] 

The  zeros  of  G(s)  are  the  complex  roots  of  bm(s). 

Definition  4  [ 16] 

The  poles  of  G(s)  are  the  complex  roots  of  b(s). 


To  delineate  the  relationship  of  these  zeros  with  the  transmission 


zeros  defined  by  Definition  1,  consider  this  expansion  of  the  determinant 
A(q,s): 

det[qlm  -  G(s)]  -  qm  -  [trace  of  G(s)]qm_1 +  (£  principal 

minors  of  G(s)  of  order  2]  qm  2  - • . .+  (-l)mdet  Q(s)»0 

(3.6,4) 

Hence,  bQ(s)  is  the  least  common  denominator  of  all  non-zero  principal 
minors  of  all  orders.  In  general,  bQ(s)  will  differ  from  p(s)  by  a  factor 
e(s),  i.e. 


p(s)  -  e(s)bQ(s) 


(3.6,5) 


Since  G(s)  is  a  square  matrix,  the  zeros  defined  by  Definition  1  will  be 
the  roots  of  the  numerator  of  det[G(s)]  provided  the  denominator  is 
adjusted  to  be  p(s).  Therefore, 


b  (s) 

det[G(s)]-am(s)=ir_ 


e(s)bm(s) 

P(s) 


(3.6,6) 


These  remarks  show  that  the  poles  and  zeros  defined  from  the  algebraic 
function  (3.6,3)  are  a  subset  of  the  poles  and  transmission  zeros  defined 
from  the  Smith-McMillan  form  of  G(s).  These  two  sets  of  poles  and  zeros 
differ  by  a  set  of  poles,  which  are  also  zeros,  that  are  introduced  by  the 
non-principal  minors  of  G(s), 


3.7.  Summary 

This  chapter  has  defined  the  transmission  zeros  of  a  transfer 
function  matrix  G(s).  The  motivation  for  this  definition  is  based  on  the 
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property  that  for  certain  inputs,  transmission  of  a  certain  frequency 
though  the  system  is  blocked.  Transmission  zeros  were  then  characterized 
by  a  rank  test  of  G(s)  (Theorem  1),  the  poles  of  an  inverse  matrix  for 
G(s)  (Theorem  2),  the  minors  of  G(s)  (Theorem  3),  and,  for  square  matrices 
G(s),  the  roots  of  complex  functions  obtained  from  G(s)  (Definition  4). 

If  Theorems  1,2,  or  3  are  used  to  compute  the  transmission  zeros  of  G(s), 
the  zeros  so  obtained  will  be  those  defined  by  Definition  1  up  to  multipli¬ 
cities.  Definition  4,  in  general,  defines  a  smaller  set  of  zeros. 

In  the  literature,  all  of  the  characterizations  of  transmission 

. 

zeros  contained  in  Theorems  1,2,  and  3  have  been  used  to  define  transmission 
zeros.  It  is  easy  to  see  that  if  any  one  of  these  three  theorems  is  used 
to  define  transmission  zeros,  the  characterization  given  by  Definition  1 
can  be  recovered  as  a  theorem.  It  should  also  be  noted  that  the  term 
"transmission  zero"  is  not  standard  in  the  literature.  Some  authors  use 
"transmission  zeros"  to  refer  to  a  larger  class  of  zeros  defined  from  state 
space  representations  of  system.  This  will  be  discussed  in  the  next  chapter. 


CHAPTER  4 


DEFINITIONS  OF  ZEROS  BASED  ON  STATE-SPACE  REPRESENTATIONS 

4.1.  Introduction 

The  investigation  of  zeros  of  a  state  space  system  is  carried 
out  in  a  manner  similar  to  the  analysis  of  transmission  zeros  when  the 
system  matrix  replaces  the  transfer  function  matrix  as  the  system  repre¬ 
sentation,  and  the  Smith  form  replaces  the  Smith-McMillan  form  as  the 
mathematical  tool  of  analysis.  The  motivation  for  defining  zeros  for 
state  space  systems  is  analogous  to  the  motivation  for  transfer  function 
matrices.  It  is  desired  that  the  zeros  of  a  state  space  system  be  a 
set  of  frequencies  whose  transmission  through  the  system  is  blocked. 

Several  analogous  results  are  obtained,  i.e.  the  zeros  of  the  system  matrix 
are  characterized  in  terms  of  the  invariant  polynomials,  rank  tests,  and 
the  minors  of  the  system  matrix. 

The  connection  between  transfer  function  matrices  and  state  space 
systems  is  well  known.  So  it  is  natural  to  ask  how  the  transmission  zeros 
defined  for  transfer  function  matrices  are  related  to  zeros  defined  for 
state  space  systems.  In  general,  state  space  representations  contain  more 
information  than  transfer  function  matrices.  It  is  not  surprising  then 
that  the  zeros  defined  for  state  space  systems  include  transmission  zeros 
as  a  subset,  but  they  are,  in  general,  a  larger  set.  This  requires  that  a 
new  type  of  zero  be  introduced.  It  turns  out  that  these  new  zeros  are  the 
eigenvalues  of  the  unreachable  and/or  unobservable  modes.  Thus  they  reflect 
exactly  the  information  that  does  not  appear  in  the  transfer  function  matrix. 
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Several  definitions  of  zeros  based  on  the  system  matrix  are 
introduced.  Each  of  these  definitions  describe  a  different  subset  of  the 
set  of  transmission  zeros  and  the  eigenvalues  of  the  unreachable  and/or 
unobservable  modes.  The  relationship  between  these  sets  of  zeros  is 
investigated  and  described  as  far  as  presently  possible. 

The  discussion  of  zeros  in  this  chapter  is  for  state  space  systems 
represented  in  the  frequency  domain.  As  was  discussed  in  Section  2.2,  when 
the  Laplace  transform  is  applied  to  a  continuous  system  (2.2.1)  and  a  z- 
transform  is  applied  to  a  discrete  system  (2.2.5),  the  functional  forms 
of  the  resulting  systems  are  the  same.  Furthermore,  if  the  transformed 
system  is  represented  by  a  system  matrix  there  is  no  loss  of  information. 
Since  the  theory  here  applies  equally  well  to  continuous  and  discrete 
systems,  state  space  systems  will  be  referred  to  by  their  system  matrix  to 
introduce  generality  into  the  notation. 

Section  2  provides  the  motivation  for  and  definition  of  a  set  of 
zeros  for  a  state  space  representation  called  invariant  zeros.  These  zeros 
are  then  characterized  by  a  rank  test  and  minors  of  the  system  matrix.  A 
slightly  different  definition  which  also  appears  in  the  literature  is  also 
discussed.  The  Smith  form  of  P(s)  is  the  main  analytical  tool  used  in 
developing  these  results. 

Section  3  explores  the  interrelationship  between  transmission 
zeros  and  invariant  zeros.  This  results  in  the  introduction  of  decoupling 
zeros.  Decoupling  zeros  are  shown  to  be  the  eigenvalues  of  the  unreachable 
and/or  unobservable  modes  of  the  system.  The  invariant  zeros  are  then  shown 
to  contain  all  of  the  transmission  zeros  and  some,  but  in  general  not  all. 


decoupling  zeros. 


41 


Section  4  introduces  and  defines  system  zeros.  Then  it  is  shown 
that  the  set  of  system  zeros  is  the  set  of  transmission  zeros  and  decoupling 
zeros,  and  contains  the  set  of  invariant  zeros  as  a  subset. 

An  example  is  carried  through  this  chapter  to  illustrate  these 
definitions . 

4.2.  Invariant  Zeros 

Zeros  defined  for  a  state  space  representation  can  be  characterized 
in  a  way  similar  to  the  way  transmission  zeros  are  related  to  frequencies 
whose  propagation  through  a  system  is  blocked.  Let  P(s)  be  a  (n+r) x  (n+m) 
system  matrix.  In  [17]  it  is  shown  that  a  system  excited  by  an  input  of 
the  form 


u(t)  -  g exp(zt) U  (t) 


(4.2.1) 


will  produce  an  identically  zero  output  if  and  only  if 


P<*> 


0 


(4.2.2) 


where  x  is  an  appropriate  initial  state.  Since  x  and  g  are  vectors, 

O  O  ' 

(4.2.2)  suggests  that  P(s)  must  lose  rank  at  s-z.  To  obtain  the  precise 
relationship,  the  Smith  form  Pg(s)  of  p(s)  is  used  to  investigate  the 
structure  of  P(s).  The  matrix  Pg(s)  is  given  by 

Pg(s)  =  diag[€1(s),62(s),...,€  (s),0,...,0]  (4.2.3) 

where  p  is  the  normal  rank  of  P(s)  (see  Section  2.3).  Again,  the  zero 
polynomial  z(s)  is  defined  as 
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z(s)€  »  n  (s) .  (4.2.4)* 

i»l  1 

The  zero  polynomial  can  be  faccored  into  linear  factors 


z(s)  «  (s-s1)(s-z,)- •  •  (s-zk)  (4.2.5) 

where  any  one  of  the  complex  numbers  z^,  i*l,...,k  may  be  repeated.  This 
gives  the  following  definition: 

Definition  5  [17] 

the  invariant  zeros  of  the  system  matrix  P(s)  are  the 
complex  numbers  z^,  i  k;  the  roots  of  the  zero  polynomial  z(s)  of 

?(s).  C 

The  name  invariant  zero  comes  from  the  fact  that  invariant  zeros 
are  the  roots  of  the  invariant  polynomials  of  P(s). 

Example  9  [  17] 

Lee  a  system  be  defined  by  the  matrices  (A.,B,C,D)  where 


r 


The  system  matrix 


0  0 
1  0 
0  3 

0  0 
0  0 
0  0 

0  0 
1  0 
0  1 

P(s)  is 


0  0  0 

0  0  0 

0  0  0 

-4  0  0 

0-10 
0  0  3 

1  0  0 

1  0  1 

0  0  1 

given  by 


0  -1 

-1  0 

1  -1 

0  0 

0  1 

-1  -1 


There  is  some  ambiguity  by  also  defining  this  polynomial  as  the 
zero  polynomial.  Hopefully,  this  will  be  justified  below. 


43 


”1 

0 

0 

0 

0 

0 

0 

1 

0 

3-1 

0 

0 

0 

0 

1 

0 

0 

0 

s  -3 

0 

0 

0 

-1 

1 

0 

0 

0 

s+4 

0 

0 

0 

0 

0 

0 

0 

0 

s+1 

0 

0 

-1 

0 

0 

0 

0 

0 

s-3 

1 

_  _  - 

1 

1 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

The  Smith  form  of  P(s),  S  (s)  can  be  calculated  as 

P 


Ps(s) 


I?  .  -  ° 

0  (s+1) (s-2) 

0  0 


The  invariant  zeros  are  z  «■ -1  and  z  =  2.  G 

Definition  5  and  properties  of  the  Smith  form  (see  (2.3,6)) 
yield  a  precise  statement  to  the  rank  characterization  of  invariant  zeros 
suggested  by  (4.2,2).  This  is  given  by  the  following  theorem. 

Theorem  4  [21] 

A  complex  number  z  is  an  invariant  zero  of  a  system  matrix  ?(s) 
if  and  only  if 

p[P(z)]<  Pa[P(  *)  ]  3 


The  proof  of  this  theorem  was  given  in  Section  2.3. 
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As  an  alternative  approach  for  defining  zeros  for  a  system 
matrix  P(s),  several  authors  have  defined  these  zeros  directly  from  the 
rank  condition  suggested  by  (4.2,2)  in  the  following  way: 

Definition  6  [22] 

A  complex  number  z  is  a  zero  of  a  system  matrix  P(s)  if 

p[P(z)]  <  n  +  min(m,r)  C3 

Davison  and  Wang  [19]  call  this  set  of  zeros  "transmission  zeros." 
As  will  be  shown  below,  this  set  of  zeros  is  not  the  same  as  those  zeros 
defined  by  Definition  1.  In  this  paper,  "transmission  zero"  shall  refer 
only  to  a  zero  defined  by  Definition  1. 

Note  that  if 

fl  fP(0]  *  n  +  min(r,m)  (4.2,6) 

n 

then  the  zeros  defined  by  Definition  6  are  the  same  as  invariant  zeros 
(Definition  5)  as  Theorem  4  shows.  However,  if 

0  [?(•)]  <  n  +  min(r,m)  (4.2,7) 

n 

then  any  complex  number  z  will  satisfy  Definition  6.  Hence,  if  (4.2,7) 
holds,  the  set  of  zeros  obtained  from  Definition  6  is  the  entire  complex 
plane.  In  general,  a  system  whose  system  matrix  P(s)  satisfies  (4.2,7) 
is  called  degenerate  [22]. 

Example  10 

Let  P(s)  be  defined  as  in  Example  9.  Since  the  Smith  form 
has  rank  8  so  does  ?(s),  i.e.  P(s)  has  full  normal  rank.  Consider  ?(-l). 
Inspection  shows  that  the  fifth  column  of  P(-l)  is  all  zeros.  Hence, 

0[P(*1)]  <  flnCP( ' )] 


Therefore,  by  Theorem  4,  z  =  -l  is  an  invariant  zero.  Similarly,  z 


can 
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be  checked.  Note  that  these  same  two  zeros  would  be  obtained  using 
Definition  6.  a 

When  transmission  zeros  were  introduced  by  Definition  1  the 
order  of  a  transmission  zero  was  also  defined  (3.2,10).  Similarly,  define 
the  algebraic  order  of  an  invariant  zero  as  follows:  Let  the  zero 
polynomial  z(s),  defined  by  (4.2,4),  be  factored  as 

h  k.  h 

z(s)  -  tr  (a-z.)  \  k  -  Z  k,  ,  k  >  0  (4.2,8) 

i-1  i-1 

where  the  complex  numbers  z^,  1*1,...,  h  are  required  to  be  distinct. 

Then  z^  is  said  to  be  an  invariant  zero  with  algebraic  order  (or  multi¬ 
plicity)  k^ .  Note  chat  Theorem  4  doesn't  permit  the  algebraic  order  of 
an  invariant  zero  to  be  calculated.  However,  for  state  space  systems, 
the  geometric  order  (or  multiplicity)  of  an  invariant  zero  is  defined  to 
be  the  rank  deficiency  of  P(z^)  [17].  Systems  characterized  by  equality 
of  these  two  multiplicities  for  all  distinct  invariant  zeros  are  said  to 
have  simple  structure.  If  these  two  multiplicities  are  different  for  any 
invariant  zero,  the  systems  are  said  to  have  non-simple  structure  [17]. 

Definition  5  and  Theorem  4  do  not  provide  very  tractable  methods 
for  computing  invariant  zeros,  even  for  systems  of  a  modist  order.  This 
motivates  the  following  theorem  which  provides  an  alternative  method  for 
calculating  zeros.  It's  proof  follows  directly  from  the  properties  of 
the  Smith  form,  Section  2.3. 

Theorem  5  [17] 

The  invariant  zeros  of  P(s)  are  the  complex  roots  of 
the  monic  greatest  common  divisor  of  all  non-zero  minors  of  maximum 
order  of  P(s) .  a 
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Note  Chat  che  algebraic  order  of  an  invariant  zero  can  be 
obtained  from  Theorem  5.  This  theorem  is  useful  when  m  =  r  and  the  system 
is  non-degenerate.  In  this  case  there  is  only  one  minor  of  maximum  order 
so  thac  the  invariant  zeros  are  the  complex  roots  of 

z(s)  =■  det[P(s)]  -  0  (4.2,9) 

provided  that  z(s)  has  been  adjusted  to  be  monic. 

The  property  expressed  by  this  theorem  is  sometimes  used  to 
define  invariant  zeros.  From  the  comments  on  the  Smith  form,  Section  2.3, 

i 

it  is  easy  to  see  that  Definition  5  would  be  a  direct  consequence  of  this 
new  definition.  Sometimes  the  zeros  of  P(s)  are  defined  to  be  the  greatest 
common  divisor  of  all  minors  of  maximum  order  of  P(s)  with  the  under- 

I 

standing  chat  if  che  minors  of  m aximim  are  all  zero  (the  system  is 
degenerate),  the  set  of  zeros  is  the  whole  complex  plane.  This  is  equiv¬ 
alent  to  Definition  6. 

Example  11 

Consider  P(s)  given  in  Example  9.  The  minors  of  order  8,  the 
maximum  order,  are  calculated  as: 


:>> 

a  - 

p(i::: 

-.3  . 

.,7,9’ 

s)  - 

p(!"' 
*  ' 

.,8 

.,6,8, 

9’s) 

p(i::: 

.,8 

•,5,7, 

3 , 9  ’ S 

-,8 

,4,6, . 

;9;s) 

p(1,3 ’ 
U,2, 

.,8 

3,  ;5 , . 

-,9;S 

(s+1) (s+4) (s-3)(s-2) 

4(s-l) (s+1) (s+4) (s-2) 
-  -  2 (s+1) (s+4) (s-2) 

■  0 

=■  4(3+1)  (s-1)  (s-2) 


u 

f  ^ 


u 


j 


lI 


i 
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p(i;2;4;8..i9;»)  -  2(s+4)(s+i)(s-2) 

P(J;2[4,’8..,9;s)  •  2(344)(s+1)(3-2) 

p/l» • • • »3  m  Q 
U,3,...,9'  ' 

-  4(s+4) (s+i) (s-2) 

*  >  •  •  •»* 

The  monic  greatest  common  divisor  of  the  non-zero  minors  is  easily  seen  as 

(s-2) (s+l) 


Hence,  the  invariant  zeros  are  z  -  2  and  z--l  which  checks  with  previous 
calculations.  O 


4.3.  Transmission  Zeros,  Decoupling  Zeros,  and  Invariant  Zeros 

It  is  natural  to  ask  if  invariant  zeros  are  related  to  trans¬ 
mission  zeros.  The  answer  is  given,  in  part,  by  the  following  theorem. 
Theorem  6  [13,  p.  Ill] 

For  reachable  and  observable  systems,  the  invariant  zeros  coin¬ 
cide  with  the  transmission  zeros.  a 

Theorem  6  is  simply  a  restatement  of  a  theorem  due  to 
Rosenbrock  which  states  that  the  invariant  factors  of  the  Smith  form  of 
P(s)  are  exactly  those  polynomials  which  appear  as  the  numerator  poly¬ 
nomials  of  the  Smith -McMillan  form  of  G(sl  provided  that  the  state  space 
system  which  gives  rise  to  P(s)  and  G(s)  is  reachable  and  observable. 

When  the  system  is  not  reachable  and/or  observable,  the 
invariant  zeros  are,  in  general,  a  larger  class  than  the  transmission 

zeros.  To  completely  characterize  all  invariant  zeros,  a  new  kind  of 
26ro  is  introduced  called  &  decoupling  zero. 


a  t~  in  It  Ion 


The  input  decoupling  zeros  are  the  comp  Lex  roots  of 


the  invariant  factors  of 


efinition  3  [  L3] 


The  output  decoupling  zeros  are  the  complex  roots  of  the 


invariant  factors  of 


'efinition  9  [  13] 


The  input-output  decoupling  zeros  are  those  decoupling  zeros 
that  are  both  input  decoupling  zeros  and  output  decoupling  zeros. 


The  following  theorem  is  Immediate  from  the  Smith  form  of  the 


matrices  in  Definitions  7  and  3.  Note  that 


because 


The  input  decoupling  zeros  are  the  such  that 


The  output  decoupling  zeros  are  the  z  €  &  such  that 


Thus ,  the  input  decoupling  zeros  are  the  eigenvalues  of  the 


unreachable  modes  and  the  output  decoupling  zeros  are  the  eigenvalues  of 


the  unobservable  modes  of  the  system 


Example  12 


Let  P(s)  be  the  system  matrix  in  Example  9.  By  Theorem  7 


z  -  -4  is  a  input  decoupling  zero  and  z  - -1  is  an  output  decoupling  zero.  a 


For  systems  which  are  not  reachable  and/or  observable  the  set 


of  invariant  zeros  contains  both  transmission  zeros  and  decoupling  zeros 


purpose 


which  of  these  zeros  appear  as  invariant  zeros 


First,  it  is  useful  to  investigate  why  decoupling  zeros  don't 


appear  among  the  transmission  zeros.  Of  course,  it  is  related  to  the 


fact  that  the  transfer  function  matrix  represents  only  the  reachable  and 


observable  parts  of  a  system.  This  can  be  made  precise  in  the  following 


Write  the  Smith  form  of  the  matrix  in  Definition  7  as 


where  L.(s)  and  R.(s)  are  unimodular  matrices.  Form  the  n  *  n  matrix 
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Now  (4.3,3)  can  be  written  as 

(sI-A.:  -B]  -  L1(3)Q1(s)[In0]R1(s) 

-  Q1(s)[«T^I:  -I]  (4.3,6) 

this  shows  that  Q^(s)  is  a  common  left  divisor  of  (sI-A)  and  (-B)  . 
Furthermore,  Q^(s)  will  contain  all  of  the  input  decoupling  zeros. 

Now  Che  matrix  Qj_(s)  is  invertable  by  construction  as  is  (sI-A).  Hence, 
when  the  transfer  function  matrix  of  the  original  system  is  formed,  it  is 
found  that 

G(s )  -  C(sI-A)"lB  +  D  -  C[(sI-A)'lQ^I(s)][Q1(s)B]  +  D 

■  C(sI-A)*1B  +  D  (4.3,7) 

This  shows  that  the  input  decoupling  zeros  "cancel  out"  when  the  transfer 
function  matrix  is  formed  so  that  they  do  not  appear  among  the  transmission 
zeros.  A  similar  analysis  can  be  carried  out  on  the  output  decoupling 
zeros.  They  represent  the  presence  of  a  common  right  divisor  between 
(si -A)  and  C. 

This  analysis  can  be  used  to  determine  which  transmission  zeros 
and  decoupling  zeros  appear  as  invariant  zeros.  Consider  a  system  matrix 
P(s).  Assume  that  r  £  m  and,  for  simplicity,  assume  chat  the  normal  rank 
of  P(s)  is  u  +  r .  First,  the  matrix  Q^(s)  is  obtained  from  (sI-A)  and 
(-B)  as  above  (see  (4 . 3  ,3) - (4 . 3 ,6) ) .  It  contains  the  input  decoupling 
zeros.  Then,  in  a  similar  fashion,  a  matrix  Q0(s)  is  obtained  from  (sI-A) 
and  C  so  that  Q^(s)  contains  the  remaining  output  decoupling  zeros; 
that  is,  Q0(s)  contains  the  output  decoupling  zeros  minus  the  input-output 
decoupling  zeros.  Now  P(s)  can  be  written  as 
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Q^s)  0 

(si -A)  -B 

Q2(s)  0 

0  1  r 

C  D 

0  I 

*  ,  V  —  ,  V  ^ 

* , 

-  Q*(s)  P(s)  Q*(s) 


(4.3,8) 


The  matrix  P(s)  can  be  thought  of  as  representing  a  new  system,  obtained 
from  the  old  system,  such  that  the  new  system  is  reachable  and  observable. 

The  minors  of  P(s)  are  described  using  the  formula  (2.3,7) 
introduced  in  Section  2.3, 


l,l2 ' ' 
J 1,J2 *  * 


E 

l^or,  <.  .  .<or  ^  n+r 
1  P 

l^r,  <. .  .<  Y  ^n+r 
1  P 


,1  -  a 
P 


(4.3,9) 


s) 


By  Theorem  5,  the  greatest  common  divisor  z(s)  of  all  non-zero  minors 

of  maximum  order  determines  the  invariant  zeros.  Set  p  *  n  +  r  to  calculate 

these  minors.  Note  that  in  this  sum,  every  minor  of  maximum  order  of  P(s) 

appears.  From  the  discussion  on  Smith  form  (Section  2.3),  it  follows 

that  the  invariant  polynomials  of  P(s)  will  appear  as  factors  of  z(s). 

Since  P(s)  is  reachable  and  observable,  by  Theorem  6  they  represent  the 

transmission  zeros.  Hence,  the  invariant  zeros  of  P(s)  contain  all  of  the 

transmission  zeros  of  the  system. 

★ 

Now  Qji(s)  has  only  one  minor  of  order  n  +  r,  so  it  appears  as 

a  factor  in  every  term  in  the  sum  (4.3,9)  and  as  a  factor  of  z(s).  Since 
★ 

Q^(s)  is  square,  its  minor  of  order  n  +  r  is  its  determinant.  From  (4.3,5) 
and  (4.3,8),  it  is  easily  seen  that 
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det[Q*(s).i  -  det[Q,(s)]  ■  or£  ^  (s)(3  2  (s)  . .  .P  n(s)  (4.3,10) 

where  or  is  a  non-zero  constant.  Therefore,  all  of  the  input  decoupling 
zeros  appear  as  invariant  zeros.  If  m  ■  r,  this  same  analysis  holds  for 
the  output  decoupling  zeros.  If  r  <  tn  this  is  no  longer  true,  i.e.  the 
invariant  zeros  will  contain  only  some  of  the  output  decoupling  zeros. 

A  concise  statement  of  these  results  is  given  in  the  following  theorem. 
Theorem  8  [23] 

Given  a  system  in  which  r«m,  the  set  of  invariant  zeros 
consists  of  the  set  of  transmission  zeros,  the  set  of  input-  decoupling 
zeros,  and  the  set  of  output  decoupling  zeros  minus  Che  input-output 
decoupling  zeros. 

If  r  <  m,  Che  set  of  invariant  zeros  consists  of  the  sec  of 
transmission  zeros,  the  set  of  input  decoupling  zeros,  and  possibly  some 
of  the  set  o’f  output  decoupling  zeros  minus  the  input-output  decoupling 
zeros . 

If  m  <  r,  the  set  of  invariant  zeros  consists  of  the  set  of  trans¬ 
mission  zeros,  Che  set  of  output  decoupling  zeros,  and  possibly  some  of 
the  set  of  input  decoupling  zeros  minus  the  input-output  decoupling 
zeros.  £ 

It  would  be  desirable  Co  characterize  those  decoupling  zeros 
that  do  not  appear  as  invariant  zeros.  One  characterization  is  given 
in  [23] . 

Example  13: 

Consider  the  system  matrix  P(s)  given  in  Example  9.  It  is 
shown  there  and  in  Example  12  chat  the  invariant  zeros  are  z  -  -l  and 

z  <■  2  ,  chat  z  -  -4  is  an  input  decoupling  zero,  and  that  z  •  -1  is  an 

output  decoupling  zero.  Note  that  m  *  2  <  3  *  r.  In  this  example,  Che 

input  decoupling  is  not  an  invariant  zero,  which  is  consistent  with 
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Theorem  8.  It  also  follows  from  Theorem  8  that  z  ■  2  Is  a  transmission 
zero.  This  can  be  varified  by  calculating  the  transfer  function  as 


G(s) 


0 

(s-3) (s-1) 
0 


_11' 

s-1 

1 

s-3 

2 

s-3 


and  its  Smith-McMillan  form,  M(s) ,  as 

1  0 

0  s-2 

0  0 

which,  by  Definition  1,  displays  a  transmission  zero  at  z  »  2.  □ 


M(s) 


1 


(s-3)  (s-1) 


4.4.  Systems  Zeros 

So  far  transmission  zeros,  decoupling  zeros,  and  invariant  zeros 
have  been  defined  for  a  system.  Since  in  general  the  invariant  zeros  don't 
include  all  of  the  decoupling  zeros,  it  is  desirable  to  provide  a  definition 
which  will  include  all  previously  defined  zeros  of  a  system. 

With  this  goal  in  mind,  consider  a  system  which  gives  rise  to  a 
system  matrix  P(s)  and  a  transfer  function  matrix  G(s).  The  minors  of  G(s) 
are  related  to  the  minors  of  P(s)  by  the  following  equation  [ 13 j . 


i i  1, . . . ,n, n+i .,..., n4i 

G(.  ,P;s)  -  P(  1  P;s)/det[sI-A]  (4.4,1) 

To  see  this  relationship,  first  consider  P(s).  It  can  be  factored  as: 


P(s)  - 


si -A  -B 

si -A  0 

J 

a 

o 

i 

In  -(sI-A)  "1B  ~J 

C  D 

0  lr 

C  lr 

0  C (sI-A) ”^B  +  D  | 

— 

r 

L_  — 1 

(4.4,2) 


' 
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Taking  determinants  of  each  side  yields: 

det[p(s)]  -  det[sl -A]det[C (sI-A) *  +  D]  (4.4,3) 


Let  c^  denote  a  row  of  C,  denote  a  column  of  B,  and  d^,  an  element 
of  D.  It  follows  from  (4.4,3)  that 


.  1 , • • • ,n,n+i 

StjOO  ■  c^Cal-A)'  Oj  +  dtj  -  P(1 . a>n+J)*>M«C.X-*]  <4.4,4) 

The  element  g^(s)  can  be  thought  of  as  a  1  X  1  minor  of  G(s).  In 
exactly  the  same  way,  this  argument  can  be  extended  to  any  minor  of 
G(s)  of  any  order,  resulting  in  (4.4,1). 

Equation  (4.4,1)  suggests  the  following  definition: 

Definition  10  [24] 

The  system  zeros  of  a  system  matrix  P(s)  are  the  complex  roots 
of  the  monic  greatest  common  divisor  of  all  non-zero  minors  of  P(s)  of 
the  form 

1,^, • . • ,n, n+i « , n . i «,..., n ! i 
PC  12  P . s \ 

U,2,...,n,n+j1,n+j2,...,n+jp’  ’ 


where  n  +  p  is  the  maximum  order  of  a  minor  of  this  form  which  is 
non-zero.  C 


Definition  10  and  Theorem  5  are  similar  in  that  both  involve 
minors  of  maximum  order  of  P(s).  However,  Theorem  5  requires  that  all 
minors  of  maximum  order  be  calculated,  where  as  Definition  10  is 
restricted  to  the  subset  of  those  minors  that  contain  the  first  n  rows 
and  n  columns  of  P(s).  When  the  number  of  inputs  and  the  number  of 
outputs  is  the  same  and  the  system  is  non-degenerate,  then  the  invariant 
zeros  are  the  same  as  the  system  zeros  because  there  is  only  one  minor 
of  maximum  order.  However,  when  one  of  these  conditions  fails  to  hold, 
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it  would  suggest  that  Definition  10  defines  a  larger  set  of  zeros  than  the 
invariant  zeros.  The  exact  relationship  is  given  by  the  following  theorem: 
Theorem  9  [24] 

a)  The  set  of  system  zeros  of  a  system  matrix  P(s)  consists  of: 

i)  the  set  of  transmission  zeros 

ii)  the  set  of  input  decoupling  zeros 

iii)  the  set  of  output  decoupling  zeros  minus  the  set  of 
input-output  decoupling  zeros. 

b)  The  set  of  invariant  zeros  is  a  subset  of  the  system  zeros. 

Proof :  This  proof  is  in  the  same  spirit  as  the  proof  of  Theorem  8. 

Consider  a  system  matrix  P(s)  with  normal  rank  n  +  k.  The  following 
result  will  be  needed  below:  If  n  +  p  is  the  maximal  order  of  a  minor  of 
the  form  in  Definition  10,  then  p  *  k.  Note  that  this  also  proves  part  (b) 
of  the  theorem,  by  Theorem  5.  The  normal  rank  is  the  order  of  the 
largest  non-zero  minor.  To  find  one  such  minor,  examine  each  column  in 
the  order  l,2,...,n+m  and  select  each  column  which  is  linearly  independent 
from  those  already  selected.  Since  si -A  is  non-singular,  the  first  n 
columns  qualify.  Exactly  k  of  the  remaining  columns  qualify.  Repeat  this 
procedure  for  the  rows  of  P(s).  Again  the  first  n  qualify  with  exactly  k 
remaining  rows.  This  minor  has  the  form  required  by  Definition  10. 

Hence,  p  *  k. 

Return  to  the  system  matrix  P(s).  It  can  be  written  as  a  product 
of  matrices  as  in  (4.3,3): 
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QL  Cs)  ol 


PCs) 


Q2(s) 


n 


?(s) 


0?Cs)PCs)Q*Cs) 


Ci.4,5) 


Recall  that  Q^CS)  contains  the  input  decoupling  zeros,  Q2(s)  contains 
the  output  decoupling  zeros  minus  the  input-output  decoupling  zeros,  and 
PCs)  represents  a  reachable  and  observable  system.  To  calculate  the 
system  zeros,  the  formula  given  in  (4.3,9)  is  used,  where  the  minors  are 
required  to  be  of  the  type  in  Definition  10. 


?c 


L, 


,  n ,  n+  i  ^  , 


,o+i 


1 ( • * ■ ,  ^ » n i  J  ^ , • . • , a  1  j 


;s) 


^  1 , . . . ,n, n i i ^ , 


^  l'  *  ’  ,<Qfn+psn+r 

UY,<. .  .<y  scrfm 
1  tifp 


Ql^ 

j.  j  .  ,  •  • . 


nfp 


,rH-p 


;s)  . 


_  y  i » • 

PC 


,a 


n-f*p 


*,  1  nfp 


Y  Y 


=■>  (4-4’5> 

■L  P 


where  n  -r  p  is  the  normal  rank  of  PCs).  In  this  sum  (4.4,6),  consider  a 

term  where  there  exists  an  3  such  that  3  >  n  and  3  4  n  +  i„  for  all  i, 

s  s  s  l 

★ 

1  S  i  s  p.  Then  that  minor  of  Q^(s)  is  zero,  as  can  be  seen  from  the  special 

*  * 
form  of  Q^(s)  in  (4.4,5).  The  same  argument  holds  for  Q7(s).  This  implies 

chat  the  sum  in  (4.4,6)  can  be  replaced  by 


1 . .  . .  ,n,rH-i  , 

P(  1 

1.. .. ,n, nt  j  ^ , 


,nfip 

P 


s) 


_  I,. 

<i«eCQl(s)]P(1 


,  n ,  irt-i  ^ 


,rt+-i 

,nfJP;s)d«cCQ2 


(s)] 

C*,4,7) 


Since  P(s)  represents  a  reachable  and  observable  system,  by  Theorem  6 

the  invariant  factors  of  P(s)  are  just  the  invariant  factors  of  its 

transfer  tunction.  Therefore,  ?(s)  can  have  at  most  p  £  min  (m,r)  invariant 

factors  that  are  not  identically  1.  It  follows  that  the  Smith  form  ?  (s) 

s 
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of  P(s)  is 

r  i  o 

i  a 

Ps(s)  =  |  (4.4,8) 

j_0  diag[6^(s)  »•••»€  (s)  ,0 , ...  ,0__ 

Using  Che  determental  divisor  properties  of  the  Smith  form  (Section  2.3) 
and  (4.4,7)  and  (4.4,8),  it  is  now  seen  chat  the  greatest 
common  divisor  of  Che  minors  in  (4.4,7)  is 

det[Q1(s)]r€1(s)...€p(s)]detrQ2(s)]  (4.4,9) 

When  this  polynomial  is  adjusted  to  be  monic,  its  roots  form  a  complete 
set  of  decoupling  zeros  (by  construction  of  Q^(s)  and  Q?(s))  and  all  the 
transmission  zeros  of  P(s).  This  completes  the  proof.  □ 

Example  14 

Consider  the  system  matrix  ?(s)  introduced  in  Example  9.  In 
Example  11,  all  minors  of  maximum  order  were  calculated  to  find  che 
invariant  zeros.  Using  Definition  10,  Che  minors  of  interest  are 

p(i,,!!!J,8,a>  -  -  7  (s+d(s+4)(s-3Hs-2) 

p([;;;;;68j8,9,s)  -  4(*-dc*+i)  (**5(3-2) 

The  monic  greatest  common  divisor  of  these  minors  is 

(s+1)  (s+4)  (s-2) 

so  that  che  system  zeros  are 

z=-l,z=-4,z=2.  S 
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Definition  10  can  be  used  to  interpret  system  zeros  in  another 
way.  Consider 


1,2, . . . ,n, i ., . .  .  ,  i 

^  ( i  2  n  i  i  > s  ^ 

li* i • • • >“>J  •  >Jp 


(4.4,10) 


which  is  a  minor  of  P(s)  that  satisfies  Definition  10.  This  minor  describes 
a  state  space  representation 


x  »  Ax  +  B  u 
P 

y  =  C  x  +  D  u 
P  P 


(4.4,11) 


where  the  matrices  B  ,  C  ,  and  D  are  obtained  from  the  rows  and  columns 

P  P  P 

indicated  in  (4.4,10);  the  input  and  output  spaces  now  both  having 
dimension  p.  This  system  is  non-degenerate  as  (4.4,10)  is  required  to  be 
non-zero.  Hence,  by  Theorem  5,  the  invariant  zeros  of  (4.4,11)  are  the 
roots  of  (4.4,10).  In  this  light,  Definition  10  says  that  the  system 
zeros  of  P(s)  are  those  invariant  zeros  which  are  common  to  all  non-degenerate 
subsystems  formed  by  selecting  p  rows  of  the  C  matrix  and  p  columns  of  the 
B  matrix;  i.e.  subsystems  of  the  form  (4.4,11). 


59 


CHAPTER  5 

PROPERTIES  OF  ZEROS 

5.1.  Introduction 

Chapter  4  introduced  the  definitions  of  several  types  of  zeros 
for  state  space  systems  and  described  their  interrelationship.  Various 
types  of  zeros  can  be  calculated  from  the  corresponding  definitions  or 
theorems  of  the  last  chapter.  However,  these  definitions  and  theorems, 
which  are  based  primarily  on  the  Smith  form,  are  not  readily  adapted  to  a 
digital  computer  and  hand  calculation  is  tedious  for  higher  order  systems. 
Therefore,  an  important  task,  is  to  present  several  algorithms  for  calculat¬ 
ing  zeros,  based  on  additional  properties  of  zeros,  that  are  easier  to 
apply  than  using  the  basic  properties,  or  either,  that  can  be  readily 
adapted  to  a  digital  computer. 

This  chapter  is  devoted  to  the  development  of  those  properties 
which  will  be  useful  in  constructing  such  algorithms.  It  should  be  noted 
chat  these  properties  are  important  in  their  own  right.  In  fact,  many  of 
these  properties  are  just  generalizations  of  properties  of  zeros  of  sca.ar 
transfer  functions .  In  particular,  it  is  shown  below  that  invariant  zeros 
are  unchanged,  bv  state  feedback,  that  they  are  the  limiting  positions  of  t-.a 
systems  poles  under  high  gain  feedback,  and  that  tr.ey  are  the  on  1*3  of  an 
inverse  system. 

This  chapter  also  examines  tne  geometrical  properties  of  zeros 
Although  these  properties  ire  a  consequence  of  ter  inruns  gi  ran  ear._er 
the  invariant  zeros  and  transmission  zeros  can  oe  defined  ann  i eveluued 
from  a  geometrical  analysis  ci  tne  state  space  r  sores  anna  r.urc  1.2  _  -.r 
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(2.2,5)  of  a  system.  It  turns  out  that  the  invariant  zeros  are  intimately 
connected  to  two  well  known  subspaces,  the  largest  (A, B) -invariant  subspace 
in  the  null  space  of  C  and  the  largest  controllability  subspace  in  the  null 
space  of  C  [12].  The  definitions  of  these  subspaces  are  generalized  to 
systems  where  Dj^O  and  the  connection  with  invariant  zeros  is  made  explicit. 
These  geometrical  ideas  also  form  the  basis  of  a  new  algorithm  (presented 
in  Chapter  7  below)  for  calculating  these  two  subspaces,  and,  hence,  the 
invariant  zeros  of  the  system. 

The  geometrical  analysis  given  here  applies  equally  well  to 
continuous  systems  (2.2,1)  or  discrete  systems  (2.2,5).  Therefore,  to  in¬ 
troduce  generality  in  the  notation,  either  of  these  two  types  of  systems 
will  be  denoted  (A,B,C,D).  Some  of  the  discussion  of  the  geometrical  proper¬ 
ties  is  done  in  terms  of  discrete  systems.  This  is  for  pedagogical  reasons 
(as  will  be  clear  below).  The  analysis,  however,  applies  just  as  well  to 
continuous  systems. 

Section  2  discusses  in  decail  the  geometrical  ideas  needed  for  the 
discussion  of  invariant  zeros.  Two  important  subspaces,  the  unknown-input 
unobservable  subspace  and  the  null-output  reachable  subspace,  are  introduced 
and  their  relevant  properties  are  presented.  It  is  shown  that  a  system 
(A,B,C,D)  can  be  placed  in  a  canonical  form  to  exhibit  the  properties  of 
these  subspaces.  This  analysis  will  be  used  again  in  Chapter  7  below. 

The  specific  relationship  of  invariant  zeros  to  the  geometric  pro¬ 
perties  of  the  system  is  given  in  Section  3.  This  is  accomplished  via  the 
canonic  form  introduced  in  Section  2.  Several  other  geometrical  properties 
of  invariant  zeros  is  3lso  discussed. 
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Section  4  describes  various  properties  of  system  zeros  and 
invariant  zeros  which  will  be  useful  in  discussing  algorithms  for  calcu¬ 
lating  zeros.  These  properties  include  the  invariance  of  system  zeros 
under  the  coordinate  changes  of  the  input  space,  state  space,  and  output 
space  and  under  output  feedback.  The  invariant  zeros  are  shown  to  be 
unchanged  by  state  feedback  and  to  be  the  limiting  positions  of  the  systems 
poles  under  high  gain  feedback. 

Section  5  reviews  the  notion  of  an  inverse  system  and  develops 
the  basic  theory  behind  this  idea.  Such  a  system  is  then  constructed  and 
it  is  shown  that  the  poles  of  an  inverse  system  are  the  invariant  zeros 
of  the  original  system.  The  geometrical  ideas  of  Section  5.2  play  a  key 
role  here. 

5.2.  Some  Geometrical  Aspects  of  Linear  Systems 

The  definition  and  properties  of  two  subspaces  relevant  to 
the  discussion  of  zeros  are  presented  in  this  section.  These  subspaces 
were  first  introduced  in  [25]  and  the  notation  and  definitions  here 
follow  that  references. 

The  first  subspace  of  interest  is  the  unknown- input  unobservable 
subspace,  denoted  by  £  throughout  this  paper. 

Definition  11  [25] 

The  vector  £  is  an  element  of  the  ith  unknown- input  unobserv- 
able  subspace,  if  there  exists  a  control  sequence 


(uo*ui . Vi* 


02 


such  that  if  x  •?  then 
o 


y 


0 


•  y 


l-l 


-  0 


By  definition,  *%  the  whole  state  space.  □ 

Unknown- input  unobservable  subspaces  have  several  interesting 
properties  [25] .  Two  basic  properties  are  given  by  the  following  lemma; 
Lemma  l  [25] 

*>  £i  3£wi  fot  *u  1 

b)  For  some  kin,  ^  *  tor  J  “  » •  •  • 

Proof:  Part  (a)  If  +  then  there  exists  a  control  sequence 

lu0,ul,u2,  . . .u^)  (5.2,1) 

such  that 


yo  "yi "  “yi*°  (5-2*2) 

But  then,  £  for  Che  *am«  control,  sequence  (5.2,1)  by  simply  dropping 
u ^  from  the  sequence. 

Part  (b).  Let 


d(£l) 


(5.2,3) 


Then 


T0«n  and 


(5.2,4> 


as  can  be  seen  from  Part  (4).  Let  k  be  the  first  i  such  that 


Clearly  It  must  be  shown  that 


The  argument  is  by  induction.  The  hypothesis  (5.2,7) 


for  j  ■  L  is  shown  in  (5.2,6).  Suppose  that  (5.2,7)  holds  for  some  X 


7)  holds  for  y  ■  1 


X  for  some  control  sequence 


By  induction  hypothesis 


This  implies  there  exist  a  control  sequence 


But  Chen  the  application  of  the  control  sequence 
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(Uo,ui»---«uk+(x-1)'  uk+(i-l))  (5.2,13) 

where  the  Last  term  in  (5.2,13)  is  the  last  turn  in  (5.2,8),  shows  that 

Therefore,  (5.2,7)  follows.  G 

This  shows  that  after  no  more  than  n  steps,  the  unknown-input 
unobservable  subspaces  are  all  the  same.  In  the  sequel,  the  term  unknown- 
input  unobservable  subspace  shall  refer  to  this  subspace. 

The  subspaces  i  =0,1,2,...  can  be  characterized  in  terms  of 
the  matrices  (A,B,C,D).  For  this  purpose,  consider  the  following  algorithm: 
Algorithm  1  [26] 

Step  1:  For  i  ■  0,  set  Mq  ■  0  where  Mq  is  a  Un  matrix. 

Step  2:  Form  the  matrix 


r 


D 


(5.2,14) 


and  carry  out  the  indicated  multiplication.  Then 


find  a  non-singular  matrix  such  that 


Vi 


i+L 

0 
’  m 


Vi 


i+l_ 

n 


(5.2,15) 


where  F  has  full  row  rank. 


Step  3:  Set  i = i+1  and  go  to  step  2. 

The  following  theorem  relates  the  subspaces  x ^  to  Algorithm  1: 


l 
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Theorem  LO  [26] 


Let  be  the  matrices  obtained  from  Algorithm  1.  Then 


£ia^Ml.];  i  *0,1,2,.., 


Proof:  Consider  any  §  =*x^  and  any  control  sequence 


(u0,u1,...>ui_L) 


(5.2,16) 


applied  to  a  system  (2.2,5).  It  follows  directly  from  the  repeated  appli¬ 
cation  of  (2.2,5)  that 


0Xi  Mcr4  M0A  1 

,  ,  cal'l  cal-2b 


CA  CB 

C  D  ‘  . 


MqAB  MqB 
CB  D 


0  0  u.  , 

L-  1 


(5.2,17) 


for  any  MQ.  Notice  that  the  blocks  in  the  first  two  rows  of  (5.2,17)  can 
be  associated  with  Fg  (5.2.14)  in  Algorithm  I.  Multiply  the  first  two  rows 
by  Sn((5.2,14))  and  use  (5.2,15).  Now  (5.2,17)  has  the  form 


o  Yi  Gr 


t-1 

_  ,i-2_ 
GjA  B 

•  •  • 

GjB 

i-2 

M  AI_3B 

MlB 

-2 

•  •  • 

D 

o  0 


!  U.  ,  . 

!  i-1 


(5.2,18) 
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It  is  seen  that  the  blocks  in  the  second  and  third  row  are 
associated  with  (5.2,14).  Again,  multiply  by  (5.2,15)  and  use 
(5.2,15).  Continue  this  process  for  i  steps.  By  the  non- singularity  of 
the  matrices  and  the  fact  that  =  0,  it  follows  that 

y0*yl~  **•  *yi-l“°  (5.2,19) 


if  and  only  if 


g/-1 

i-2 

GjA  ...  G^B 

F1 

ar 

* 

G2Ai_2  . 

•  y  ■  •  f2 

0 

uo 

Gi 

Fi 

. 

• 

Mi 

0 . 

.  0 

Vi 

— 

_ 

L-  _! 

(5.2,20) 


If  §  is  an  element  of  £  ,  then  by  Definition  11,  (5.2,19)  holds  for  some 
control  sequence  (5.2,16).  For  that  control  sequence  it  follows  from 
(5.2,20)  that 


M.C  a  0  (5.2,21) 

Conversely,  assume  chat  (5.2,21)  holds.  Then  since  the  matrices  F^  have 
full  row  rank,  there  exists  some  control  sequence  (5,2,16)  such  that  (5.2,20) 
is  satisfied.  This  implies  that  (5.2,19)  holds  and  by  Definition  11, 

Thus  the  theorem  is  proved.  23 

Algorithm  1  immediately  yields  the  first  property  of  the  sub¬ 
spaces  which  is  fundamental  to  what  follows; 


Lemma  2  [26] 

The  subspace  is  unaffecced  by  state  feedback. 

Proof:  Consider  a  system  defined  by  the  matrices  (A+BK , B , C+DK ,  D ) ,  from 
(5.2,15)  it  is  seen  that 


(5.2,22) 


Thus,  by  Theorem  10  the  subspace  £  are  the  same  for  the  system  (A,B,C,D), 
with  or  without  feedback.  a 

That  the  subspaces  are  not  changed  by  state  feedback  suggests 
that  they  might  be  very  powerful  for  describing  the  structure  of  a  system. 
Indeed,  its  first  important  consequence  is  embodied  in  the  following 
theorem. 

Theorem  11  [26] 

Consider  state  feedback  of  the  form 


Then  the  system 


A+BK  B  x 

i 

C+DK*  D  I  v 

-  J  LA 


after  the  appropriate  similarity  transfo 


form: 


rmations ,  takes  on  the  canonical 


where 


t  0 


S _ 

C1  0 


Bj_  0 


B3  B4 
D1  0 


and  7‘fF.+1J  = 


\ 

V. 

1 


Proof:  From  (5.2,15) 


(5.2,24) 


VA*+™  +  W 


Cx  +  Du 


M  .  ,x 
n+1 


(5.2,25) 


for  all  (x, u) .  Note  that 


by  hemma  1  and  Theorem  10.  Consider  rhe  feedback 
Substituting  (5.2,23)  into  (5.2,25),  lt  ls  s„„  th 


(5.2,26) 


defined  in  (5.2,23). 


is  seen  that  for  all  vCXfF  1 

^  n4. 1  * 
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and  x  €  £n 


Mn(A.x+  Bu) 
Cx  +  Du 


0 


(5.2,27) 


(By  Lemma  2,  this  substitution  doesn't  changed  ).  It  follows  that 

n 

(Ax+Bu)  is  an  element  of  This  leads  to  the  representation  (5.2,24) 

in  the  appropriate  bases.  Q 

The  canonic  form  (5.2,24)  shows  that  £^  is  (A, B)- invariant  (See 
Section  2.1).  In  fact,  £n  is  a  generalization  of  Wonham's  maximal  (A,B)- 
invariant  subspace  in  the  null  space  of  C  [12].  The  relationship  is  made 
precise  by  the  following  definition  and  theorem. 

Definition  12  [26] 

A  subspace  7  is  a  null-output  (A, B)- invariant  subspace  if  for 
every  x€^  there  exists  some  u  such  that  (Ax  +  Bu)  and  (Cx  +  Du)  *  0. 
Theorem  12  [26] 

A  maximal  null-output  (A, B)- invariant  subspace,  denoted!  , 
exists,  and  £n » £  where  £  is  (A  +  BK.  ) -invariant. 

Proof:  From  Definitions  11  and  12,  it  is  seen  that  Vr  £  for  every  « 

n  * 

a  null-output  (A, B) -invariant  subspace.  Conversely,  take  u=K  x  and 

consider  x££  .  From  canonic  form  (5.2,24),  it  is  obvious  that  £  satis- 
u  n 

fies  Definition  12.  Q 

ie 

The  subspace  £  has  been  very  useful  in  geometric  control  theory 
[12].  It  will  play  a  fundamental  role  in  what  follows. 

The  second  subspace  of  interest  is  the  null-output  reachable 

space. 
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Definition  13  [25] 

The  vector  ?  is  an  element  of  the  ith  null-output  reachable 

space,  if  there  exists  a  control  sequence 


(u0’ul’ ‘ ' * ,Ui-l) 


such  that  if  Xq  *  0  then  x^  *  |  and 


Q 


Like  2^,  has  nesting  properties  and  a  limit  condition  which 
are  given  in  the  next  lemma.  The  proof  is  similar  to  Lemma  1  and  is 
given  in  [22]  . 

Lemma  3  [25] 

a)  RiC*i+l  i  *0,1,2,... 

b)  For  some  ksn,  *  l\+ j  f°r  j  Q 

The  null-output  reachable  subspaces  are  related  to  another 
fundamental  subspace,  the  maximal  reachability  subspace  in  the  null  space 
of  C.  To  establish  this  connection,  define  the  reachable  subspace  >1  of 
the  pair  (A,B)  as 


«-<A|«[B]  >-»[8]  +Aft[B]  +  ...  +An'1ft[B]  (5.2,28) 


if  dfi.)  *n.  If  it  is  said  that  (A, B)  is  reachable.  Certain  sub¬ 

spaces  of  the  reachable  space  turn  out  to  be  of  interest  in  geometric 
control  theory.  In  general,  a  subspace  >1  is  a  reachability  subspace  it 
there  exist  maps  K:  U  and  J:  such  that 
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a  =  <(A+BK)  lft[BJ]>  (5.2,29) 

The  reachable  subspace  is  a  reachability  subspace  as  can  be  seen  by  taking 
K  =  0  and  J  »  I  .  Note  that  reachability  subspaces  are  (A, B) -invariant 
subspaces  by  (5.2,28)  and  the  definition  of  (A, B) -invariant  subspaces 
(Section  2.1).  The  particular  reachability  subspace  of  interest  here  is 
defined  next. 

Definition  14  [26] 

The  subspace  ft  is  a  null-output  reachability  subspace  if  there 
exist  matrices  K  and  J  such  that 

=  <(A+BK)  ift[BJ]> 

(C-t-DK)R  =  0 


Dft[J]  =0  □ 

The  relationship  between  2.  ,  &  ,  and  the  null-output  reachability 

n  n 

space  is  given  by  the  following  theorem. 

Theorem  13  [26] 

"k 

A  maximal  null-output  reachability  space,  denoted  si  ,  exists  and 

ft*  -  ftn  n  £ 

n  a 

It  is  generated  by  K  (5.2,23)  and  J  where 

*'Ji  ■ 

Proof:  Consider  any  null-output  reachability  subspace  ft.  It  follows  from 

Definitions  13  and  14  and  Lemma  2  that  Also,  by  reachability, 


■ . . ■ 
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(A+BK)ftc  ft,  and  (C+DK)ft  -  0  imply  that  ftc£*  by  Definition  12  and  Theorem 
12.  It  follows  that 


ftC*nn£n  (5.2,30) 

If  it  can  be  shown  that  ft^fl  &n  also  satisfies  Definition  14,  then  (5.2,30) 

shows  the  ftnfl  £n  is  the  required  maximal  null-output  reachability  space 

by  part  (b)  of  Lemma's  1  and  2.  First,  ft  ft  2  is  characterized  in  terms 

n  n 

of  canonic  form  (5.2,24): 

<A4,B4>>  (5.2,31) 

Note  that  the  right  side  is  a  subset  of  the  left,  as  can  be  seen  directly 

from  (5.2,24).  Now  consider  any  x.  ,€ft  H  £  .  From  (5.2,24)  it  follows 

J.+  in  n 

that 

*i+l  3  0  and  ^  1  0  (5.2,32) 

With  reference  to  (5.2,24),  it  is  seen  that 


ft  nx 

n  n  l 


(5.2,33) 


This  in  turn  implies  that 


x,  =  0  and  v,  =  0 
i  i 


(5.2,34) 


This  argument  can  be  used  to  show  that 


I 

I 

I 

I 

I 

I 

1 

I 

T 

i 

I 

I 

i 

I 

1 

I 

I 
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x j  *  0 ,  v j  —  0  j  — 


Hence,  it  has  been  shown  that 


a 

n 


V  <A4Ib4>} 


(5.2,35) 


(5.2,36) 


which  establishes  (5.2,31).  Finally,  from  the  construction  of  canonic 
form  (5.2,24)  it  follows  that  (5.2,31)  satisfies  Definition  14  for 


K  -  K*  and  ft[J] 


(5.2,37) 

Q 


Since  ft  is  a  (A, B) -invariant  subspace  for  K  (Definition  12),  it 
follows  that  canonic  form  (5.2,24)  would  admit  an  even  more  explicit  repre¬ 
sentation  in  the  appropriate  state  and  control  bases. 

Corollary  1  [27] 

Canonic  form  (5.2,24),  in  the  appropriate  state  and  control 
bases  takes  on  the  form 


xi+l 

A1 

0 

0 

B1 

0 

A 

xi+l 

= 

*3 

A5 

A6 

*3 

B5 

xi+l 

a3 

0 

V 

§3 

0 

_/i_ 

H 

CJ 

1 

0 

0 

D1 

0 

(5.2,38) 


; 


where  (A^,B^)  is  reachable.  □ 

Interperted  geometrically,  canonic  form  (5.2,33)  exhibits  a  direct 


sum  decomposition  of  state  space 


74 


^  0K/  sm,  4 

z  » 


(5.2,39) 


where 


Z  -  R 


X9%  =  £' 


—  /«  y< 

/w  is  isomorphic  to  /£ 


* 

~  £  * 
£  is  isomorphic  to  /ft 


(5.2,40) 


The  invariant  zeros  can  now  be  described  in  terms  of  canonic  form  (5.2,38) 
and  the  stace  space  decomposition  (5.2,40).  This  is  undertaken  in  the  next 
section. 


5.3.  Geometric  Properties  of  Invariant  Zeros 

Recall  that  the  invariant  zeros  are  the  roots  of  the  invariant 
polynomials  of  the  system  matrix  P(s)  (Definition  5).  The  invariant  zeros 
are  closely  related  to  canonic  form  (5.2,38)  developed  in  Section  5.2.  This 
connection  is  formally  established  by  the  methods  of  Chapter  4  applied  to 
the  results  of  the  geometric  analysis  in  Section  5.2. 


The  following  lemma  will  be  needed  in  the  proof  of  the  main  theorem 


below. 


Lemma  4  [28] 


Consider  a  system  with  system  matrix  P(s).  If 


SP^  •  '  MM 


I 

I 

I 

I 

I 

I 

I 

1 

i 

I 

I 

I 

I 

1 

1 

l 

I 

I 

I 
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where  T  is  defined  in  Algorithm  1  (Section  5.2),  then  there  exists  a 
n 

polynomial  matrix  W(s)  such  that 

W(s)P(s)  -  I 


This  implies  that 


P(s)  * 


where  in  this  work  means  "has  the  same  Smith  form  as".  □ 

Proof:  See  Reference  [28].  From  (5.2,15)  it  is  seen  that 

7t[rn+1]  *  0  implies  that 

’  °  -£n+l  =**  <5‘3’L) 

from  Theorem  10  and  12,  and 

^[Fn+1]  -  0  (5.3,2) 

By  Lemma  4,  such  a  system  has  no  invariant  zeros.  The  idea  of  the  following 
analysis  is  to  identify  a  subsystem  of  a  given  system  with  these  properties. 

Given  a  system,  bring  it  into  canonic  form  (5.2,38).  As  was 
discussed  in  Section  5.2,  this  is  done  by  applying  feedback  of  the  form 

u  -  K*x  +  v  (5.3,3) 

and  transforming  the  state  and-  input  spaces  (Corollary  1).  Let  these 
transformations  be  given  by 


I 
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2  *  Tx  and  u  *  Gv 


(5.3,4) 


where  T  and  G  are  non-singular  matrices  of  the  appropriate  size.  The 
transformed  system  has  a  system  matrix  P(s)  that  is  given  by 


sI-Al 

0 

0 

-Bi 

1 

O 

-a3 

Sl-A5 

-A6 

A 

-b3 

-B5 

-A3 

0 

sI-A? 

-i3 

0 

C1 

0 

0 

D1 

0 

— 

(5.3,5) 


The  transformed  system  is  obtained  by  substituting  (5.3,3)  and  (5.3,4)  into 
the  original  state  equations.  This  gives  the  following  relationship  between 
P(s)  and  P(s): 


T 

0 

si -A 

-B 

1 

1 1 
i  n 

0 

t‘l  0 

0 

rjL 

C 

1 

D  ! 

i  * 

1  K 

I 

m 

0  G 

3  P(s) 


(5.3,6) 


This  can  be  checked  by  straight  forward  calculation.  The  matrices  pre-and 
post-multiplying  P(s)  can  be  interpreted  as  elementary  row  and  column 
operations.  It  follows  from  the  comments  is  Section  2 . 3  on  the  Smith  form 
that 

P(s)*P(s)  (5.3,7) 

By  using  elementary  row  and  column  operations,  (5.3,5)  can  be  written  as 
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(5.3,8) 


When  Algorithm  1  is  applied  to  the  subsystem  (A^,B^,C^,D^) ,  it  is  found 
that 


(5.3,9) 


This  follows  from  the  construction  of  canonic  form  (5.2,24).  (See  proof  of 
Theorem  11).  It  also  follows  from  the  discussion  of  decoupling  zeros  in 
Section  4.3  and  the  fact  that  (Ag,B^)  is  reachable,  that 

[sI-A5  :  -B5]  *  [I  0]  (5.3,10) 


Hence,  by  Lemma  4  and  equation  (5.3,10),  it  follows  that 


»  P(s) 


(5.3,11) 


Now,  it  is  clear  from  Definition  5  that  the  following  theorem  has  been  proved. 
Theorem  14  [28] 

The  invariant  zeros  are  the  roots  of  the  invariant  polynomials  of 
(sI'Ay)  in  canonic  form  (5.3,8).  G 


The  proof  of  Theorem  14  relied  on  the  system  matrix  characteriza¬ 


tion  of  invariant  zeros  given  in  Definition  5.  However,  a  description  of 
invariant  zeros  can  be  formulated  from  a  purely  geometrical  point  of  view 
as  follows  [29].  Let  R  be  a  feedback  matrix  such  that 

(A+BK)£*c  S*  (5.3,12) 

Z  * 

Denote  by  A.^  the  map  (A+BK)  induced  in  /ft  .  The  invariant  zeros  are 

* 

y  * 

the  roots  of  the  invariant  factors  of  the  map  A^  restricted  to  /ft  . 
(Compare  to  the  remark  at  the  end  of  Section  5.2). 

There  is  also  a  geometric  connection  with  the  system  matrix. 
Consider  a  system  matrix  P(s)  with  full  normal  rank.  By  Theorem  4,  z  is  an 
invariant  zero  of  P(s)  if 


P[P(z)]  <  •3n[P(-)l 


This  implies  that 


rd P(Z)]  *  0 


(5.3,13) 


(5.3,14) 


Therefore,  there  exist  non-zero  vectors  x  and  g,  elements  of  the  state 

o 

space  and  input  space  respectively  such  that 


P(z) 


0 


(5.3,15) 


The  vectors  x^  and  g  are  called  the  invariant  zero  direction  vectors  for  z, 
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or  zero  directions,  for  shore  [17].  Suppose  that  the  system  represented 
by  P(s)  is  excited  by  an  input  of  the  form 

u(t)  =  g  exp(zt) ,  (5.3,16) 

that  z  is  not  an  input  decoupling  zero,  and  z  has  order  1.  If  the  initial 
condition  of  the  system  is  xq,  then  the  response  in  state  space  will  be 
given  by 


x(t)  *  x^  exp(zt) 


(5.3,17) 


but  the  output,  y(t),  will  be  identically  zero  [17].  (Compare  to  the 
introduction  in  Section  4.2). 

Now  suppose  that  the  matrix  P(s)  has  p  distinct  invariant  zeros, 

z^,...,z  ,  of  order  1,  for  simplicity.  Then  for  each  z^  there  exists  a 

state  zero  direction  x*-.  These  state  zero  directions  x^  have  a  number  of 

o  o 

interesting  properties.  For  instance,  the  state  zero  directions  are  un¬ 
affected  by  state  feedback  of  the  form  u=Kx  +  v  [3],  and  the  subspace 

it 

£  * 

spanned  by  the  s tate  zero  directions  is  isomonphic  to  /ft  as  can  be  seen 
from  the  geometric  interpretation  of  Theorem  14  above.  For  a  more  complete 
discussion  of  zero  directions  see  [3],  [17]  and  [30]. 


5.4.  Properties  of  Zeros 

The  following  properties  of  zeros  are  important  in  the  construction 
of  algorithms  for  calculating  zeros.  The  first  theorem  characterizes  the 
basic  properties  of  system  zeros  (Definition  10). 


Theorem  15  [24] 


a)  System  zeros  are  unaffected  by 

1)  Input  space  transformations 

2)  Output  space  transformation 

3)  State  space  transformations 

4)  Linear  output  feedback,  u=Fy+v 

b)  The  system  zeros  of  a  system  (A,B,C,D)  are  the  same  as  those 

T  T  T  T 

of  the  dual  system  (A  ,C  ,B  .D'1'). 

Remark:  Since  transmission  zeros  (Definition  I),  invariant  zeros  (Definition 

5),  and  decoupling  zeros  (Definitions  7-9)  are  all  subsets  of  system  zeros 
(Theorem  9),  it  follows  that  this  theorem  applies  to  these  subsets,  too. 

Proof  of  Theorem  15:  Consider  a  system  matrix  P(s).  Let  T,  G,  V  be  non¬ 
singular  transformation  matrices  such  that 

x  =  Tx1 

u  =  Gu'  (5.4, 1) 

y 1  ■  Vy 

where  the  equations  in  (5.4,1)  define  the  primed  variables.  Applying  these 
transformations  to  the  original  system,  a  new  system  is  obtained,  represented 
by  the  system  matrix  P(s)  where 
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(5.4,2) 


as  can  be  varified  by  direct  calculation.  This  shows  that  P(s)  can  be 
obtained  from  P(s)  by  elementary  row  and  column  operations.  Hence,  P(s) 
has  the  same  Smith  form  as  P(s),  and,  therefore,  the  same  invariant  zeros. 
If  P(s)  has  no  decoupling  zeros,  the  invariant  zeros  are  just  the  trans¬ 
mission  zeros  (Theorem  6)  which  shows  the  transmission  zeros  are  unchanged 
by  the  transformations  in  (5.4,1).  Similar  arguments  show  that 


(5.4,3) 


and 

[sI-A  |  -B]  *  [sI-l"LAT  |  -TBG]  (5.4,4) 

So  the  decoupling  zeros  are  also  unchanged.  Combining  these  results  with 
Theorem  9,  1-3  of  part  (a)  of  the  theorem  follow. 

Now  consider  linear  output  feedback  of  the  form 

u  *»  Fy  +  v  (5.4,5) 

where  F  is  a  mxr  matrix.  The  closed  loop  system  can  be  represented  by 
the  system  matrix  P1 (s)  where 


By  using  elementary  row  and  column  operations,  the  following  reduction  of 
P'(s)  is  obtained: 


Therefore,  the  invariant  polynomials  of  P' (s)  are  the  same  as  the  invariant 
polynomials  of  P(s).  Hence,  the  invariant  zeros  (Definition  5)  of  the 


system  after  the  application  of  feedback  (5.4,5)  are  the  same  as  Che  in¬ 
variant  zeros  of  the  original  system.  As  remarked  above,  it  follows  that 
the  transmission  zeros  are  also  unchanged.  Similarly,  it  follows  that 
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I 

I 


r 


si -A 

-B 

0 

0 

si -A 

-B 

0 

0 

C 

D 

*r 

0 

0 

0 

Zr 

0 

0 

Lm 

-F 

I 

m 

0 

0 

0 

I 

m 

and 


(5.4,9) 


"  — “ 

—  — 

si -A  -B  0 

si -A  0  0 

C  D  -Ir 

o  o  ir 

0  I  -F 

m 

0  10 
m 

C  D  0 

COO 

(5.4,10) 


This  shows  that  the  input  decoupling  zeros  (Definition  7)  and  the  output 
decoupling  zeros  (Definition  8),  respectively,  are  unchanged.  Again,  by 
Theorem  9,  (4)  of  part  (a)  of  the  theorem  follows. 

These  arguments  can  easily  be  used  to  show  that  the  system  zeros 
of  P(s)  are  also  the  system  zeros  of  the  dual  system  PT(s)  where 


T 

P  (s) 


si -A 

-bt 


(5.4,11) 


The  difference  is  that  the  input  (output)  decoupling  zeros  of  (A,B,C,D) 

T  TXT 

become  output  (input)  decoupling  zeros  of  (A  ,C  ,B  ,D  ). 

Now  consider  state  feedback  of  the  form 


u  =  Kx  +  v 


(5.4,12) 


where  K  is  a  mxn  matrix.  Since  the  observability  properties  of  a  system 
are  not  invariant  under  state  feedback,  in  general  the  system  zeros  will 
change  with  the  application  of  (5.4,12).  However,  the  invariant  zeros 


^  '  ... -  .. 


. . 
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remain  the  same. 

Theorem  16 

The  invariant  zeros  are  unaffected  by  state  feedback  of  the  form 
in  (5.4,12). 

Proof:  This  theorem  follows  directly  from  Lemma  2  and  Theorem  12,  which 

ic 

say  that  Z  is  unaffected  by  state  feedback,  and  the  geometrical  inter¬ 
pretation  of  invariant  zeros  from  Theorem  14.  Q 

For  an  alternative  proof,  see  [3] . 

The  application  of  feedback  of  the  form  (5.4,12)  to  a  system 
(A, B ,C , D)  generates  a  new  system  (A+BK,B,C+DK,D) .  It  is  possible  to  call 
the  original  system  the  open  loop  system  and  the  new  system,  the  closed 
loop  system.  Then  Theorem  16  says  that  the  open  loop  invariant  zeros  and 
the  closed  loop  invariant  zeros  are  the  same,  which  is  analogous  to  the  well 
known  classical  result.  A  similar  statement  could  be  made  about  output 
feedback  and  system  zeros. 

Another  property  of  zeros  of  a  scalar  transfer  function  is  that 

under  high  gain  feedback,  the  poles  of  the  system  tend  toward  the  zeros, 

or  toward  infinity  as  the  gain  tends  toward  infinity.  There  is  also  a 
generalization  of  this  result  for  multivariable  systems  [31] [32].  Consider 
a  system 

A  *  Ax  +  Bu  (5.4, 13a) 

y  -  Cx  +  Du  (5.4,13b) 

where  the  system  is  non-degenerate;  i.e.  the  system  matrix  has  full  rank. 
Form  a  closed  loop  system  by  using  output  feedback  of  the  form 

ML  


where  P  is  a  scalar  and  K  is  an  arbitrary  matrix  with  rank  min(r,m) 


Substitute  (5.4,14)  into  (5.4,13b)  and  solve  for  u 


Substitute  (5.4,15)  into  (5.4,13a)  to  obtain  the  closed  loop  system 


In  [32],  the  following  theorem  is  proved 


Given  a  system  (5.4,13),  form  a  closed  loop  system  using  (5.4,14) 


Then  if  m  = r,  the  finite  eigenvalues  of 


coincide  with  the  invariant  zeros  of  (5.4,13).  If  m/r,  then  for 


almost  all"  matrices  K,  the  invariant  zeros  of  (5.4,13)  will  be  contained 


in  the  finite  eigenvalues  of  (5.4,17)  as  P— 35 


Of  course,  this  theorem  holds  for  discrete  systems,  too 


Consider  systems  of  the  form 


86 


where  u€Rm,  x  €  Rn,  y  €  Rr  and  che  matrices  (A,B,C,D)  are  defined  appro¬ 
priately.  This  system  can  be  solved  in  a  straight  forward  way  to  give 


yi  3  (j£o  CAJ"lBut-j)+CAljec 


(5.5,?) 


Introduce  the  following  notation: 


[0,i]  ;  u. 


y[o,i] 


y0 


(5.5,3) 


Form  the  (i+l)r  x  (i+l)m  matrix  Q  and  the  (i+l)rxn  macrix  Q 

[°.i]  ^  _  1*0.11 

ICAO'S  CAl'2S  ...  CB  D  j  IcAl 


Q 

[0,i] 


CAL'2B  CA1-3B  ...  D  0 


I  CB 


[x»i] 


(5.5,4) 


Finally  define  the  (i+1) r  x(n+(i+l)m)  matrix  by 


Qi  =  ^Q[xQ,i] I  Q[0,i]^ 


(5.5,5) 


Now  for  any  input  u^q  ^  and  any  initial  condition  x^  the  output  sequence 


[O.i] 


is  given  by 


Note  that  (5.5,2)  is  just  the  first  row  of  blocks  in  (5.5,6) 


The  matrix  Q.  can  be  thought  to  represent  a  input-output  map  9 


parameterized  by  x,.,  which  sends  input  sequences  u 


infinite  sequences  of  points  in  R  .  Define  the  input  space  U  and  the  output 


A  left  inverse  is  any  operator  9.:  Y-U  such  that 


for  all  input-output  pairs  (u 


exists  when  it  is  possible  to  uniquely  determine  the  input  sequence  which 


produced  any  known  output  sequence.  A  right  inverse  is  any  operator 


Y-U  such  that 


A  more  specific  definition  of  the  type  of  inverses  to  be  discussed  for 
(5.5,1)  will  be  given  below. 
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for  all  input-output  pairs  (u^q  i ]  * ^ [ 0  i]  ^  ant*  ^or  A  r^8ht 

inverse  exists  when  it  is  possible  to  find  some  input  sequence  u^q  ^ 
for  any  given  output  sequence  y^Q  . 

Example  15 

Suppose  that  (5.5,1)  has  zero  initial  conditions.  Then  its 
input-output  map  is  just  its  transfer  function  matrix  G(s)  (2.2,4). 

Now,  it  is  clear  chat  (5.5,1)  is  left  invertible  if  G(s)  has  normal  rank 
m.  Similarly,  (5.5,1)  has  a  right  inverse  if  its  associated  transfer 
function  matrix  has  normal  rank  r„  u 

As  these  preliminary  remarks  show,  the  problem  of  left  and 
right  inverses  is  approached  from  the  input-output  mapping  point  of  view. 
Therefore,  the  analysis  below  will  consider  only  systems  (5.5,1)  which 
are  controllable  and  observable.  It  will  also  be  assumed  that 


—  — 

B 

C1 

*  n 

T 

[°J 

D 

0  ;  i.e . 


(5.5,10) 


that  there  are  no  redunant  inputs  and  outputs.  This  will  be  assumed  through¬ 
out  this  section. 

Left  Inverses 

The  discussion  of  left  inverses  will  consider  two  cases;  when 
Xq  is  known  and  when  x^  is  unknown.  Suppose  that  x^  is  known  but  not 
equal  to  zero.  Then,  from  (5.5,6), 


[0,i]  7[0,i] 


-  Q 


tvijV 


(5.5,11) 
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where  y  is  che  known  output  sequence,  is  the  output  sequence  of 

(5.5,1)  for  the  same  input  sequence  but  with  zero  initial  conditions. 
Hence,  if  Xq  is  given,  without  loss  of  generality,  it  can  be  taken  to 
be  zero.  This  is  done  below. 

Introduce  the  following  definition: 

Definition  15  [ 33] 

A  system  (5.5,1)  has  an  L-delay  left  inverse  if  for  any  non¬ 
negative  integer  i,  the  input  segment  u  can  be  uniquely  determined 

l ,  i  J 


by  the  output  segment  y 


[0, i+L] * 


The  immediate  goal  is  to  establish  the  conditions  under  which 
an  L-delay  left  inverse  exists  for  a  system  (5.5,1). 

Lemma  5  [ 33] 

A  system  (5.5,1)  has  an  L-delay  left  inverse  if  and  only  if  for 


an  output  segjnent  y 


[0,L]‘ 


(a)  if  Xq  is  given,  then  Uq  can  be  uniquely  determined 

(b)  if  Xq  is  not  known,  then  Xq  can  be  uniquely  determined 
Proof:  First,  assume  Xq  =  0.  Necessity  follows  directly  from  Definition  15. 


Conversely,  suppose  an  output  segment  y  ^  uniquely  determines  Uq.  Then 
the  effect  of  Uq  on  the  entire  output  sequence  can  be  subtracted  out: 


„  i- 1 
CA  B 


[0,i]  "  1 0, i] 


ca;l‘2b  cal'3.  .ce  d 


u  =  CB 
0  D 


u[l,i]  (5.5,12) 
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where  (5.5, L2)  is  obtained  from  (5.5,6)  with  xn  *  0.  Observe  that 


y[l,i]  *Q[0,i-l]U[l,i] 


(5.5,13) 


3y  hypothesis,  if  i-lal,  (5.5,13)  can  be  solved  to  uniquely  determine  u^. 
Now  induction  establishes  part  (a). 

Now  assume  that  Xq  is  unknown  and  that  (5.5,1)  has  an  L-delay 
left  inverse.  Then  u^q  ^  can  be  uniquely  determined  from  the  output 
segment  y^  by  Definition  15.  In  particular,  u^q  ^  can  be  found. 

Then  (5.5,6)  yields 


y[0,n]  '  Q[0,n]U[0,n]  “Q[xQ,n]X0 


(5.5,14) 


But  Qf  .  is  the  observability  matrix.  (See  (5.5,4)).  Since  (5.5,1)  is 

l^Q.nj 

assumed  to  be  observable,  P[Qr  ,]  = n.  This  implies  (5.5,14)  can  be 

LXQ>nJ 

solved  uniquely  for  x^. 

Suppose  Xq  is  uniquely  determined  by  y^  .  Then  the  effect 

of  xn  on  yf  .  can  be  subtracted  off  as  noted  in  (5.5,11).  This  gives 
U  l  0,  LJ 

Che  equacion 


yo  "  Cxo  3  Duo 


(5.5,15) 


From  (5.5,11)  the  following  relationship  is  obtained  by  dropping  the  blocks 
associated  with  the  first  time  unit: 


y[l,i]  “Q[0,i-l]U[l,i]  +Q[BuQ,i-l]Bu0 


(5.5,16) 


so  that  y  .  is  the  output  sequence  to  (5.5,1)  on  [1,»)  with  an  input 

l  *  9  *  4 

sequence  Up  and  initial  condition  x^  * Buq*  By  hypothesis,  x^  can  be 
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determined  from  the  segment  y^  l+j.]  '  assumption  (5.5,10),  it  follows 

that  x^=Buq  an<i  (5.5,15)  taken  together  uniquely  determine  u^.  An 
induction  argument  completes  the  proof.  Q 

The  following  theorem  gives  conditions  for  che  existence  of  an 
L-delay  left  inverse  in  terms  of  Algorithm  I  (Section  5.2).  Recall  that 

is  the  i-th  unknown-input  unobservable  subspace,  Definition  11.  Further¬ 
more,  throughout  the  rest  of  this  section,  shall  represent  the  matrix 
F^  defined  in  Algorithm  1,  (5.2,15). 

Theorem  18 


a)  If  Xq  is  unknown,  an  L-delay  left  inverse  for  a  system 
(5.5,1)  exists  if  and  only  if  £L+^ *  0. 

b)  If  Xq  is  known,  an  L-delay  left  inverse  for  a  system  (5.5,1) 
exists  if  and  only  if  P[FL+jJ  =®« 


Proof:  By  equation  (5.5,6) 


y[o,i] 


(5.5,17) 


Compare  (5.5,17)  to  equation  (5.2,17).  They  are  exactly  the  same  if  the 
row  of  blocks 

[MqA\  M(}Ai'1B,...,M0AB,  MQB]  (5.5,13) 

is  attached  to  the  top  of  and  MQxi+j.  atcached  to  the  top  of  the 
vector  y j q  ^ .  If  the  proof  of  Theorem  10,  where  (5.2,17)  appears,  (5.2,17) 
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S  such  that 


5Q,  -  Qt 


under 

the  assumption 

X 

o 

H 

0. 

It  follows  chat  Q. 

n 

(5.2,2 

0);  i.e.  there 

exists  a 

non-singular  matrix 

G1Ai‘:B  ... 

G1B 

F1 

c,*1'2 

G0A  B  •  •  • 

F, 

0 

** 

• 

* 

"  ^[x0,ij^[0,i]] 

Gi+1 

“i+1 

Fi+i 

0  .... 

•  • 

0 

(5.5,19) 

‘ 

where  the  matrices  ,  F^ ,  are  produced  by  Algorithm  1. 

To  prove  part  (a),  assume  that  Xq  is  unknown.  Firsc,  suppose 
(5.5,1)  has  an  L-delay  left  inverse.  If  Xq  can  be  uniquely  determined  for 
all  Xq,  it  follows  from  (5.5,19)  that  ?i[M,  +  ^]  “0,  By  Theorem  10, 

Jll+i  *  Conversely,  suppose  £L+1  » 0.  This  implies  M^]  -  0. 

Then  (5.5,19)  shows  that  x^  can  be  uniquely  determined  from  y^  .  Hence, 
by  Lemma  5,  an  L-delay  left  inverse  for  system  (5.5,1)  exists. 

To  show  part  (b) ,  assume  that  x^  »  0.  Then  equations  (5.5,17) 
and  (5.5,19)  reduce  to 


Sy[0,L]  *Q[0,L]U  [0,L] 


Suppose  chat  F^^  has  rank  m,  i.e.  is  a 
Then  u^  ia  uniquely  determined.  By  Lemma  5, 
left  inverse. 

Now  suppose  that  system  (5,5,1)  has 
Fl+^  doesn't  have  full  rank.  Note  chat 


(5.5,20) 

square,  non-singular  matrix, 
the  system  has  a  L-delay 

an  L-delay  left  inverse  but 


KBESBOiaMI 


(5.5,21) 


I 

I 

I 

1 

1 

I 


L 

•  m 
[ 


L 

I 

I 

1 

I 

I 

I 
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XlrL)  =>  fl[F  1 


(This  is  easily  seen  from  Che  structure  of  columns  that  the  F^'s  are 
associated  with) .  Hence,  if  thun  *I[Fj  *0,  i  ■  1,2  , . . .  ,L. 

It  follows  chat  Uq  can't  always  be  uniquely  determined,  contradicting 
Lemma  5 .  O 


Part  (b)  of  Theorem  18  was  first  proved  by  Silverman  and  Payne 

(34], 

A  number  of  known  results  follow  immediately  from  Theorem  18. 


Corollary  2 

The  following  statements  are  true  for  a  system  with  unknown 
initial  conditions: 


(a)  There  exists  a  least  integer  q^  such  chat  if  the  system 
(5.5,1)  has  a  q-delay  left  inverse  and  a  q-delay  left  inverse  for  every 
qa  V 


(b)  q  s n  if  q  exists. 

v  o  o 

(c)  The  following  statements  are  equivalent  for  Che  system 

(5.5,1): 

(1)  (5.5,1)  has  a  q-delay  left  inverse. 


(2)  £*  -  0 

(3)  There  is 
followed  by  zero  inputs 


no  x,  and 
0 

such  chat 


no  non- zero  input  sequence 
the  output  sequence  ^ 


U[0.n] 

is  identically 


zero  for  all  t. 


Proof:  To  see  part  (a),  recall  chat  is  monotonically  nonincreasing 

(Lemma  1,  Section  5.2).  If  d(£^)  *  0  for  some  q,  then  there  exists  a 
first  integer,  q^ ,  such  that  d(£^  )  *0  and  d(£^)  *0  for  ail  qiq^.  The 
result  follows  by  Theorem  18. 

Part  (b)  is  also  a  consequence  of  Lemma  1  and  Theorem  18.  Since 
£n  *  ,  j  »1,2,...  (by  Lemma  1),  if  £^0,  then  the  system  (5.5,1) 

doesn't  have  a  left  inverse  (Theorem  18).  Hence,  if  a  q-delay  left  inverse 
exists  tor  a  system  (5.5,1),  then  q^  ^ n  by  part  (a). 

In  part  (c),  the  equivalence  of  (1)  and  (2)  is  a  consequence  of 

it 

parts  (a)  and  (b)  and  Theorem  12  wnich  states  chat  £  Ftom  Defini¬ 

tion  11,  it  is  seen  that  (3)  simply  states  that  » 0.  By  Theorem  12, 
this  is  equivalent  to  £  *  0.  □ 

The  corollary  characterizes  to  what  extent  past  values  of  x^ 
affect  the  present  value  of  the  output.  The  existence  condition  in  part 
(o'*  was  first  established  by  Bengtsson  [2]  for  the  case  D*0. 


Corollary  3 


The  follows  statement  are  true  for  a  system  (5.5,1)  with  known 


initial  conditions. 


(a)  If  a  system  has  an  L-delay  left  inverse,  then  there  exists 
a  least  integer  Lq  such  that  the  sytem  has  a  L^-delay  left  inverse  and  an 
L-delay  left  inverse  for  every  LaL^. 

(b)  LQ  s;  n  if  LQ  exists  . 

(c)  The  following  statements  are  equivalent. 
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(1)  A  system  has  a  L-delay  left  Inverse. 

(2)  ptQ^0  Ljl  -  P tQ ^ 0  L j  1  *m.  Furthermore,  LQ  is  the 

first  integer  for  which  this  is  true. 


(3)  There  exists  no  input  segment  u^  ^  ,  not  equal  to  zero, 
followed  by  all  zeros  which  produces  an  identically  zero  output  sequence. 


(4)  *1  -  0. 


Proof:  The  properties  of  this  corollary  follow  from  the  properties  of  the 

matrix  and  its  reduced  form  ^  defined  in  (5.5,19). 

Part  <,a)  follows  from  the  fact  that  the  first  column  of  blocks  of 

Mr.  .,  has  a  monotonies  1 ly  non-decreasing  column  rank.  If  P[F  ]  »m 
l  U  *  *-•  J  Lq+I 

(the  column  rank  of  the  firsc  column  of  blocks),  then  ^[F  ]  »m  and  LiL.+l. 

La  U 

Theorem  18  completes  the  proof  of  part  i,a). 

The  Cay  ley  -  Hamilton  theorem  applied  to  the  first  column  of  blocks 
proves  part  (b). 

To  prove  part  t,c)  consider  che  matrix  Qjj  .  This  matrix  shows 

.  e . 


the  rank  of  both  and  l-1]’ 

*'V,L]1  •  i=l  P['l> 


LfL 

y 


(5.5.22) 


#"WU1  "ill 


so  that 


pfQ[0,L)^  'pfQ[0,L-lj]  *,pfFL+l) 


l5.5,:3) 


Now  Theorem  18  establishes  the  equivalence  of  (1)  and  ^C).  It  is  also 


Ho 


clear  chat  the  last  remark  follows  from  part  ia3. 

To  show  the  equivalent  of  (l-)  and  (33,  suppose  the  system  is 
left  invertible.  Then,  by  Theorem  13  and  part  (b\  F^  is  a  square  non¬ 
singular  matrix  for  some  l1*!'.  Hence,  a  non- re ro  u  will  produce  a  non- 

o 

zero  output.  Conversely,  suppose  (33  but  the  system  is  not  invertible. 

Then  it  follows  chat  p(F^]<m,  i»0,l,...,n.  Then,  it  is  possible  to  find 
a  non-zero  input  segment  nj  followed  by  all  zeros  which  would  produce 
an  identically  zero  output.  This  contradiction  establishes  the  equivalence 
of  (l)  and  (3). 

*  s 

That  (3)  implies  (,43  follows  from  the  fact  chat  >i  >-  n  (See 

equation  (5.2,303).  From  Definition  13,  it  follows  that  ^33  implies 
*  * 

il  "0.  Hence,  il  “0.  How  assume  that  >\  »0  but  chat  the  svscem  has  no  L- 
n 

delay  left  inverse.  3y  (33,  there  exists  a  non-zero  input  segment  u^  nj 
followed  by  all  zeros  which  produces  an  identically  zero  output.  Clearly, 

u  €~[D]  i  ■  0, 1 , . . . , n  (5.5,24) 


It  follows  that 


i  *  0 ,  l . n  (5.5,25) 

since 

“tB]nntD]-0  i5.5,2b3 


by  assumption.  Therefore,  this  input  segment  generates  a  state  sequence 

x^,  t*l . n+1.  By  Definition  11,  these  x^  are  elements  of  A^.  For  if 

x.,  *x.,  then  the  control  sequence 

0t’ 


ui+j  0sJSn-i 


0  n-i  <  J  £  n 


wtlL  produce  an  identically  zero  output.  Since 


(5.5,27) 


-  Buq  +  0 


(5.5,28) 


it  follows  that  £^^0.  But,  by  construction,  x^  is  also  an  element  of 
Definition  13.  The  input  segment 


0  OSjin-2 


u0  J=*n-1 


shows  this.  Now,  by  Theorem  13, 


(5.5,29) 


a  -a  ft  z  t  $ 

n  n 


(5.5,30) 


This  contradiction  shows  that  (4)  implies  (1).  O 

Comparison  of  Corollary  3  with  Corollary  4  shows  that  there  are 
similar  results  for  systems  with  known  and  unknown  initial  conditions. 

The  main  difference  is  that  if  the  initial  condition  is  unknown,  it  is 
necessary  to  reconstruct  the  state,  where  as  if  the  initial  condition  is 
known  this  is  not  necessary.  This  leads  to  the  tighter  restrictions  on 
systems  with  unknown  initial  conditions.  (£  =0  implies  ft  = 0,  but  not 

conversely) . 

The  integer  Lq  defined  in  part  (a)  of  Corollary  3  is  called  the 
inherent  integration  of  the  system  [35].  It  should  be  noted  that  this 


number,  Lq»  Is  not  the  same  as  the  integer  that  was  introduced  in  part  (a) 


of  Corollary  2.  Silverman  and  Payne  have  shown  that  LQ s qQ  [34]  but 
this  is  all  that  can  be  claimed  for  left  invertible  systems.  Sain  and 
Massey  have  shown  [35]  that  if  G(s)  is  a  transfer  function  matrix  for  a 

A 

left  invertible  system  with  inherent  integration  Lq  and  G(s)  is  the  transfer 
function  matrix  for  that  left  inverse,  then 

rm 

G(s)G(s)  *  (5.5,31) 

s  o 

The  number  Lq  characterizes  the  extent  to  which  past  values  of  the  control 
affect  the  present  value  of  the  output. 

Part  (c)  of  Corollary3  gives  3  equivalent  characterizations  of 
systems  with  left  inverses.  Sain  and  Massey  [35]  proved  (2)  and  (3); 

Wonham  and  Morse  [36]  proved  (4)  for  the  case  D=Q.  Corollary  3  and  Theorem 
18  together  show  the  equivalence  of  the  characterizations  of  left  inverse 
systems  given  by  Sain  and  Massey  [35]  and  Silverman  and  Payne  [34].  This 
equivalence  was  also  noted  in  [37]  and  [38]. 

Right  Inverses 

The  right  inversion  problem  for  (5.5,1)  can  be  analyzed  in  much 
the  same  way  as  left  inverses  were  analyzed  above.  The  right  inversion 
problem  occurs  when  it  is  desired  to  find  an  input  sequence  for  (5.5,1) 
which  will  produce  some  specified  output  sequence.  Note  that  this  problem 
is  slightly  different  in  that  the  specified  output  sequence  is  not  known 
to  be  in  the  image  of  the  operator  3  (5.5,7).  This  leads  to  some  differences 
in  the  way  the  initial  condition,  Xq,  is  handled.  Three  cases  are  distin¬ 
guished:  1)  Xq  is  unknown,  2)  Xq  is  known  and  fixed,  and  3)  SCq  is  known, 

but  can  be  chosen  to  match  the  specified  output  sequence.  To  motivate  the 


rest  of  che  discussion  of  initial  conditions,  consider  the  following 
definition: 

Definition  16 

A  system  (5.5,1)  has  an  L-delay  right  inverse  if  for  every 
non-negative  integer  i,  an  input  segment  u^q  ,  which  produces  the  given 
output  segment  y^  ^  ,  can  be  determined  from  the  output  segment 

y[0,i+L]  13 

Definition  16  differs  from  Definition  15  in  rwo  ways.  First, 

the  input  sequence  in  Definition  16  is  not  required  to  be  unique.  Secondly, 

Definition  16  implies  that  the  output  sequence  y^  ^  is  in  the  unage  of 

the  operator  3.  Since  3  is  parameterized  by  x^,  the  class  of  output 


sequences  y 


[0,i] 


for  which  it  is  possible  to  obtain  an  input  sequence 


u^q  is  dependent  on  x^.  To  see  this  explicitly,  recall  the  matrices 
(5.5,4)  and  (5.5,19),  which  relate  the  input  and  output  sequences 
as  follows: 


y[0,i]  =  3y[o,i] 


3  Q[x0,i]X0+Q[0,i]U[0,i] 

(5.5,32) 


where  S  is  defined  in  (5.5,19).  Consider  again  the  three  cases  listed 
for  the  initial  condition  Xq.  First,  suppose  that  the  initial  condition 
is  unknown.  Then  Xq  must  be  determined  if  (5.5,32)  is  to  be  solved  for 
u^q  .  This  implies  that  *0  where  is  defined  in  (5.5,19). 

This  is  quite  a  strong  condition  but  further  treatment  of  the  case  is 


deferred  until  later.  Note  that  if  x~  can  be  determined,  then  the  calcula- 
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cion  of  a  suitable  u^g  depends  on  Che  Mcrix  Q^g  ^ . 

Secondly,  suppose  chac  Xg  is  known  and  can  be  chosen  arbitrarily. 
Then  from  Che  boctom  row  of  blocks  in  (5.5,02)  (see  (5.5,19)) 


y[0.i]  ■  Wo 


(5.5,33) 


muse  be  solvable  for  some  Xg  where  y^g  ^  is  the  appropriate  subveccor 
of  y^Q  .  This  unplies  chac  y^g  ^  €>l[M.  +  i]  wnich  puts  a  conscrainc  on 
che  class  of  sequences  y ^ g  for  which  (5.5,32)  is  solvable.  Again,  if 

(5.5,33)  is  satisfied  Chen  Che  question  of  che  existence  of  a  right  inverse 
rests  on  che  matrix  Q^g  . 

Finally,  suppose  xQ  is  known  and  fixed.  Then  (5.5,33)  uniquely 
specifies  y^  ^  and  all  output  sequences  y^  ^  muse  satisfy  (5.5,33)  Co 
be  in  che  image  of  If  this  requirement  is  met,  che  existence  of  a 
suitable  input  sequence  u^-q  ^ j  depends  on  Q^.g  .  j  as  above. 

Two  remarks  are  in  order.  First,  the  condition  in  (5.5,33)  can 
be  decomposed  into  constraints  on  certain  subvectors  of  each  output  vector 
yL  [34]  but  che  details  will  not  be  given  here.  Secondly,  suppose  chat 
r  £ m  in  (5.5,1)  and  that  3[D]  = r.  Then  che  application  of  Algorithm  I  to 
chis  system  immediately  yields  3 0  for  all  i.  In  this  case,  the  constrait 
(5.5,33)  doesn't  exist  and  an  appropriate  input  sequence  u^g  ^  can  be  found 
for  every  output  sequence  y^g  ^ j  (as  will  be  shown  below). 

Now  suppose  chac  che  proposed  output  sequence  meets  the  require¬ 
ments  imposed  by  the  initial  condition,  and  that  the  effect  of  the  initial 
condition  has  been  subtracted  out,  as  in  (5.5,11).  The  possible  existence 
of  an  appropriate  input  sequence  u^g  is  still  to  be  determined.  This 
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question  is  answered  by  the  following  lemma  and  theorem. 


Lemma  6 


A  system  (5.5,1)  has  an  L-delay  right  inverse  if  and  only  if 


Uq  can  be  determined  from  y^  • 

Proof:  Exactly  the  same  as  Lemma  5,  part  (a). 


The  following  theorem  characterizes  L-delay  right  inverses. 


Theorem  19 


(a)  A  system  (5.5,1)  has  an  L-delay  right  inverse  if  and  only 


if  p[FT4.il  Br- 


(b)  If  a  system  has  an  L-delay  right  inverse,  then  there  exists 
a  least  integer  LQ  such  that  the  system  has  an  LQ-del ay  right  inverse  and 
an  L-delay  right  inverse  for  every  L*Lq. 

(c)  If  Lq  exists  then  LQsn. 

(d)  The  system  (5.5,1)  has  an  L-delay  right  inverse  if  and  only 
if  ptQ[0,L]~  PIQ[0  L-l]^  = r‘  Furthermore,  Lq  is  the  least  integer  for 
which  this  is  true. 

Proof:  This  theorem  follows  from  the  equation 

y[0,L]  =S(y[0,Lf  Q[xq,L]X0)  =SQ[0,L]U[0,L]  -VLJVL] 

(5.5,34) 

where  S  and  Lj  are  defined  in  (5.5,19).  Note  that  the  matrices 
appear  in  Q[0jL]. 

To  show  part  (a),  suppose  that  P[FL+1]  * r  for  some  L.  Then 
(5.5,34)  shows  that  it  is  possible  to  determine  a  un.  By  Lemma  6,  the 


■ 


system  has  a  L-delay  right  inverse.  Now  suppose  the  system  has  an  L-delay 
right  inverse  but  P[F^+iJ  ^r>  Then  from  (5.5,34),  it  is  easy  to  see  that 
Uq  can't  be  determined  for  every  permissible  output  sequence.  Hence  by 
Lemma  6,  a  L-delay  right  inverse  doesn't  exist.  This  contradiction 
establishes  part  (a). 

Part  (b\  follows  by  noting  that  the  row  rank  of  the  top  row  of 
blocks  in  ^on-decreasing.  The  arguement  is  similar  to  the  one 

used  to  prove  part  ^a)  of  Corollary  3. 

Part  tc"*  is  established  by  applying  the  Cay ley-Hamilton  theorem 
to  the  top  row  of  blocks  of  and  using  part  (b). 

Part  (d)  follows  directly  from  part  (a).  The  details  arc  similar 
to  the  proof  of  (2)  of  part  (c)  in  Corollary  3.  O 

Hie  criteria  for  the  existence  of  a  right  inverse  in  part  taN 
was  first  presented  by  Silverman  and  Payne  [34];  the  criteria  in  parts 
(b)-(d)  by  Sain  and  Massey  [35].  Definition  lb  for  a  system  with  a  known 
initial  condition  was  first  Introduced  in  [39]  and  called  funct ional 
reproducibility. 

Consider  again  the  question  of  a  right  inverse  for  a  system 
(5.5,1)  with  unknown  initial  conditions.  Suppose  that  such  a  system  does 
have  a  L-delay  right  inverse.  By  the  remarks  proceeding  Lemma  b,  it  follows 
Chat  ?([M1+1]  -0.  By  Theorem  10,  this  implies  “  0*  Now  by  Theorem  IS, 
this  system  also  has  a  L-delay  left  inverse.  A  standard  result  says  that 
in  this  case  the  inverse  is  unique  and  the  distinction  between  right  and 


left  need  not  be  made. 

There  is  a  duality  between  left  and  right  inverses.  Consider 
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the  dual  system  of  (5.5,1) 


x 


i+l 


y 


i 


T  T 

»  A  x^  +  C 

T  T 

•  B  x^  +  D 


(5.5,35) 


If  Xq-0,  the  solution  to  this  system  is  given  by 


y[0,i]  "  Q  [0,L]U[0,L] 


(5.5,36) 


where  Qrrt  ,  is  defined  in  (5.5,4).  Now  the  duality  between  Corollary  3 
10,  LJ 

and  Theorem  19  is  completely  obvious. 

Theorem  20  [ 35] 


If  the  initial  condition  of  (5.5,1)  is  zero,  then  the  system  (5.5,1) 
has  an  L-delay  left  inverse  if  and  only  if  its  dual  system  (5.5,35)  has  an 
L-delay  right  inverse.  □ 

Example  16 


In  Example  15,  it  was  pointed  out  that  a  transfer  function  matrix 

G(s)  has  a  left  inverse  when  its  normal  rank  in  m.  In  this  case,  it  is 

T 

clear  that  the  dual  system  G  (s)  has  a  right  inverse. 

For  completeness,  it  is  noted  that  if  the  system  (5.5,1)  has 
zero  initial  conditions,  then  its  transfer  function  matrix  G(s)  has  normal 
rank  P[Fn+1l.  See  [34]  for  details.  □ 

Minimal  Order  Left  Inverses 

It  has  already  been  shown  that  for  transfer  function  matrices, 
the  transmission  zeros  of  a  transfer  function  matrix  are  the  poles  of  an 
inverse  transfer  function  matrix  (Section  3.4).  In  what  follows,  the  left 
inverse  of  a  system  (5.5,1)  will  be  constructed  to  show  that  this  result 
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generalizes  to  state  space  systems.  It  will  be  assumed  that  the  left 
inverse  system  will  have  the  form 


A 

a’i+l  "A  a'i'f  NI(P)>’i 
«i  "C  UL'i+  N2(p)yi 


(5.5,37) 


where  N^(s)  and  N^(s)  are  polynomial  matrices  and  p  is  the  unit  delav.  A 
minimal  left  inverse  is  any  left  inverse  with  representation  (5.5,37)  such 
that  has  minimal  dimension  [2]. 

Bengtsson  (2]  has  given  a  construction  of  a  minimal  left  inverse 
for  systems  with  D“0  and  unknown  initial  conditions.  To  study  systems 
with  D^O,  it  is  necessary  to  further  delineate  the  structure  of  the  D 
matrix.  By  applying  input  and  output  transformations,  the  system  (5.5,1) 
can  be  reduced  to 


x 


i+1 


-  Axi  +  [Bx  B2) 


(5.5,38a) 


(5.5,38b) 


where  is  a  square  non-singular  matrix.  Let  C,,  be  a  r^xn  matrix.  The 
following  three  lemmas  will  be  needed  in  the  construction  of  a  left  inverse 
for  (5.5,38). 
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Lemma  7 


tIc 

If  £  -  0  for  the  system  (5.5,38),  then  £  -0  for  the  system 


(a,b2,c2,0). 

Proof:  Suppose  that  £*  -  0  for  (5.5,38)  but  £*  f1  0  for  (A,  B^,  ,C2,0) .  By 
Theorem  12,  £  *^n*  If  £n  +  0,  there  exists  an  x^  i  0  and  an  input  segment 
u[0  n-L]  suc^  chat  c^e  output  segment  is  identically  2ero  (Defini¬ 

tion  11).  This  input  segment  generates  as  state  segment  x^,  i*l,...,n. 

Now  consider  the  full  system  (5.5,38)  with  x^.  The  equation 


y0«°  -CjX^  +  DjUq 


(5.5,39) 


can  be  solved  for  u^  since  is  square  and  non-singular.  Furthermore,  the 


equation 


x'  -  Ax'  +  B2u'  -  Ax'  +  B^uq  +  B2uq 


(5.5,40) 


yields 


Vo  '  V*  -  Bl“o 


(5.5,41) 


By  assumption  (5.5,10),  the  special  structure  of  (5.5,38)  shows  that 
/ < [ B 2 ]  -0.  It  follows  that  B0  has  a  left  inverse  B*  so  that  (5.5,41)  can 

A 

be  solved  for  UqJ  i.e. 


"0  *  "o  •  B2Bl“o 


(5.5,42) 


I 
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v 


0 


(5.5,43) 


where  uQ  is  obtained  from  (5.5,39)  and  u^  is  given  by  (5.5,42).  By  induc¬ 
tion,  this  procedure  will  generate  an  input  segment  v^q  n_i]  w^ich  produces 
the  state  segment  x|,  i = 1, . . . ,n,  starting  from  Xq,  such  that  the  output 

A 

is  identically  zero.  (Note  that  y^*C2x^“0  by  choice  of  x^.  Equation 

(5.5,39)  shows  that  y^  is  zero).  Hence,  by  Definition  11,  Xg€£nJ*0,  where 

is  defined  for  the  full  system  (A,B,C,D).  By  Theorem  12,  this  is  a 
* 

contradiction  as  £  =0  by  hypothesis.  □ 

Note  that  this  proof  actually  shows  that  £^  for  the  system 
(A,B9,C0,0)  is  the  same  as  £^  for  the  full  system  (5.5,38). 

Lemma  8 

* 

If  £  =  0  for  the  system  (5.5,38),  then  there  exists  n  x r^  matrices 
N.  such  that 

i 

qo  k 

Z,  nN  .  C„A  -  I 
k«0  q  -k  2  n 

o 

q 

CASE'S  ■  ° 


where  q^  is  defined  in  Corollary  2. 

★ 

Proof:  Since  £  "0  for  the  system  (5.5,38),  £^  *  0  as  was  shown  in  the 

o 

proof  of  Corollary  2.  By  the  remark  following  the  proof  of  Lemma  7, 

£  +1  -  0  for  the  system  (A,B2,C2>0).  Consider  the  matrix  Q  (5.5,4)  de- 

V  _  qo 

fined  for  (A,B2»C2,0).  This  matrix  can  be  reduced  to  Q  as  was  done  in 

qo 

the  proof  of  Theorem  18  (5.5,17-19).  By  Theorem  10,  +^]  *0  where  M 

_  _  q°  V 

occurs  in  Q  .  From  the  structure  of  Q  ,  it  is  seen  that  ?2[M  *  0 

qo  qo  V 
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N  y  «N  C„x. 
<ln  1  2  L 


Nq0-lyi+l  “Nq0-l(C2Axi  +  C2B2Ui) 


1 N  .C.Ax. 

vl  2  i 


(5.5,46) 


by  Lemma  8.  This  gives 


«  y,+N  ,y  .  =  (N  +pN  ,)y  *  (N  C.+N  ,C„A)x. 

^  ^  q°  ^  q°  q°  (5.5,47) 


If  this  procedure  is  continued  for  q^+1  steps,  it  is  found  that 

q0  Vk  »  k 

(Ik-op  Vyi  "N(p)yi  -  Ok-oV-k0*  )xi =  xi 


(5.5,48) 


LI 


again  by  Lemma  8.  Now,  from  the  system  (A,B2>C2>0) 

(pI-A)xt  “B2ui 

Substituting  (5.5,48)  into  (5.5,49)  yields 


(5.5,49) 


LJ 


Ui“B2  (PI-A)N(P>yi 


(5.5,50) 


+  -f 

where  B0  is  a  left  inverse  of  B0.  Note  that  B2  exists  as  7![B9]  =0  from 
the  requirement  7J[B]  0  7i[D]  -0  for  the  system  (5.5,38).  Equations  (5.5,48) 
and  (5.5,50)  are  the  required  results,  as  (5.5,46-50)  hold  for 
all  ilO.  D 

From  Lemma  9,  the  complete  left  inverse  of  (5.5,38)  is  easily 
constructed.  From  (5.5,38b) 


yi"cixi+Vi 


(5.5.51) 


or 


0 

u 

ii 


! 

I 

B 


0 

0 

fi 

I: 

1 . 


0 

D 

II 

0 

B 

!. 

I 

I 

I 
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VDIlVDIlcixi 


(5.5,52) 


since  is  square  and  non-singular.  Combining  equations  (5.5,50)  and 
(5.5,52),  it  follows  that 


-1 

-Dl  C1N(p) 

*i 

■f 

A 

B2(pI-A)N(p) 

*i 

_  _ 

(5.5,53) 


-  N  (p)  yi 


Equation  (5.5,53)  holds  for  all  iiO.  Hence,  (5.5,53)  combined  with  (5.5,48) 
constitute  a  left  inverse  for  (5.5,38).  The  following  theorem  has  been 
proved : 

Theorem  21 

If  a  system  (5.5,38)  with  unknown  initial  condition  x^  has  a 
qQ-delay  left  inverse,  then  there  exists  a  polynomial  matrix  N(s),  defined 
by  Lemma  8,  and  a  second  polynomial  matrix  N(s),  defined  by  (5.5,53)  such 
that 


■N(p)yi  -  [0  N(p)]yi 
u L  "N(p)yt 


(5.5,54) 


where  p  is  the  unit  delay,  holds  for  all  input  and  output  sequences  and  all 
Xq!  i.e.  (5.5,54)  is  a  left  inverse  for  (5.5,38).  □ 

Note  that  this  left  inverse  is  minimal  as  Che  dimension  of 


in  (5.5,37)  is  zero. 
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This  construction  was  given  by  Bengtsson  [2]  for  systems  with 
D  sO.  Lemma  8  also  appeared  there  for  the  special  case  q^n.  The  left 
inverse  system  (5.5,54)  can  be  interpreted  as  a  bank  of  delays,  N(p), 
followed  by  a  dynamical  system.  From  the  construction  of  N(p)  given  in 
the  proof  of  Lemma's  8  and  9,  it  is  seen  that  qg  is  the  minimum  number  of 
delays  needed  to  implement  this  inverse  system.  Hence,  the  construction 
of  the  left  inverse  given  here  is  for  a  larger  class  of  systems  and  it  also 
identifies  the  maximum  order  of  a  polynomial  in  N(p);  i.e.  the  maximum 
number  of  delays  needed  to  implement  the  left  inverse  system. 

Now  consider  a  system  (5.5,38)  with  zero  initial  conditions.  If 

•k 

this  system  has  a  left  inverse,  by  Corollary  3  R  *0.  By  Theorem  11,  there 
exists  a  feedback  matrix  K  such  that  (5.5,38)  can  be  placed  in  the  following 
canonic  form: 

*i+l  *  (A+BK)xi  +Bvi 


y^^  *  (C+DK)xi  +Dv^  * 

u^  *  Kx^  +  v^^  *  KjXt  +  K2xi  +  vi  (5.5,55c) 

In  general,  this  transformation  will  require  a  change  of  basis  in  state 
space.  Without  loss  of  generality,  it  will  be  assumed  that  this  has  been 


An  0 


A21  A22 


x. 

i 


B11  B12 


B21  B22 


v, 

l 


(5.5,55a) 


yi 


C11  ° 


C21  ° 


x. 

i 


D,  0 


0  0 


v. 

l 


(5.5,55b) 


Ill 


done.  In  comparing  (5.5,55)  Co  canonic  form  (5.2,24),  note  Chac  il  "0 
for  (5.5,55)  implies  chac  ■  0  where  B^  is  defined  in  (5.2,24).  This 
implies  chac  B^*[B^  B^]  and  B.,  ■  [B.,^  ]  where  B^  and  B.,  are  defined 

in  (5.2,24).  Of  course,  is  a  square  non-singular  macrix  defined  in 
(5.5,38). 

The  firsc  seep  in  constructing  a  left  inverse  for  the  system 
(5.5,55)  is  defined  in  the  following  lemma: 

Lemma  10 

If  the  system  (5.5,55)  has  a  left  inverse  when  Xq  ■  0,  there  exist 
polynomial  matrices  N(p)  and  N(p)  such  thac 

xt  -  N(p)yt 


7i 

Di1  ■DilcnN(p) 

*1 

V.  * 

as 

i 

a 

4" 

_vi_ 

0  B12(pI-A)N(p) 

yi 

-  -J 

where  B^  is  a  left  inverse  of  B^ ,  and  p  is  a  unit  delay. 

Proof:  The  input-output  operator  of  (5.5,55)  is  found  by  neglecting  the 
unobservable  part  of  (5.5,55);  i.e.  it  corresponds  to  the  input-output 

operator  of  the  subsystem  (AU,[BU  B^J,  Cu,  [Dj_  0]).  Since  £*  is 

★ 

maximal  for  (5.5,55),  it  follows  that  £  for  (Au>  [B11  B12l,  cu> 

[D  0])  is  zero.  Now,  if  it  can  be  shown  thac  a  left  inverse  exists  for 
Chen  Theorem  21  can  be  applied  to  (Aj^,  [B^  B^l  ,  C^,  (D^  0]) 

Co  give  the  desired  result.  To  show  the  existence  of  a  left  inverse  for 
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Bj^»  it  will  be  shown  chat  =  0.  Suppose  that  7J[B^2]  J4  0.  Since 

2?[B]  0  ?2[D]  *  0  by  assumption,  it  follows  that  ??[B^2]  f\  7?[B?2]  “0.  Chose 

a  non-zero  input  sequence  u^q  such  that  u^  €/2[B^2)  for  all  j.  This 

input  sequence  generates  a  non-zero  state  sequence  x^^  since  u^  £7J[ B22]  ^or 

all  j.  The  structure  of  (5.5,55)  shows  that  x  ^  =  0  for  all  j  and  that  the 

★ 

output  is  identically  zero.  By  Definition  13,  it  follows  that  ft  #0  for 
(5.5,55),  a  contradiction.  Thus,  ?J[B^7  J-°,  B12  has  a  left  inverse,  and  the 
lemma  is  proved.  □ 

Before  a  minimal  left  inverse  for  (5.5,55)  can  be  given,  the 
following  two  results  are  needed. 

Lemma  lx  [2] 

The  pair  (K2,A,,2)  is  completely  observable  where  and  A^2  are 
defined  as  in  (5.5,55). 

Proo f ;  If  the  pair  (K2,A22)  is  not  completely  observable,  there  exists  a 
A22-invariant  subspace  Y  in  .  If  W  is  a  basis  matrix  for  Y,  this 

implies  that 


A22W  -  WQ 

for  some  matrix  Q  and 


Define 


k2w  =  o 


0 


w 


(5.5,56) 


(5.5,57) 


W  * 


(5.5,58) 


11J 


so  Chat  W  has  n  rows.  From  the  special  structure  of  (A+BK)  in  (5.5,55), 
it  foLLows  that 


(A+BK)W  -  WQ 

-  AW  +  BKW 

-  AW  +  B1C,W 


(5.5,59) 


-  AW 


From  the  first  and  last  lines  of  (5.5,59) 

AW  -  WQ 

which  shows  that  ¥  is  A-invariant.  It  also  follows  that 

(C+DK)W  -  0 

-  CW  +  DKW 

-  CW+DK^W 

-  CW 


(5.5,60) 


(5.5,61) 


From  the  first  and  last  lines  of  (5.5,61) 

CE  -  0  (5.5,62) 

which  shows  that  ¥  is  in  the  null  space  of  C.  It  follows  that  ¥  is  an 
unobservable  subspace  of  (A,B,C,D)  which  contradicts  the  observability 
assumption  on  (C,A).  Hence,  (l^.A^-,)  is  a  completely  observable  pair.  Q 
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Theorem  22  [2] 

Denote  the  characteristic  polynomial  of  A^,0  by  a^(s).  For  any 

A 

arbitrary  left  inverse  of  (5.5,55)  with  xQ  ■  0,  let  a(s)  be  the  characteristic 

a  m 

polynomial  of  A  in  representation  (5.5,37).  Then  °L(s)  divides  a(s). 

Proof:  Let  of  the  system  in  (5.5,55).  Since  this  system  is  completely 

controllable,  there  exists  an  input  segment  Ujq  ^  ^  such  that  §  ■  x^  for 
some  finite  k.  Consider  the  input 


OsiSk-1 


ksi<tt 


where  K  is  the  feedback  matrix  in  (5.5,55).  For  i  ik, 
(5.5,55)  is  given  by 


(5.5,63) 


the  solution  of 


Xi+1  *  (A  +  BK)xi 


y.  -  (C  +  DK)x. 


u .  *  Kx. 
i  l 


with  the  initial  condition 


(5.5,64) 


(5.5,65) 


since  §  €  JC  .  The  special  form  of  the  system  given  by  (5.5,55)  shows  that 
for  i  2  k 
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i-k~ 

ui  =  K2A22  *k 

(5.5,66) 

yi  *  0 


Now  this  u^  must  be  produced  as  an  output  from  any  left  inverse  of  (5.5,55) 

with  yi  as  an  input.  Let  this  inverse  have  a  representation  (5.5,37)  with 

*  »  —  _  „  „ 

A:  and  &  and  let  A  and  C  be  the  observable  subsystem  of  (A,C).  Since  y.  *0 

for  iik,  from  (5.5,37)  it  follows  that 

ui=CA1"kwk  (5.5,67) 

ic 

Now  was  an  arbitrary  vector  of  £  .  This  implies,  from  equations 

(5.5,66)  and  (5.5,67),  that  there  exists  a  matrix  W  such  that 


^2A22^  =  ®  ^i_kw 


(5.5,68) 


for  isk.  It  follows  that 


K2 

c 

K2A22 

C  A 

• 

• 

sb 

• 

*2A22 

• 

clj-k 

w-q2w 


(5.5,69) 


where  (j-k)  amax(d(A22) ,d(A)) .  The  pair  (K2,A22)  is  observable  by  Lemma  11 

* 

so  that  P[Q^J  *s  where  d(£  )=s.  From  (5.5,69),  it  follows  that  P[W]  =s. 

As  the  pair  (C,  A)  is  completely  observable,  Q2  has  full  column  rank  and, 
hence,  a  left  inverse,  Q+.  Thus,  (5.5,69)  yields 
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“  w  (5.5,70) 

Extending  (5.5,69)  one  more  step, 

Q1A22  =  Q2AW  (5.5,71) 

Multiply  on  the  left  by  Q*  and  use  (5.5,70): 

WA-22  -  AW  (5.5,72) 

This  expression  shows  that  y*ft[W]  is  A  -  invariant  and  nat 


A 


22 


«  A 


V 


(5.5,73) 


Thus,  the  characteristic  polynomial  of  A00,  °L(s)  divides  the  characteristic 
polynomial  of  A,  Q(s).  Since  (C,A)  is  the  observable  subsystem  of  (C,A), 
it  follows  that  &(s)  divides  Q'(s).  Hence,  QL(s)  divides  »(s)  and  the  theorem 
is  proven.  □ 

This  theorem  provides  a  lower  bound  for  the  dimension  of  a 
minimal  left  inverse.  This  bound  is  d(£  ). 

The  next  theorem  gives  an  explicit  representation  to  a  minimal 
left  inverse. 

Theorem  23 

Suppose  that  a  system  (A,B,C,D)  with  zero  initial  conditions  has 
a  left  inverse.  Then  a  minimal  order  left  inverse  of  order  d(£*)  is  given 
by 


Xi+1  “  A22Xi  +  A2lN(p)'y:i  +  tB2l  B22^N(p)yi 


mi 


I 

I 
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I . 


li 

[ 

E 

1 

I 


UL  “  K2xi  +  KiN(^P^t  +  N (p)y L 

where  the  notation  Is  taken  from  (5.5,55),  the  polynomial,  matrices,  N(s) 
and  N(.s),  are  defined  by  Lemma  10,  and  p  Is  the  unit  delay. 

ft 

Proof:  It  follows  from  Theorem  22  that  d(£  )  Is  the  minimum  order  of  any 

left  inverse.  For  an  input  sequence  define 

Vi  “  Ul  “  **L  (5.5,74) 

where  K  is  a  feedback  matrix  such  that  the  system  can  be  placed  in  the 
form  of  (5.5,55).  Place  the  system  in  the  form  of  (5.5,55)  and  adopt  the 
notation  of  (5.5,55).  From  Lemma  10,  there  exist  matrices  N(s)  and  N(s) 
such  chat 


*t  -  N(p)yt 
VL  *  N(p)yL 


(5.5,75) 


Substituting  (5.5,75)  into  (5.5,55c)  yields 


Ui  "  KlXl  *  K^xi+  vi 

-  K;,xi  +  KjN(p)yt  +  N(p)yt 


(5.5,7b) 


where  p  is  the  unit  delay.  From  (5.5,55a),  it  is  seen  that  x^  satisfies 


Xi+l  "  *22Xi  +  A2lXi  +  B2Vi 


A22Xi  +  A2lN^p)yt  +  B2N(P'yi  (5.5,77) 


x0  "  0 


J 
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Equations  (5.5, 76)  and  (5. 5, 77)  are  the  desired  result.  □ 

Consider  now  the  construction  of  a  right  inverses.  Suppose  a 

system  (5.5,1)  with  zero  initial  conditions  has  a  right  inverse.  In  general, 
★ 

ft  r  0  so  that  by  Theorem  11,  there  exists  a  feedback  matrix  K  such  that 
(5.5,1)  has  the  form 


*i+1  “  (A  +  BK)xi  +  Bvt  ■ 


An  0 


A2i  A22 


0 


B3  B4 


(5.5,78a) 


yi  -  (C+DK)xi+Dv.  =  [CL  0] 


u 

0 


— _ — 

—  — 

LI 

*i 

+  [D  0] 

(5.5,78b) 

_xi_ 

LJ 

ui  -  Kxi  +  vi  ■  K1xi  +  K2xt  +  vi 


(5.5,78c) 


(It  is  assumed  that  the  necessary  change  of  bases  has  been  done). 
Lemma  12 

The  system  ((A  +  BK), 


L  B3_J 

system  (5.5,76)  has  a  right  inverse. 


B, 


,  (C+DK),D)  has  a  right  inverse,  if  the 


Proof:  The  special  structure  of  the  system  in  (5.5,78)  permits  the  following 
calculation: 


CA  B  •  [C^  0] 


A11 

I 

o 

i 

B1 

0 

a2i 

A22 

B3 

B4 

[ c i A i i B i  01(5.5,79) 


This  shows  Chat  the  block  columns  in  the  matrix  Qjq  (5.5,4)  have  at 
most  linearly  independent  columns  where  is  the  number  of  columns 
in  B^.  From  equation  (5.5,19),  it  is  seen  that  P [ J  is  the  number  of 
linearly  independent  columns  in  that  column  of  blocks,  i.e.  P[F^]  for 
(5.5,78)  is  at  most  m.  .  Since  (5.5,78)  has  a  right  inverse,  it  follows 

Pi 

from  Theorem  19  that  m^"r.  Hence,  the  system  ((A  +  BK),  ,  (C+DK),D) 

B3 

has  a  right  inverse.  □ 


Lemma  13 


The  system  ((A+BK),  „  ,  (C+DK),D),  obtained  from  (5.5,78),  has 

Hi 


a  left  inverse. 


Proof:  For  the  system  (5.5,78), 


ft  -  <  A„ 


(5.5,80) 


See  equation  (5.2,31).  Since  ft  is  maximal  for  (5.5,78)  ,  it  follows  that 

*  Pll  -  B~  _ 

ft  -0  for  ((A+BK),  ,  (C+DK),D).  By  Corollary  3,  ((A  +  BK),  ,D) 

l?3J  -BSJ 

has  a  left  inverse.  a 


Lemma  13  provides  a  method  for  constructing  a  right  inverse  for 


(5.5,78).  If  a  system  has  both  a  left  and  right  inverse,  then  these  in- 

P 

verses  are  the  same.  By  Lemma's  12  and  13,  this  is  true  for  ((A+BK),„ 

lb 

(C+DK),D)  and  its  inverse  is  given  by  Lemma  10. 


’L?J’ 


To  obtain  the  inverse  of  (5.5,78),  form  two  matrices,  S^  and  S^ , 


such  that 
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Now  (5.5,78c)  can  be  written  as 


(5.5, Sla) 


(5.5, 81b) 


ui  ’  Kxi  +  Vi  *  Klxi  +  K2xi  +  Vi  +  s2^i  (5.5,  82) 

Substituting  from  Lemma  10, 

ui  “K1N(p)^t  +  K2xi  +  S1N(p)y.+S2vi  (5.5,  83) 

where  (from  5.5,  78)  xi  satisfies 

Vl  "  A22xi  +  A21xi  +  Vi  +  Vi  (5 . . . 5,  84) 

mA22Xi+A21N(p)y'i  +B3N(P)yi+B4v. 

again  by  Lemma  10.  Equations  (5.5,83)  and  (5.5,84)  constitute  a  right  in¬ 
verse  for  (5.5,78).  These  results  are  summarized  in  the  next  theorem. 
Theorem  24 

Suppose  that  a  system  (5.5,1)  with  zero  initial  condition  has  a 
right  inverse.  Then  a  right  inverse  is 

*1+1  '  A22*i  +  A21N  +  B35<P>  *1 +  Vi 
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ui  -K2xi +KiN<P)yki  +siN(p)yi  +  S2vi 

where  the  matrices  N(s)  and  N(s)  are  given  by  Lemma  10,  p  is  the  unit 
delay,  and  is  an  arbitrary  input  sequence.  0 

The  construction  of  the  left  and  right  inverses  of  a  system 
(5.5,1)  as  in  Theorem's  23  and  24,  clarifies  the  role  of  invariant  zeros  in 
inverse  systems.  Consider  left  inverses.  The  poles  of  the  minimal  left  in¬ 
verse  are  the  eigenvalues  of  A2 , .  Recall  that  A.,.,  was  obtained  from  canonic 

•k 

form  (5.2,24)  and  that  ft  *0  for  a  left  invertible  system.  By  Theorem  14, 
the  invariant  zeros  are  just  the  eigenvalues  of  A,,2>  invariant  zeros 

of  (5.5,1)  are  the  poles  of  its  minimal  left  inverse,  if  it  exists.  By 
Theorem  22,  the  invariant  zeros  are  contained  among  the  poles  of  any  other 
left  inverse. 

Now  suppose  that  the  system  (5.5,1)  has  a  right  inverse  and  it  is 
given  by  Theorem  24.  As  noted  in  Section  5.3,  the  invariant  zeros  are  the 
uncontrollable  modes  of  (A22>B^).  Again,  the  invariant  zeros  are  found 
among  the  poles  of  the  inverse  system.  This  is  a  generalization  of  the  well 
known  classical  result.  These  results  are  summarized  in  the  following 
theorem. 

Theorem  25 

Consider  the  system  (5.5,1)  with  zero  initial  conditions. 

a)  If  (5.5,1)  has  a  left  inverse,  then  the  invariant  zeros  are 
contained  among  the  poles  of  the  left  inverse  system.  If  this  system  has 
minimal  order  then  the  invariant  zeros  are  exactly  the  poles  of  the  left  inverse. 
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b)  If  (5.5,1)  has  a  right  inverse  system  and  is  constructed 
according  to  Theorem  24,  then  the  invariant  zeros  appear  among  the  poles 
of  the  right  inverse  system.  q 

These  results  were  noted  by  Bengtsson  [2]  for  left  inverses  of 
systems  with  D  *  0. 

Actually,  a  slightly  stronger  statement  can  be  made  for  systems 
with  a  left  inverse.  Theorem  14  states  that  the  invariant  polynomials  of 
the  system  matrix  for  (5.5,1)  are  the  same  as  the  invariant  polynomials  of 
A22'  In  Section  4.2,  it  was  noted  that  the  invariant  zeros  determine  whether 
a  system  has  simple  or  non-simple  structure.  This  structure  is  determined 
by  the  invariant  polynomials,  and  this  structure  carries  over  into  the  pole 
configuration  of  the  left  inverse  system. 

Of  course,  the  stability  of  the  inverse  system  is  also  governed 

★ 

by  the  location  of  the  invariant  zeros.  For  a  right  inverse,  if  ft  r 0 
some  of  the  poles  of  the  right  inverse  system  are  not  invariant  zeros.  It 
is  exactly  those  poles  which  are  controllable  from  the  arbitrary  input. 

Hence,  if  the  instability  of  a  right  inverse  is  not  due  to  a  right  half 
plane  invariant  zero,  the  system  can  be  stabilized  by  using  the  extra  free¬ 
dom  in  the  control.  These  stability  questions  of  inverse  systems  are  also 
addressed  in  [40],  although  the  role  of  invariant  zeros  is  not  explicit. 

This  reference  also  provides  a  method  for  constructing  a  reduced  order  in¬ 
verse  system  based  on  the  work  in  [34], 

Finally,  it  is  noted  that  in  [41]  it  is  shown  for  a  class  of  systems 
that  the  invariant  zeros  of  the  inverse  system  are  actually  the  poles  of  the 
original  system;  an  interesting  generalization  of  a  classical  result. 
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CHAPTER  6 

COMPUTATION  OF  ZEROS 

6.1.  Introduction 

Several  algorithms  have  appeared  in  the  literature  for  calcu¬ 
lating  zeros  ([ 22] , [3 1] , [ 32] , [42] , [43) , [44] )  based  on  the  properties  of 
zeros  described  in  Chapter  5.  These  algorithms  provide  efficient  methods 
for  calculating  zeros  or  they  are  readily  adaptable  to  a  digital  computer. 
All  of  these  algorithms  utilize  state  space  representations  of  systems 
and  they  all  calculate  invariant  zeros.  However,  there  are  certain  limita¬ 
tions  and  distinctions  in  each. 

Section  2  discusses  the  generalized  eigenvalue  method.  This 
algorithm  requires  that  the  system  be  nondegenerate.  It  can  also  be  used 
to  calculate  decoupling  zeros. 

Section  3  presents  two  methods  based  on  high  gain  feedback.  The 
first  is  a  straightforward  application  of  Theorem  17.  It  requires  that 
the  system  be  nondegenerate.  The  second  algorithm  applies  the  idea  of 
high  gain  feedback  to  obtain  a  geometric  method  of  calculating  invariant 
zeros . 

The  Davison-Wang  method  is  introduced  in  Section  4.  The 
invariant  zeros  are  observed  to  be  the  limiting  positions  of  the  eigen¬ 
values  of  the  matrix 


A  YB 
_C  yd. 


(6.1,1) 


where  Y  is  a  scalar  that  approaches  infinity.  The  system  is  required  to 
be  nondegenerate. 
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Section  5  discusses  a  certain  class  of  systems  for  which  a  right 
inverse  system  is  easily  constructed.  For  these  systems,  the  results  of 
Section  5.5  on  inverse  systems  are  considerably  simplified.  The  invariant 
zeros  are  then  calculated  as  the  uncontrollable  modes  of  the  right  inverse 
system. 


6.2.  Generalized  Eigenvalue  Approach 


As  a  motivation  for  this  algorithm,  consider  a  nondegenerate  state 
space  system  with  an  equal  number  of  inputs  and  outputs.  The  invariant 
zeros  are  given  by 


sI-A  -B 

det  =  0. 

-  c  D. 

Because  P(s)  loses  rank  when  (6.2.1)  is  satisfied,  it  follows  that 


(6.2,1) 


(6.2,2) 


s  A  x  =  B  x  . 


(6.2,3) 


This  last  equation  defines  a  generalized  eigenvalue  problem.  The 
numerical  solutions  to  this  problem  have  been  studied  and  the  results  can 


be  applied  here  to  obtain  the  invariant  zeros  of  a  system  defined  by  the 
system  matrix  in  (6.2,1).  Patel  [42]  first  suggested  this  approach  and 
proposed  a  numerical  algorithm  based  on  [45] .  Since  then  it  has  been 
suggested  that  the  QZ-algorithm  be  used  to  solve  (6.2,3)  [43].  The 
numerical  advantages  of  the  QZ-algorithm  are  emphasized  in  this  reference. 
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This  analysis  has  assumed  Chat  the  system  matrix  in  (6.2,1)  has 
normal  rank  n+m  as  this  is  a  requirement  for  the  numerical  algorithm  given 
in  [45].  If  this  algorithm  is  used,  this  requirement  must  be  checked 
beforehand.  Alternatively,  if  the  QZ-algorithm  is  applied  to  (6.2,3), 
it  will  automatically  detect  degeneracy.  In  either  case,  this  approach 
will  not  identify  invariant  zeros  (Definition  5)  of  a  degenerate  system. 
However,  it  is  naturally  suited  to  calculate  zeros  defined  by  Definition  6. 
See  comments  in  Section  4.2. 

In  [43] ,  this  method  is  extended  to  systems  with  an  unequal 
number  of  inputs  and  outputs  by  "squaring  up"  the  system.  Suppose  that 
m>r.  Generate  two  pseudo-random  matrices  and  such  that 

0 

n 

0  0  ,  B  ■ 

0  0 

where  B  is  a  (n+m) x  (n+m)  matrix.  It  can  be  seen,  via  Theorem  4,  that  the 
invariant  zeros  of  the  original  system  are  contained  in  the  invariant 
zeros  of  the  new  system.  To  find  the  invariant  zeros  of  the  original 
system,  first  calculate  the  invariant  zeros  of  (6.2,4).  Then  generate  two 
new  matrices,  E2  and  F0 ,  and  calculate  the  ’’nvariant  zeros  again, 
replacing  E^  and  F^  by  E,,  and  F, ,  respectively.  Then  the  invariant  zeros 
of  the  original  system  are  almost  surely  contained  in  the  intersection  of 
these  two  sets  of  invariant  zeros.  If  r>m,  this  problem  can  be  reformu¬ 
lated  in  an  analogous  way,  or  the  dual  system  can  be  considered. 

The  decoupling  zeros  can  be  calculated  using  this  generalized 
eigenvalue  technique,  too.  Recall  that  the  input  decoupling  zeros  are 
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chose  complex  numbers  s  such  that 

s[  1  ;  0]  -  [A  ;  B] 


(6.2.5) 


loses  rank  (Theorem  7).  Augmenting  these  matrices  by  pseudo-random 
matrices  and  of  the  appropriate  size  gives: 


s 


(6.2,6) 


This  again  can  be  used  as  a  generalized  eigenvalue  problem  and  solved  as 
in  the  case  of  invariant  zeros  with  unequal  number  of  inputs  and  outputs. 
The  output  decoupling  zeros  are  handled  in  a  similar  fashion. 

Of  course,  once  the  invariant  zeros  and  input  and  output 
decoupling  zeros  are  known,  the  transmission  zeros  and  the  system  zeros 
can  be  recovered  using  the  relationships  in  Theorem  9. 

More  recently,  this  problem  has  been  cast  in  a  more  general 
theoretical  framework  [9],  It  has  been  noted  that 


P(s) 


sI-A  -B 

C  D 


(6.2,7) 


is  a  singular  pencil  of  matrices  [11]  and  that  its  finite  divisors  are 
just  the  invariant  zeros.  The  solutions  to  [6.2,2)  are  just  the  roots  of 
these  finite  divisors.  However,  in  general,  the  number  of  solutions  to 
(6.2,3)  is  less  than  the  dimension  of  A.  Therefore,  it  is  proposed  in 
[46]  that  a  sequence  of  transformations  be  applied  to  (6.2,7)  that  will 
reduce  P(s)  to  a  block  canonic  form.  One  of  the  blocks  will  contain 
exactly  the  invariant  zeros  (finite  divisors)  of  (6.2,7).  Then  the  QZ- 
algorithm  can  be  applied  with  greater  efficiency.  It  turns  out  that  the 


necessary  transformations  are  exactly  the  structure  algorithm  [34]  as 
modified  in  [40].  Furthermore,  these  transformations  will  reduce  a 
nonsquare  P(s)  to  a  square  matrix  P(s)  so  that  the  CZ-algorithm  can  be 
applied  directly.  The  details  are  found  in  [46]. 

6.3.  High  Gain  Feedback 

Several  algorithms  for  computing  invariant  zeros  have  appeared 
that  are  based  on  the  idea  that  zeros  are  the  limiting  positions  of  the 
systems  poles  under  high  gain  feedback  ([ 32] , [ 31] , [47] , [48] ) . 

The  first  algorithm  is  obtained  directly  from  Theorem  17.  It 
requires  that  the  system  be  nondegenerate  so  that  the  algorithm  first 
checks  this  condition.  If  the  system  is  nondegenerate,  then  the  invariant 
zeros  are  contained  among  the  finite  eigenvalues  of 

A  +  B (~  Ir-KD)-1KC  (6.3,1) 

as  p-00.  Here,  K  is  an  arbitrary  mx  r  matrix  with  rank  min(m,r).  If 
m  =  r,  K  can  be  chosen  to  be  the  identity  matrix.  Then  the  invariant 
zeros  coincide  with  the  finite  eigenvalues  of  (6.3,1).  If  m^r,  then  for 
"almost  all"  K  the  invariant  zeros  will  be  contained  in  the  finite 
eigenvalues  of  (6.3,1).  Hence,  choose  a  K  and  a  suitably  large  P  and 
calculate  the  eigenvalues  of  (6.3,1).  If  m#r,  the  finite  eigenvalues 
must  also  satisfy  the  definition  of  an  invariant  zero. 

Again,  the  nondegeneracy  condition  implies  that  for  systems 
for  which  this  algorithm  is  applicable,  Definitions  5  and  6  define  the 
same  set  of  zeros.  It  will  not  identify  invariant  zeros  (Definition  5) 
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for  degenerate  systems.  Hence,  it  is  more  suited  for  use  with  Definition 
6.  See  [32]  for  details. 

The  second  algorithm,  however,  does  calculate  invariant  zeros 
(Definition  5)  for  degenerate  systems.  Suppose  that  D *  0  so  that  (6.3.1) 
reduces  to 

A  +  PBKC.  (6.3,2) 

If  m * r  (number  of  inputs  equals  the  number  of  outputs),  then  it  is 
shown  in  [31]  that  the  invariant  zeros  tend  to*the  roots  of  the  equation 

det[ sNM  -  NAM]  -  0  (6.3,3) 


where  N  and  M  are  the  full -rank  left  and  right  annihilators  of  B  and  C, 
respectively,  such  that 


NB  *  0 
CM  *  0 


(6.3,4) 


when  p-"*.  Furthermore,  if  the  product  CB  has  full  rank,  it  is  shown  [31] 
that  N  and  M  can  always  be  selected  so  that 


NM  -  I.  (6.3.5) 

In  this  case,  the  invariant  zeros  are  just  the  eigenvalues  of  NAM.  If 
CB  does  not  have  full  rank  or  m^r  then  suitable  modifications  are  made 
in  the  algorithm.  The  details  are  in  [31]. 

In  [47] ,  the  following  interpretation  is  given  to  the  high  gain 
feedback  system  (6.3,2).  Suppose  that  m*r  and  choose  K  to  be  the 
identity  matrix.  Then  as  p— <*>  some  of  the  eigenvalues  of  (6.3,2)  go  to 
infinity  while  others  tend  to  finite  values  (the  invariant  zeros).  Hence, 
there  is  a  natural  separation  of  the  poles  of  the  closed-loop  system  into 


slow  and  fast  modes  as  p  increases.  This  heuristic  argument  shows  that 


this  system  is  naturally  suited  for  singular  perturbation  analysis  which 
is  done  in  [47]  and  [49] .  The  results  are  essentially  the  same  as  for  the 
NAM  algorithm,  however,  this  analysis  does  produce  a  specific  procedure  to 
calculate  the  matrices  N  and  M  when  CB  has  full  rank.  See  [47]  or  [49] 
for  details. 

Another  specific  construction  of  the  matrices  N  and  M  is  given  in 
[48]  for  systems  with  equal  numbers  of  inputs  and  outputs.  There  the  NAM 
algorithm  is  extended  to  systems  with  D/0.  This  provides  a  method  for 
calculating  system  zeros  as  well  as  invariant  zeros.  Recall  that  the 
set  of  system  zeros  is  the  intersection  of  the  sets  of  invariant  zeros 
obtained  from  certain  subsystems  of  the  original  system.  (See  the  remark 
in  Section  4.4  following  Example  14.)  So  simply  apply  the  algorithm  in 
[48]  (or  any  other  algorithm  for  computing  invariant  zeros)  to  all 
subsystems  of  the  form  (4.4,11).  The  system  zeros  are  then  the  invariant 
zeros  common  to  all  these  subsystems. 


6.4.  The  Davison-Wang  Method 

The  following  method  was  proposed  by  Davison  and  Wang  [22] . 
Consider  a  system  with  an  equal  number  of  inputs  and  outputs  which  is 
nondegenerate.  Define  the  matrix  S (y)  as 

(a  yb] 


S(Y)  = 


(6.4,1) 


C  YDJ 


It  is  shown  that  as  Y  becomes  arbitrarily  large,  p  eigenvalues  of  S  (y) 
become  arbitrarily  close  to  the  invariant  zeros  of  the  original  system. 
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To  apply  this  algorithm,  first  determine  if  the  system  is  degenerate.  If 
it  is  not,  then  choose  a  suitably  large  Y  and  calculate  the  eigenvalues  of 
S(Y).  The  invariant  zeros  are  those  eigenvalues  of  S(Y)  which  satisfy  the 
definition  of  invariant  zeros.  Under  the  assumption  of  nondegeneracy, 
this  could  be  either  Definition  5  or  6.  However,  as  with  the  generalized 
eigenvalue  approach,  this  algorithm  is  naturally  suited  to  Definition  6. 

If  the  system  has  an  unequal  number  of  inputs  and  outputs,  then 
S(Y)  can  be  augmented  by  pseudo-random  matrices  and  yF^  to  make  S(Y) 
square.  The  procedure  is  exactly  the  same  as  discussed  for  the  generalized 
eigenvalue  problem,  Section  6.2. 

6.5.  Calculating  Invariant  Zeros  from  Inverse  Systems 

The  final  method  for  computing  invariant  zeros  is  found  in  [44] . 
It  is  based  on  the  results  of  Section  5.5  in  which  invariant  zeros  are 
related  to  the  poles  of  the  inverse  system.  Specifically,  suppose  that 
a  left  or  right  inverse  system,  as  given  in  Theorems  23  or  24,  can  be 
constructed.  Then  the  invariant  zeros  could  be  computed  from  the  dynamics 
of  the  system.  Of  course,  this  construction  can  always  be  carried  out  as 
suggested  by  the  proofs  of  the  theorems.  However,  in  some  cases  this 
construction  is  accomplished  with  much  less  effort. 

First  consider  a  system 

xi+l  "  ^*i  +  Bui  (6*5,  la) 

+  Dut  (6.5,  lb) 

in  which  mir  and  D  has  full  rank.  From  Algorithm  1  (Section  5.2),  it  is 
seen  that  p [ F^ ]  ■  r .  Therefore,  by  Theorem  19,  (6.5,1)  has  a  right  inverse. 
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If,  in  addition,  r«m  then  p  [  ]  ■  m  and  (6.5,1)  also  has  a  left  inverse. 

Using  the  notation  of  Lemma  10  and  Theorem  23,  the  fact  that  p[D]  ■  r 
implies 

7i  ■  yt  (6.5,2) 

y t  -  o. 

Now  by  Algorithm  I,  if  P [ F^]  “r,  it  follows  that  M^O.  Hence, 

-TJIMq]  =2  (6.5,3) 

where  ^  is  the  state  space.  By  Lemma  1  and  Theorem  12 

XQ  -  Zl  »  £*  *  Z.  (6.5,4) 

Therefore,  to  find  K,  the  feedback  matrix  which  puts  (6.5,1)  into  canonic 
form  (5.2,24),  it  is  necessary  to  solve 

C  +  DK  -  0.  (6.5,  5) 

First  assume  m-r.  Then  D  1  exists  and  (6.5,5)  yields 

K  -  -D‘lC.  (6.5,6) 

From  canonic  form  (5.2,24) 

-  A  +  BK  *  A-BD_1C.  (6.5,7) 

Under  the  assumptions  above,  from  Lemma  10  it  follows  that 

N(p)  -  D'1  (6.5,8) 

N (p)  -  0. 

Now,  Theorem  23  gives  the  left  inverse  for  (6.6,1),  assuming  that  m*r  and 
D  has  full  rank, 
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x 


i+1 


u 


i 


(A-BD*1C)xi  +  BD'^i 
-D_1Cxi  +  D_1yi. 


(6.5,9) 


Hence,  the  invariant  zeros  are  the  eigenvalues  of  (A-BD  C)  which  is  one 
of  the  results  given  in  [44]  . 

Now  suppose  that  m>r  and  D  has  full  rank.  Equations  (6.5,2)- 
(6.5,4)  still  hold  and  the  required  matrix  K  must  still  solve  (6.5,5). 
The  procedure  given  in  [44]  is  this:  Write  (6.5,5)  as 


Cxi  +  Du^  “  yi.  (6.5,10) 

A  general  solution  to  this  equation  is 

ui  *  Ulm-D+D)L-D+C]x.  +  D+yt  (6.5,11) 

where  L  is  an  arbitrary  mxn  matrix  to  provide  degrees  of  freedom  and 

D+  -  D+(DD+)'1.  (6.5,12) 


Substitute  for  in  (6.5,1a) 

Xi+1  =  (  (A-BD+C)  +B(Im-D+D)Llxi  +  BD+yi . 

To  find  the  invariant  zeros,  calculate  the  eigenvalues  of 


(A-BD+C)  +  B(Im-D+D)L 


(6.5,13) 


(6.5,14) 


for  two  values  of  L,  say  zero  and  L  "larger."  The  invariant  zeros  are 
the  eigenvalues  of  (6.5,14)  which  are  the  same  for  both  values  of  L. 

To  see  how  this  fits  into  the  theory  of  inverse  systems,  perform 
an  input  space  transformation  to  bring  (6.5,1)  into  the  form 


i+1 

*  Axi  +  [Bj 

b2] 

*=  Cx.  +  [D 

0] 

“i 


(6.5,15) 


L 


where  D  is  a  rxr  nonsingular  matrix  and  B  is  partitioned  compatibly. 
A  straightforward  calculation  shows  that  (6.5,12)  yields 


(6.5,16) 


In  particular,  K * -D+,  where  D+  is  given  in  (6.5,16),  solves  (6.6,6). 


Lemma  10  again  yields 


N (p)  *  D  1 
N(p)  -  0. 


(6.5,17) 


Now  Theorem  24  gives  a  right  inverse  for  (6.6,15) 

xi+l  *  (A-BDTC)xi  +  Bjlflyt  +  B,u.. 


(6.5,18) 


Again,  a  straightforward  calculation  using  (6.5,15)  and  (6.5,16)  shows 


fl  ol 

B(I  -D  D)L  -  [B,  BJ  I  -  r  L 

m  1  2  L  L0  0| 


[o  b21l. 


(6.5,19) 


So,  L  can  be  interpreted  as  a  feedback  matrix  between  (A-BD  C)  and  B,,  in 


(6.5,18).  Since  £*  «*  for  (6.5,15),  it  follows  from  (5.2,31)  that 


ft*  *  (  (A-BD+  C)  ftlB2l  >  . 


(6.5,  20) 


Then  the  invariant  zeros  of  (6.5,15)  are  the  uncontrollable  modes  of 
( (A-BD+C) ,B2)  by  Theorem  14.  By  choosing  those  eigenvalues  of 
( (A-BD^C)  +B(Im-D+D)L)  which  are  invariant  to  choices  of  L,  the  eigenvalues 
of  the  uncontrollable  modes  of  ((A-BD+C) ,B,)  are  selected.  These  are 
exactly  the  invariant  zeros  of  (6.5,15). 
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Finally,  suppose  r>m  for  (6.5,1)  and  D  has  full  rank.  Then, 
the  invariant  zeros  can  be  calculated  by  applying  the  procedure  of 
(6 .5  ,10)- (6 .5 ,14)  to  the  dual  system  (A^,CT,B^,DT) . 

The  method  given  in  [44]  can  be  applied  to  systems  (6.5,1)  for 
which  Ds  0.  It  is  required  that  B  and  C  have  full  rank  and  that  CB  have 
full  rank. 

Define  a  similaritv  transformation 


*1  '  **1 


where 


n-r 


lcx  c2l . 


(6.5,21) 


(6.5,22) 


Since  T  must  be  nonsingular,  must  be  nonsingular.  As  C  has  full  rank, 
this  can  always  be  accomplished  by  permuting  the  state  variables  if 
necessary.  The  application  of  (6.5,21)  to  (6.5,1)  yields 


“i+1 


A11  A12 


V21  h.2 


2i  + 


(6.5,23) 


y±  -  I  ir  o]  zi . 

In  [50],  it  is  shown  that  the  transfer  function  of  (6.5,23)  can  be  written 
as 


G(s)  =  Q~ 1 (s )W  (s) 


where 


-1, 


Q (s  )  *  si  -  A-^  -  (^^""^22)  A2i 


W(s)  *  Bj_  +  A12(sI-A22)  Bj. 


(6.5,24) 

(6.5,25a) 

(6.5,25b) 


1 
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The  system  matrix  of  (6.5,23)  is  given  by 


P(s) 


sT-An  -A12 


_A21 

I 


SI’A22 

0 


'  B„ 


(6.5,26) 


Using  elementary  row  and  column  operations  it  can  be  brought  into  the  form 


I 

r 


0 


0 


0 


•  (6.5,27) 


By  Definition  5,  the  invariant  zeros  of  (6.5,23)  are  the  roots  of  the 
invariant  polynomials  of  the  matrix  (6.5,27).  On  the  other  hand,  the 
lower  right  hand  corner  of  (6.5,27)  can  be  interpreted  as  the  system 
matrix  of  the  system  (A22,B2,A12,B1^ ’  w^ose  transfer  function  is  given  in 
(6.5,25b).  Hence,  the  invariant  zeros  of  (6.5,23)  are  the  same  as  the 
invariant  zeros  of  the  reduced  order  system  (A^B,,  >A12>BL)  and  the  trans 
mission  zeros  of  (6.5,23)  are  the  same  as  the  reduced  order  system  whose 
transfer  function  matrix  is  given  by  (6.5,25b)  (Theorem  6). 

To  actually  calculate  the  invariant  zeros  of  >  B2  ,A12  ’ 
note  that  the  transformation  in  (6.2,22)  yields 


CB 

TB  * 

.2 

Since  CB  was  assumed  to  have  full  rank,  B^  has  full  rank.  Therefore,  the 
invariant  zeros  of  (A22>B->>  A^2,B1^  can  calculated  using  (6.5,7),  the 
procedure  given  in  (6 .5 , 10) - (6 . 5 , 14) ,  or  this  procedure  applied  to  the  dual 
system  (A22 ’^12  * ) Bl^ * 


(6.5,28) 
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Mor,  recently,  thi8  procedure  £ot  ^  of 

extended  to  the  cl.ee  of  all  tnverclbie  systems.  rhle 


was  accomplished  by  algorithm  „hlch  u  , 
Structure  Algorithm  [34];  i.e.  by  augmenting 
full  rank.  The  details  are  in  [51], 


variation  of  Silverman's 
the  D  matrix  to  make  it  have 
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CHAPTER  7 

A  NEW  ALGORITHM  FOR  CALCULATING  INVARIANT  ZEROS 

7.1.  Introduction 

This  chapter  presents  an  algorithm  for  calculating  invariant 
zeros.  The  algorithm  is  based  on  the  geometrical  properties  of  linear  time- 
invariant  systems;  specifically.it  calculates  canonic  form  (5.2.24)  which 
displays  explicitly  £  *,  the  maximal  null-output  (A, B)-invariant  subspace. 
Once  this  canonic  form  is  available,  the  invariant  zeros  are  easily 
calculated  via  Theorem  14. 

The  algorithm  is  developed  using  the  geometrical  techniques 
introduced  in  Section  5.2.  The  algorithm  proceeds  by  identifying  the 
subspaces  X^  for  i *  0,1,2,....  This  is  accomplished  by  applying  a 
sequence  of  similarity  transformations  and  feedback  matrices  to  the  system 
(A,B,C,D).  The  algorithm  terminates  when  ^  “^i+i  and  c^e  system  is 
in  canonic  form  (5.2.24).  This  form  then  can  be  further  reduced  to  canonic 
form  (5.2.38). 

Section  2developes  the  algorithm,  called  Algorithm  II,  for  calcu¬ 
lating  canonic  form  (5.2.24).  Algorithm  II  is  closely  related  to  Wonham's 
vector  space  algorithm  for  calculating  £*  [ 12 1 ,  Silverman's  structure  a 
algorithm  [34],  and  Algorithm  I  (Section  5.2).  This  relationship  is 
discussed  in  detail.  Other  methods  for  calculating  X*  are  also  briefly 
discussed . 


Section  3  explains  how  the  canonic  form  produced  can  be  used  to 
calculate  invariant  zeros.  Since  Algorithm  II  can  be  used  on  any  system 


( A ,  B ,  C ,  D ) ,  the  invariant  zeros  can  be  calculated  for  any  system  (A,B,C,D). 
This  method  of  computing  invariant  zeros  is  compared  to  other  methods  which 
were  discussed  in  Chapter  6.  Algorithm  II  is  then  summarized  and  several 
examples  are  given  to  illustrate  the  use  of  Algorithm  II. 


7.2.  An  Algorithm  for  Calculating  X* 

Theorem  11  states  that  any  system  (A,B,C,D)  can  be  placed  in 
this  canonic  form: 


xi+l 

A1 

0 

B1  0 

xi+l 

- 

*3 

A4 

B3  \ 

t-4 

( 

C1 

0 

D1  0 

xi 


(7.2.1) 


However,  this  theorem  does  not  provide  the  required  state  feedback  matrix, 
or  the  input  and  state  space  transformation  matrices.  The  purpose  of  the 
following  algorithm  is  to  provide  a  method  for  placing  a  completely 
arbitrary  system  (A,B,C,D)  in  canonic  form  (7.2.1). 

The  algorithm  proceeds  in  cycles;  the  i-th  cycle  identifying 
the  subspace  (Definition  11).  By  Lemma  1  when  ^  ■!*,  and 

the  system  will  be  in  canonic  form  (7.2.1).  Recall  that  Algorithm  I  also 
calculated  I*.  This  algorithm  is  based  on  its  properties. 

Algorithm  II 
1st  Cycle 

Consider  the  system  (A,B,C,D).  Since  it  is  desired  to  calculate 
Xl*  appl?  Steps  1  and  2  of  Algorithm  I  to  (A,B,C,D).  From  (5.2.15) 
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Vo 


(7.2.2) 


The  idea  is  to  express  this  calculation  in  terras  of  a  similarity  transforma¬ 
tion  on  the  original  system.  To  this  end,  find  a  nonsingular  input  space 
transformation  W  such  that 


DWo  -  (  (dV  0)  -  D1  (7.2.3)* 

where  (D*-)'  is  a  rxm^  matrix  which  has  full  column  rank.  Denote  B1,  ■  BW^ 
and  u^"Wou..  This  operation  corresponds  to  postmul tiplication  of  (7.2.2) 
by 


(7.2.4) 


Since  is  nonsingular,  the  row  rank  of  is  not  affected  and  the 
algorithm  is  unchanged. 

Now  define  a  nonsingular  output  space  transformation  S  such  chat 


S 

o 


(D1)' 


(7.2.5) 


-1 


where  D  is  a  x  m^  nonsingular  matrix.  Then 


y,  ■  S  v  * 
J  i  o^  i 


_1 

>’i 

-1 


S  Cx.  +  S  D1u]‘ 
o  i  o  i 


1  2  1 
■  C  +  D“u^ 


xi  + 


D1 


_1 

ui 


lui 


(7.2.6) 


In  this  section,  superscripts  indicate  indices,  not  exponents. 


■Hi 


■ 
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-a  ~i 

where  C  is  a  m^  x  n  matrix,  C  is  a  r^  x  n  matrix,  and  r^-r-m^.  Furthermore, 

y]  and  u}  are  partitioned  compatibly.  Clearly,  the  matrix  S  defined  in 
11  o 


(7.2.4)  could  be  used  in  (7.2.2).  In  that  case,  by  comparing  (7.2.2)  and 


(7.2.5),  it  follows  that  *5  . 


-1 


Remark  1 :  It  will  be  required  below  that  C  have  full  row  rank.  If  this 


is  not  true,  then  there  exists  a  nonsingular  matrix  such  that 


K'°l 


(c1)’ 


(7.2.7) 


where  (C*)  '  has  full  row  rank  and  E1  is  a  r.xr,  matrix.  Let 

o  11 


m. 


(7.2.8) 


and  define  S^“EoSo.  Clearly,  D“  is  not  affected  when  is  replaced  by 


S^.  Now  it  is  seen  that  M  (C*)'.  Notice  that 

7J[M1]  «  7UC1]  *  (C1) ']  . 


(7.2.9) 


Assume  now  that  C  has  full  rank  with  the  understanding  that  if  this  is 
not  true,  the  necessary  modifications  can  be  made  as  per  this  remark.  £ 
To  obtain  canonic  form  (7.2.1),  it  is  necessary  to  find  K*; 
that  is, a  feedback  matrix  that  will  make  S*  (A+BK*) -invariant.  This  is 
accomplished  by  finding  a  matrix  K^  in  each  cycle  and  calculating  A+BK,. 
When  the  algorithm  terminates,  K*  can  be  recovered,  if  necessary.  To  this 
end  let 


,—l .  - 1~1 
-(D)  C  . 


(7.2.10) 


Then  form  the  mx  n  matrix  K  as 

o 


u 


y 


a 
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(7.2.11) 


Compute 


A1  -  A  +  BLK 


(7.2.12) 


C  2  12  C1  D1  0  K 

,  “  C  «  C  +  DK  *  .  + 

c2  °  c1  0  0  0 


By  Lemma  2,  this  leaves  X*  unchanged.  (Note  that  C*"3^) 

Now  the  system  is  in  a  form  such  that  X^  can  be  displayed 
explicitly.  If  C^=0,  or  does  not  exist  (D  has  full  row  rank),  then  from 
(7.2.10)  and  (7.2.12),  it  follows  that  C2=0.  Since 


-  WC1]  -  ^[C2] 


(7.2.13) 


this  implies  that 


«  mMx]  -  % 


(7.2.14) 


where  £  is  the  whole  state  space.  Furthermore,  as  X  *7?[M  ]  =  M,  ]  =  £,  , 

o  o  i  1 

by  Lemma  1  it  is  seen  that 


X 

o 

112  2 

Hence,  the  system  (A  ,B  ,C  ,D  )  given  by 


(7.2.15) 


A1  -  [A1]  B1  -  [bJ:  bJ2] 

C2  -  [0]  D2  *  [D2l  0  ] 


(7.2.16) 


is  in  canonic  form  (7.2.1)  with  KqmK*.  Therefore,  the  algorithm  terminates. 

If  C^^O,  then  £Q^£^.  To  display  £^  explicitly,  define  the 
similarity  transformation 


z.  m  T  x 
i  o  i 


(7.2.17) 


where  Tq  is  the  nxn  matrix 
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~2  ~2 

C11  C12 


(7.2.18) 


-2  r  ~2  ~2  1 

C  "  [C11  C12 


and  Tq  is  nonsingular. 


Remark  2 :  If  Tq  is  to  be  nonsingular,  must  be  nonsingular.  By 

-2  ~2 
Remark  1,  C  has  full  row  rank.  Therefore,  if  is  singular,  it  can  be 

made  nonsingular  by  a  simple  permutation  of  the  state  variables.  This 

can  be.  expressed  as  a  similarity  transformation  of  the  type  (7.2.17), 

however  the  details  are  omitted  for  clarity  of  presentation.  0 


From  (7.2.5),  (7.2.17),  and  (7.2.18)  it  is  seen  that 


(7.2.19) 


2  1-12  1  32-1 

Calculate  A  'TAT  ,  B  =T  B  ,  and  C  =C  T  . 


The  system  matrices  now 


have  the  form 


2  2 
A11  A12 


B11  B12 
B21  B22 


(7.2.20) 


where  is  a  r^x  ^  matrix.  Since  Tq  is  nonsingular, it  follows  from 


(7.2.13)  that 


where 


fttMj.]  =  me2]  *  7?[c3] 


[  Ir  0]  . 

rl 


(7.2.21) 


(7.2.22) 


L43 


Now,  it  is  obvious  from  (7.2.20)  that 


TUMj]  -  fl[C  1  -  Xj  -  j  • 

xi 


(7.2.23) 


This  completes  the  1st  cycle  of  the  algorithm.  To  summarize 
the  development  thus  far,  motivate  the  next  step  of  the  algorithm,  and 
lend  further  insight  into  the  system  structure,  consider  again  (7.2.20). 
For  i  ■  0 


-1 

2  -1 

.2  1 

.  r.2  -1 

yo 

"  Allyo 

+  A12xo 

+  hlu0 

1 

2  -1 

.2  1 

,  1  -1 

X1 

"  ^1*0 

+  A22xo 

+  ®21uo 

-1 

-1-1 

y0 

-  D  uo 

-1 

-1 

yo 

"  V 

(7.2.24) 


Now  Definition  11  says  that  is  the  set  of  initial  conditions  for 

(7.2.24)  for  which  there  exists  a  control  u ^  such  that  v^"0.  The 

o  •  o 

structure  of  the  equations  in  (7.2.24)  immediately  imply  u^  •  0  since 
is  nonsingular  and 


(7.2.25) 


which  was  the  result  obtained  in  (7.2.23).  Substitution  of  these  results 
into  (7.2.24)  yields 


-1  21  2  -1 

y  1  "  A12xo  +  *12% 

1  2  1  ,  1  -1 

X1  A22xo  +  B22Uo 


(7.2.26) 
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As  vet  u1  is  unspecified.  Indeed,  an  arbitrary  choice  of  u^ 
o  o 

will  not  affect  However,  note  that  vj  will  appear  in  the  system 

output  in  the  next  time  unit.  Hence,  the  definition  of  will  require 

that  yj^O.  The  equations  in  (7.2.26)  can  be  thought  of  as  a  reduced 
2  12  2 

system  (A-,.,  .B,,,,  ,A.p  ,B.p)  for  which  it  is  desired  to  identify  £j.  This  is 

exactly  the  problem  which  was  solved  by  the  first  cycle  of  the  algorithm. 

The  second  cycle,  then,  is  Just  the  analysis  of  the  first  cycle  applied  to 
2  12  2 

the  subsystem  (A^0 , , Aj n . B^0 ) .  However,  since  the  ultimate  goal  of  this 
algorithm  is  to  produce  canonic  form  (7.2.1),  it  is  necessary  to  embed  the 
transformations  of  this  subsystem  in  larger  transformations  which  can  be 
applied  to  the  whole  system  such  that  the  already  defined  structure  is 
preserved.  The  necessary  calculations  for  the  next  cycle  are  given  to 
illustrate  this  procedure. 

2nd  Cycle 


Find  a  (n-m^ ) x  (n-m^ )  nonsingular  matrix  Wj  such  that 


B12W1  “  1  <b12)'  01  (7.2.27) 

3  * 

where  (B,»)  is  a  r.  X  im,  matrix  that  has  full  column  rank.  This  can  be 

1  d.  1  d, 

2  1 

embedded  in  an  input  space  transformation  u‘“WjU^  where  the  mvm  matrix 
Wj  is  given  by 


W 


1 


(7. 2. 28'! 


Applying  (7.2.28)  to  (7.2.20)  gives 

“2  (B?~ ) '  0 


B 


b:w 


1 


“11 

LBn 


12 J 

H3  K3 

R22  B23J 


(7.2.291 
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2  2 

Note  that  D  -  D  as  can  be  seen  directly  from  (7.2,28)  and  (7.2,20), 
Next  find  a  mj  xnij  non-singular  matrix  such  that 


i1<BJ2V  - 


(7.2,30) 


— 4^ 

where  is  a  non-singular  m->  matrix. 

In  terms  of  the  system  (7.2,20),  (7.2.30)  can  be  implemented  by 
applying  the  similarity  transformation  z^-V^z^  where  the  nxn  non-singular 


matrix  is  defined  as 


z2  -  V  r1 
i  li 


S1 

1 

o 

yi 

m 

Al 

yi 

0 

I 

X3 

x3 

— 

n-\J 

i 

^1 

_l 

(7.2,31) 


Calculate  A^V^V*1,  B^VjB3,  and  C3«C2V^L. 
now  have  the  form 


3  4 

The  matrices  A  and  B 


S11 

A12 

-4 

B12 

0 

*11 

;  b4  - 

*ii 

0 

0 

_A21 

4_ 

l 

B2l 

»3 

B22 

S23_ 

which  applies  to  C  , 

also 

applies 

to  A3 

(7.2.32) 


If  A3,,  does 

have  full  row  rank,  the  necessary  adjustments  should  be  made).  Define 
-3 

the  row  rank  of  to  be  r>* 

From  (7.2,32),  the  feedback  matrix  can  now  be  calculated 


(following  (7.2,10))  as 
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—A.  - 1^3 

<h2)  \2 


mixri 

0 

0 


0 

K, 


(7.2,33) 


where  is  a  m^  x(n  -r^  matrix  and  is  a  mxn  matrix.  (Recall  that 
3  4 

A12  and  B12  3re  t*le  "C"  and  "D"  n'atr^ces  of  the  subsystem  (7.2,26)  to 
which  the  2nd  cycle  of  the  algorithm  is  being  applied).  Calculate 

A^=A^  +  B^K^  and  note  that  ■=  +  D^K^. 

4 

Now  the  matrix  A  has  the  form 


73 


A  J 

11 

0 

73 

73 

A11 

A12 

3 

.  4 

A21 

A22 

;  that 

II 

CM 

4  4 

_3  ; 

■  , 

c  ,D 

(7.2,34) 


If  4^2  ~  0>  then  it  follows  that  £2=£  ^=£  by  the  same  argument  that  gave 


algorithm  terminates. 


~3 

If  A12  4  0,  then  proceed  as  above  to  identify  £  .  Define 


o 

1 

0 

_ 

k\ 

3 

L2 

°"1 

|  I 

(7.2,35) 


where  is  a  nxn  non-singular  matrix.  The  comments  in  Remark  2  show 
that  this  is  always  possible.  Form  the  similarity  transformation 


1 


l 

I 

l 

[ 

i: 


o 


o 


1 0 

lo 


0 

0 

D 

0 


U 


T  z~ 
llZL 
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~2 


(7.2,36) 


—2  ~3  1  54-154  3 

where  Calculate  A  “TjA  ,  B  "T^B  ,  and  note  that  C  and 

2 

D  are  unchanged  by  this  transformation.  The  system  now  has  the  form 


*la 

yi+l 

4 

0 

0 

B4 

B11 

I4 

B12 

0 

*lb 

yi+l 

4 

42 

0 

Bu 

0 
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-2 

yi+l 

m 

4 

4 

A23 

4 

B22 

B23 

2 

Xi+1 

A31 

4 

A33 

4 

4 

4 

-1 

yi 

0 

0 

0 

D1 

0 

0 

-1 

yi  _ 

0 

0 

0 

0 

0 

— 

... 

~2 

yi 

2 

xi 

-1 

ui 

-2 

ui 

-2 

ui 


(7.2,37) 


By  inspection 


£2  " 


(7.2,38) 


The  algorithm  now  enters  the  third  cycle  with  the  new  subsystem  (Aj3>Bj3» 

A5  B5  ) 

A23’  23' * 

Clearly,  this  algorithm  must  terminate  after  a  finite  number  of 
steps  (at  most  n) .  Then  the  system  matrices  have  the  following  form: 


I 

I 

I 

I 

1 

I 

I 

1 

I 

1 

1 

I 

I 

I 
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Bj,J+l  “ 

where  B ^  is  a  mj+^xnij+^  matrix.  This  is  exactly  the  canonic  form 

desired.  (For  example,  by  comparing  (7.2,1)  to  (7.2,39)  it  is  seen  that 

A4‘Al+l,t+L  “d  B4'Bi+l,i+2>-  ° 

At  this  point  it  is  useful  to  discuss  the  relationship  of 

Algorithm  II  to  the  procedure  for  obtaining  the  Generalized  Hessenberg 
Representation  (GHR)  [52].  The  relationship  can  be  seen  from  the  trans¬ 
formation  matrices  T^,  i*0,l,..  defined  in  Algorithm  II  (7.2,18).  Note 
that  the  matrices  T^  used  in  Algorithm  II  have  the  same  form  as  the 
transformation  matrices  used  in  obtaining  the  GHR  of  a  system  except  that 
the  matrices  T^  are  formed  using  only  a  block  of  the  "C"  matrix,  not  the 
entire  "C"  matrix  as  in  obtaining  the  GHR.  Furthermore,  that  portion  of 
the  "C"  matrix  that  is  used  in  T^  is  exactly  the  portion  which  cannot  be 
cancelled  using  state  feedback.  Hence,  the  GHR  of  the  system  can  be  ob¬ 
tained  by  using  Algorithm  II  when  it  is  assumed  that  B  ■  0  and  D  ■  0.  The 
unobservable  subspace  can  be  identified  directly  from  the  GHR  of  the  system. 
On  the  other  hand,  Algorithm  II  identifies  the  largest  subspace  which  can 
be  made  unobservable  using  state  feedback. 

Note  that  by  duality,  the  transformations  that  give  the  GHR  can 


B 


j.J+l 

0 


(7.2,41) 


be  used  to  identify  the  unreachable  subsystem  of  a  given  system.  This 
motivates  the  following  extension  of  Algorithm  II  which  further  reduces 
(7.2,37).  By  Corollary  1,  canonic  form  (7.2,1)  can  be  reduced  to 
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(7.2,42) 


This  is  just  the  decomposition  of  (A^,B^)  in  (7.2,1)  into  its  reachable  and 
unreachable  subsystems.  The  following  algorithm,  based  on  the  transforma¬ 
tions  to  obtain  the  GHR,  is  proposed  to  accomplish  this  decomposition.  The 
algorithm  is  a  sequence  of  similarity  transformations  applied  to  the  sub¬ 
system  (A^,B^).  As  in  Algorithm  II,  these  similarity  transformations  can 
be  embedded  in  a  larger  similarity  transformation  which  can  be  applied  to 
the  full  system  (7.2,42).  For  notational  simplicity  this  will  be  left  to 
the  reader. 

Algorithm  II  (continued) 

Let  A -A,  and  B  ■  B.  and  consider 
4  4 


wi+l  “  Aw^  +  Bw^ 


Define  a  transformation  such  that  where 


B1  ° 
B2  1 


(7.2,43) 


(7.2,44) 


where  is  a  square  non-singular  matrix.  This  implies  that  B^  must  be 
non-singular.  If  B  has  full  column  rank,  this  can  be  achieved  by  a  permu¬ 
tation  of  the  state  variables.  If  B  doesn't  have  fullrank,  there  exist 
linearly  dependent  controls  with  respect  to  the  full  B  and  D  matrices  in 
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(7.2,42).  There  is  no  loss  of  generality  to  exclude  this  case.  Hence,  a 
non-singular  always  exists. 

Now  consider 


w 


i+1 


T  T 

'A  Wi*  yi  “B  Wi 


(7.2,45) 


The  first  step  in  finding  the  GHR  for  (7.2,45)  is  the  application  of  the 

1  T 

similarity  transformation  w.  *  R^  w^.  The  procedure  for  obtaining  the  GHR 

T 

of  (7.2,45)  will  generate  a  sequence  of  transformations  R^  and  transform 
(7.2,45)  to  explicitly  display  its  unobservable  subsystem.  By  duality, 
the  transformations  w^  ^ » R^w^  will  transform  (7.2,43)  so  as  to  explicitly 
display  its  unreachable  subsystem.  The  system  (7.2,43)  after  transformation 
will  have  the  form 


A11  A.l2  .......  .  AL>i 

1 

M 

_ l 

A21  A22 

0 

• 

.  B  ■ 

CM 

CO 

< 

O  • 

•  •  • 

•  •  • 

_ 0  *  *  •  *  0  Ai,i-1  Ai,i 

- 1 

°  i 

(7.2,46) 


If  A.  .  ,  “0,  then  A  .  represents  an  unreachable  subsystem.  If  A.  .  +  0, 
1 )  i  1)1  1)1*1 

then  the  pair  (A,B)  is  reachable.  In  any  case,  (7.2,46)  combined  with 
(7.2,39)  gives  canonic  form  (7.2,42).  (For  example,  if  A.  .  ,  *  0,  then 

1 ,  L-  1 


Algorithm  II  is  closely  related  to  three  other  algorithms  which 
have  appeared  in  the  literature:  Algorithm  I  [26]  (Section  5.2),  Silverman's 

•ft 

Structure  Algorithm  [34],  and  Wonham's  algorithm  for  calculating  £  [12]. 
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The  relationship  between  Algorithm  I  and  Algorithm  II  was  detailed  in 
the  development  of  the  first  cycle  of  Algorithm  II  and  will  not  be  dis¬ 
cussed  further. 

Using  the  notation  of  Algorithm  I,  the  Structure  Algorithm  can 
be  described  as  follows: 

Structure  Algorithm: 

Step  0:  Set  i  »  0,  AQ  =*  0,  FQ=D,  and  GQ  =*  C 

Step  1:  Determine  any  nonsingular  S^  such  that 


S. 

i 


j  A  .B 
1 


A.  A 

1 


Step  2: 


where  F...  has  full  row  rank. 
r+1 

Set  i * i+1  and  go  to  Step  1. 


In  [26]  the  following  relationship  between  Algorithm  I  and  the  Structure 
Algorithm  is  proved: 


M 


i 


A 


i 


(7.2,47) 


This  implicitly  establishes  the  connection  between  Algorithm  II  and  the 
Structure  Algorithm.  Actually  a  closer  .^nnection  can  be  established 
but  the  details  will  not  be  given  here.  Suffice  it  to  notice  that  A^*C 
(7.2,6)  and  that  A2  is  related  to 


(7.2,34)  through  similarity  transformations.  Also  note  that  when  the 


Structure  Algorithm  is  applied  to  a  system  (A,B,C,D),  at  the  i-th  step 


Hence,  (7.2,49)  is  not  equivalent  to  the  original  system  in  the  sence  that 


9)  can't  be  obtained  from  (A,B,C,D)  by  change  of  basis  in  input,  out 


put,  and  state  spaces  or  by  the  application  of  state  feedback.  In  this 


respect.  Algorithm  II  and  the  Structure  Algorithm  differ 


Algorithm  I  and  Algorithm  II  above  generate  a  sequence  of  sub 


It  can  be  shown  [26]  that  the  subspaces  X.  satisfy  the  follow 


(The  proof  is  similar  to  the  proof  of  the  following  theorem).  This  algorithm 


is  a  generalization  of  Wonham's  algorithm  for  the  calculation  of  X  .  The 


following  theorem  relates  (7.2,50)  to  Algorithm  II 


'heorem  2  6 


The  subspaces  X.  satisfy  the  algorithm 


for  appropriate  choices  of  L. .  These  matrices  can  be  obtained  from  the 


feedback  matrices  K.  calculated  in  Algorithm  II 


Proof:  Algorithm  I  (5.2,15)  shows  that 


Let  u  * L. x  where  L 


Equation  (7.2,51)  can  be  rewritten  as 


Condition  (2)  can  be  written  as 


for  some  z  satisfying  M  z  =  0.  Since  x  also  satisfies  condition  (1),  the 
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first  part  of  the  theorem  follows  for  V*"  =  ^[M^  =  £^. 

Algorithm  II  provides  for  the  explicit  representation  of  the 
subspaces  £^.  It  is  easily  seen  that  the  matrices  can  be  obtained 
from  the  feedback  matrices  by  a  straight  forward  calculation.  For 
example,  is  calculated  as 

L^W^+WqKqTq1^1  (7.2,54) 

□ 

* 

Two  other  methods  for  the  computation  of  £  have  appeared  recently. 
* 

In  [53] ,  £  is  computed  by  finding  compatible  sets  of  eigenvectors  in 

where  P(s)  is  the  system  matrix  and  is  the  eigenvalue  which 
corresponds  to  the  eigenvector  being  computed.  In  [46]  a  variation  of 
Algorithm  II  appears  as  a  numerical  procedure  applied  to  the  system  matrix. 
The  theory  is  drawn  from  Kronecker's  theory  on  the  structure  of  pencils  of 
matrices.  The  presentation  there  is  for  system  with  D  =  0  and  doesn't  include 
the  geometric  interpretation  given  here. 

7.3.  Calculating  Invariant  Zeros 

Algorithm  II  in  Section  7.2  provides  a  method  for  calculating 
canonic  form  (7.2,42).  Once  this  canonic  form  is  available,  the  invariant 
zeros  (Definition  5)  can  easily  be  calculated  from  it  by  using  Theorem  14. 

The  invariant  zeros  are  simply  the  eigenvalues  of  the  submatrix  A^.  The 
invariant  zeros  can  also  be  easily  computed  from  canonic  form  (7.2,1).  In 
this  case,  they  are  the  eigenvalues  that  are  associated  with  the  unreachable 
modes  of  (A^,B^).  This  can  be  seen  from  (5.2,31);  it  is  also  clear  from  the 
proof  of  Theorem  14. 


15  6 


Using  Algorithm  II,  the  invariant  zeros  can  be  computed  for  any 

system  (A,3,C,D).  In  particular,  the  invariant  zeros  of  a  degenerate 

system  can  be  calculated  via  Algorithm  II.  Algorithm  II  differs  in  this 

respect  from  other  computational  methods  such  as  the  Generalized  Eigenvalue 

Method  or  the  Davison-Wang  method.  (The  latter  compute  zeros  defined  by 

Definition  6).  Of  course,  any  procedure  whose  theoretical  basi3  rests  on 
* 

£  (such  as  the  method  by  inverse  systems)  will  compute  the  same  zeros  as 
does  Algorithm  II.  Algorithm  II,  however,  places  no  restrictions  on  the 
matrices  (A,B,C,D)  and,  therefore,  is  more  general  than  certain  other 
methods . 

The  procedure  for  finding  invariant  zeros  is  summarized  in  Fig.  1 
the  following  flow  chart. 


STbP  1 


Fig.  1.  Calculating  Invariant  Zeros  Via 
Algorithm  II. 

Algorithm  II  is  condensed  and  summarized  in  the  following  step  by  step 
procedure. 


I 

I 

I 

I 

I 

I 

I 

I 

l 

[ 


1 

E 

E 

E 

E 

l 

l 

I 


Step  1 
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Given  an  arbitrary  system  (A,B,C,D), 


Define  Oq-0,  60“0,  and 


C  D 


A  B 


Find  non-singular  input  and  output  space  trans format  Ions ,  and  S( 
tlvely,  such  that 


P' 

r0 


S0  ° 


0  I 


ln  0 


0  W, 


"c1  Dl  0 

cl  0  0 


0  0  0 
a  »; 


m0  m  "  m0 


r  -  r„ 


—  1  **1 
where  D  is  a  square  non-singular  matrix  and  C  has  full  row  rank. 


£to£j. 


Calculate  Kq  *  -(D1)’^  and  form  the  mxn  matrix 


K0" 


If  r^-0,  !"Jt  To"la’  Otherwise  £orm  thc  m*trix 


c~l 


I. 


n-r 


respec- 
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Calculate 


Step  3 


If  r0»0,  stop.  Otherwise,  set  i -1  and  go  to  Step  3. 


Find  non-singular  transformation  matrices  and  such  that 


where  Di+1  is  a  square  non-singular  matrix  and  C1+1  has  full  row  rank. 


Form  the  nxn  non-singular  matrix 


Calculate 


i+1 


_ 

— 

I 

0 

m 

0 

T_i 

p[ 

T^1  0 

Ki 


X 

X 

X 

X 

X 

X 

X 

0 

Q 

X 

Di+1 

0 

X 

I 

ri 

0 

X 

0 

0 

X 

0 

0 

X 

0 

0 

Ai+l 

A1*1 

ai+1 

Bu 

.1+1 

X 

An 

A12 

X 

B12 

X 

Ai+1 

A2I 

22 

X 

Bi+1 

B21 

Bi+1 

“22 

If  -  0,  stop.  Otherwise,  set  i * i+1  and  go  to  Step  3.  a 

The  algorithm  to  further  transform  the  system  from  (7.2,1)  to 
(7.2,42)  is  given  next.  The  notation  is  in  terms  of  a  pair  (A,B).  If  this 
pair  is  a  subsystem  of  a  larger  system  (as  above),  the  transformation  matrices 
should  be  embedded  in  a  larger  transformation  matrix  which  can  be  applied  to 
the  full  system. 

Given  a  pair  (A,B);  A  is  nxn  and  B  is  nxm. 


Step  1 

Define 


PQ  -  [B  A] 


Find  a  non-singular  input  space  transformation  such  that 
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P '  «  p 
*0  *0 


w„ 


[B*  0  A] 

xo  m"xo 


where  B'  has  full  column  rank. 

Step  2 

From  the  nxn  transformation  matrix  Rq  as 


R0  “ 


B2 


n-\. 


Calculate 


—  — 

I 

m 

. 

Ro_ 

0 

0 


A' 

A11 

A21 


A  1 
A12 


a  > 

22 


-  P, 


Set  Tq  -  \Q  and  i  -  1. 


Step 


Set 


X  X 

X 

X  X 

X 

i 

.  i 

m 

i 

i 

_x  *21 

A2_2_ 

x  B 

A 

Step  4 


Find  a  non-singular  matrix  W ^  such  that 


bV 


[Bi+1  0] 


.1+1 


where  B  has  full  row  rank.  If  X^-0,  stop.  Otherwise  form  the  nxn 


non-singular  matrix. 


Set  +  i  *  i+1  and  go  to  Step  3. 

The  following  examples  illustrate  the  use  of  Algorithm  II  to 
compute  invariant  zeros.  Note  that  these  examples  also  Illustrate  the 
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geometrical  ideas  presented  in  Section  5.2. 
Example  17 

Consider  the  system 


Calculate  the  invariant  zeros  of  this  system  by  first  finding  canonic  form 
(7.2,1)  via  Algorithm  11. 

Step  1 

Note  that  if 


then 


1 

1 


0 

0 


is  in  the  required  form.  Furthermore,  S^*^.  Then 


Note  that  rQ  =  1  and  mQ  =  1. 
Step  2 


Since  r„  ^  0,  form  the  transformation  matrix 


Identify 


Cycle  2 
Step  3 


Step  4 

Step  5 


Since  r 


[-10  2]  ;  b: 


Since  D 


2 


has  full  row  rank, 


=  P 


2’ 


-  -3  -1—3 

k2  =  -(dj)  V 


*  [5  -L] 


K 


2 


0 

0 


0 

5 


0 

1 


o,  T0 


Then 


L68 


Find  Che  invariant  zeros  based  on  the  following  pairs  of  inputs  and  outputs; 


a)  ui  and  y 


b)  ui  and  y^ 


where  a,  “ 


c)  u.  and  y.  where  y. 


a)  u .  and  y^ 


Steo  1 


Define 


0  0 


0  0 


1  0 


1  2 


1  1 


Since  D»0,  the  algorithm  is  simplified.  It  is  seen  chat  w  *  T  c^-r 

0  *,2  ’  ^ » 


rn  *  2  ;  and  a.  ■  0. 


Step  2 


Again,  D*0  implies  Kn“0.  Form  the  transformation  matrix  T 


0 
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and 


0  0  6  -2 
0  0-42 


Since  r  *  0,  T,  » I  .  Then 
1  1  n 


x 

0 

4 


D2 


J21 


The  algorithm  terminates  as  ^  = 0.  By  writing  out  P?  and  identifying  the 
original  system  matrices,  it  follows  that 


i+1 


x.  + 
i 


3  1 

0  -1 


1  1 
2  2 


1  0 
0  1 


0  0 
0  0 


This  is  clearly  in  canonic  form  (7.2,1).  From  P,, ,  it  is  seen  that  Br,  »  0 

★ 

which  implies  that  >1  =0.  By  Theorem  14,  the  invariant  zeros  are  the 
eigenvalues  of 
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-4 

-6 


1 

l 


These  are  easily  calculated  as  z^  •  -1  and  z,  »-2. 
b)  u|  and  yt 

Since  Che  C  matrix  is  unchanged  from  part  (a),  is  the  same 
except  Che  last  column  is  deleted: 


—  ■■ 

— 

1  0 

0  0 

0 

0  1 

0  0 

0 

ft 

*-* 

1 

-14  4 

3 

O 

l 

-4  2 

0 

t 

O 

-6  1 

1 

_0  -3 

-10  1 

2_ 

Step  3 

Identify 


Step  4 

Note  chat  is  in  the  required  form,  so  that  S.  "I,,  W  "I,,  and  V  -  I  . 

12  11  In 


Hence  •  P^. 


Immediately  identify 


c  -  [  -14  4]  ;  D  -  [3] 


C"  -  [-4  2] 


X 

_ X  _ 

X 

X 

X 

X 

X 

— 

X 

X 

0 

0 

3 

X 

0 

0 

D2 

X 

1 

0 

0 

X 

X1 

Ax 

0 

0 

X 

-1 

0 

0 

X 

2 

A12 

B11 

X 

0 

-2 

2_ 

X 

4 

a22 

4 

Cycle  2 
Step  3 

Identify 


C2  -  [0]  ;  D2  -  [0] 
A2  -  [-2]  ;  B2  -  [2] 


0 

u 

u 

0 


) 


Step 


Since  D 


2  —3 

is  in  the  required  form,  it  follows  that  C“»C  *>[0] 


and  ■  0. 

Step  3 

Clearly  ^  ■  0  and  *  0  implies  T,,  a  I  .  Hence,  the  system  is  in 


the  required  form. 
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1  0  0 


0  10 


2  * 

From  P0  it  is  seen  Chat  *0  so  that  »l  -0.  Hence,  this  system  has  one 
invariant  zero  at  z»-2.  Notice  that  this  system  has  two  unreachable 
modes  but  only  one  of  them  is  associated  with  an  invariant  zero. 


c)  uL  and  yL 


Let  the  system  be  given  by 


x  + 

•1  1 


yL  *  [-1  0  -3  -2]xi 


where  the  first  two  state  variables  have  been  permuted  to  ensure  the  non¬ 
singularity  of  Tq  and  the  inputs  have  been  permuted  to  ensure  consistency 
of  notation  below. 


Step  1 


Step  2 


Form  Pq.  Note  that  C-C  and  rQ  -  1 


Note  that  Kq *  0  as  D  *  0.  Form  as  before.  Then 


I 

I 

I 

I 

I 


r 

L 


Since  r  *0,  T,  H  I  ,  and 
1  ’  1  n 


1. 
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0  0-42 

0  0  0  0 


1 

0 

0 

0 

0 

0 

-4 

0 

0 

o 

-1 

0 

-1 

-1 

0 

0 

0 

1 

-1 

0 

-10 

3 

1 

1 

-3 

0 

-18 

5 

2 

2 

1 

2 

11 

2 

12 


0 

,2 

22 


D2 


321 


0 

j2 

522 


i: 

r. 

o 

i 

r 

L 

L 

I 

i 

I 

I 


Since  *  0  the  algorithm  terminates.  However,  note  that  B2?  ^0  so  that 
★ 

ft  j*0.  Now  the  invariant  zeros  could  be  calculated  from  the  unreachable 
2  2 

modes  of  B22^’  However,  here  the  system  will  be  further  reduced  to 

canonic  form  (7.2,42). 

Step  1 

Identify 


1 

-1 

0 

0 

po- 

1 

0 

-10 

3 

_2 

0 

-18 

5 

Immediately,  \ 

C 

Step  2 


1  and  p0  “  P0* 


R0  ' 


2 


0  1 


Identify 


Step  3 


Then 
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CHAPTER  8 
CONCLUSION 

This  work  has  three  major  areas  of  emphasis.  The  first  is  a 
survey  of  selected  literature  on  zeros  for  linear  time -invariant  multi- 
variable  systems.  The  second  is  a  survey  of  algorithms  for  calculating 
zeros.  This  includes  the  properties  of  zeros  on  which  these  algorithms 
are  based.  The  third  area  of  emphasis  is  the  presentation  of  a  new 
algorithm  for  the  calculation  of  invariant  zeros. 

Transmission  zeros  were  defined  for  a  transfer  function  matrix 
using  the  Smith-McMillan  form.  The  motivation  for  this  definition  comes 
from  the  fact  that  these  transmission  zeros  are  frequencies  whose 
transmission  through  the  system  is  blocked;  a  generalization  of  zeros  of 
a  scalar  transfer  function.  These  transmission  zeros  have  other  properties 
(i.e.  alternative  definitions)  which  can  be  thought  of  as  generalizations 
of  properties  of  zeros  of  a  scalar  transfer  function.  For  example,  the 
transmission  zeros  are  the  poles  of  an  inverse  system  (if  it  exists). 

All  of  these  properties  are  closely  interrelated  through  the  Smith- 
McMillan  form. 

A  similar  analysis  is  carried  out  for  zeros  defined  for  state 
space  systems  by  using  the  system  matrix  and  the  Smith  form  of  the  system 
matrix.  In  this  way  invariant  zeros  are  introduced.  The  transmission 
zeros  are  contained  in  the  set  of  invariant  zeros;  however,  the  analysis 
is  complicated  by  the  appearance  of  decoupling  zeros.  The  decoupling 
zeros  turn  out  to  be  just  the  eigenvalues  associated  with  the  uncontrollable 
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and/or  unobservable  modes  of  the  system.  The  invariant  zeros,  in  general, 
contain  all  of  the  transmission  zeros  but  may  not  contain  all  of  the 
decoupling  zeros.  A  final  set  of  zeros  is  introduced,  called  system 
zeros.  This  set  contains  all  of  the  transmission  zeros  and  decoupling 
zeros.  The  invariant  zeros  are  a  (sometimes  proper)  subset  of  system  zeros. 

As  with  transmission  zeros,  invariant  zeros  have  several  properties 
which  can  be  considered  generalizations  of  properties  of  zeros  defined  for 
a  scalar  transfer  function.  In  fact,  the  motivation  for  defining  invariant 
zeros  is,  again,  to  identify  those  frequencies  whose  propagation  through 
the  system  is  blocked.  In  addition,  the  invariant  zeros  are  shown  to  be 
unaffected  by  state  feedback,  to  be  the  limiting  positions  of  the  system 
poles  under  high  gain  feedback,  and  to  be  contained  in  the  set  of  poles  of 
an  inverse  system  (when  it  exists).  System  zeros  (and  invariant  zeros)  are 
also  unaffected  by  input,  output,  and  state  space  transformations  as  well 
as  output  feedback. 

The  geometrical  properties  of  zeros  are  also  discussed  in  detail. 

It  is  shown  that  invariant  zeros  are  related  to  £*,  the  maximal  null  output 
(A, B) -invariant  subspace  and  to  ft*,  the  maximal  null  output  reachability 
subspaces.  These  two  subspaces  and  invariant  zeros  play  a  key  role  in  the 
discussion  of  inverse  systems.  They  also  provide  the  theoretical  basis 
for  a  new  algorithm  for  computing  invariant  zeros. 

It  turns  out  that  the  definitions  of  the  various  zeros  are  not 
very  convenient  for  actually  computing  them,  either  by  hand  or  by  digital 
computer.  These  properties,  then,  provide  the  basis  for  several  algorithms 
for  calculating  zeros.  Sometimes  these  properties  can  be  applied  directly. 


For  instance,  for  a  certain  class  of  systems,  the  inverse  system  can 
easily  be  computed.  High  gain  feedback  can  also  be  applied  directly.  This 
property  has  also  lead  to  a  number  of  other  algorithms  including  the  NAM 
algorithm.  Finally,  invariant  zeros  have  been  computed  using  a  generalized 
eigenvector  method  in  combination  with  a  QZ-algorithm.  This  method  has 
also  been  cast  in  the  more  general  setting  of  computing  the  Kronecker 
structure  of  a  pencil  of  matrices. 

In  Chapter  7,  a  new  algorithm  for  the  calculation  of  invariant 
zeros  is  introduced.  This  algorithm  is  a  sequence  of  transformations  on 
the  system.  Then  a  feedback  matrix  can  easily  be  identified  which  will 
place  the  system  in  a  canonic  form  such  that  the  subspaces  X*  and  A*  are 
displayed  explicitly.  From  this  canonic  form,  the  invariant  zeros  can 
easily  be  calculated.  The  fact  that  this  algorithm  actually  calculates 
X*  and  ft*  suggests  that  its  use  is  not  restricted  just  calculating  invariant 
zeros.  In  fact,  it  is  useful  in  carrying  out  the  construction  of  inverse 
systems  as  presented  in  Chapter  5.  This  will  be  explored  in  a  future  paper. 
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